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A Numerical Method of Evaluating the Velocity Potential and 
the Minimum Drag Warping of Arbitrary Supersonic Wings 


By Max ANLIKER, San Diego, Calif., USA?) 


Symbols 
a coefficient; 
b wing span; 
b, box length in streamwise direction; 
Aa root chord; 
" arbitrary function; 
t,4,k subscripts; 
m,n integers, also nondimensional coordinates and subscripts; 
p pressure; 


Ap=p,—f; pressure difference between upper and lower wing surfaces; 
q dynamic pressure; he i 


G subscript; 

w downwash ; Ae 
x, Y cartesian coordinates; “hy 
Re ee ro. cintensiohal partest dinates; | 

ey, | nondimensional cartesian coordinates; 

4 A(x, y) disturbed flow region within the forward Mach cone of the point 

i (x, 9); 

BC, lift coefficient; 

BC, drag coefficient; 

D pressure drag; * 
K conversion factor; i 
M Mach number; M 
ay generalized aerodynamic force; | ae at be 
ERS. zx velocity potential influence coefficient; th 
S wing area; abe 
EU free-stream velocity; » 
e 1) CONVAIR, A Division of General Dynamics Corporation. 2 
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Z(x, Y) mean chord line; 

a angle of attack; 

B / M2 —1; 

(age integers, also nondimensional coordinates and subscripts; 
ao 

: velocity potential; 

Ay=w,,—%, velocity potential difference between upper and lower wing sur- 

faces; 

A Lagrangian multiplier; 

E,7 cartesian coordinates; 

&1, 1 ‘nondimensional cartesian coordinates; 

P fundamental mode. 


Introduction 


For a supersonic aircraft the minimization of the pressure drag has become i 
an essential problem within the design procedure and calls for practical methods 
of solution. Various authors [1-6]?) have shown that for a prescribed lift the 
pressure drag of thin wings can effectively be reduced by an appropriate ; 
warping. For certain planforms the optimum drag to be gained by warping has 
been derived independant of any set of angle-of-attack distributions with the — 
so-called harmonic method [7-13]. However, the problem of the rigorous 
evaluation of the corresponding optimum warping still remains to be solved. 
A practical approach to this problem is to approximate the warping by a series 
of functions, e. g., a power series, in the stream- and spanwise coordinates. Then, ~ 
for a wing with fixed dimensions the drag can be expressed in terms of the | 
unknown coefficients, which, finally, are determined by minimizing the drag 
and satisfying aerodynamic and geometric constraints. r 

In the case of sonic and supersonic wings with pointed tip wings KAINER — 
has shown that satisfactory answers for the optimum warping are obtained if 
the latter is approximated by the first ten terms of a power series in the 
stream- and spanwise coordinates. He made use of analytical expressions for 
the pressure distribution which he derived earlier for such wings, assuming 
linearized potential flow and neglecting wing-body interference. : 

Comprehensive design studies now call for an extension of KAINER’s work — 
to arbitrarily-tapered wings with trailing-edge sweep. However, in view of the 
excessive labor, the generalization of his analytical expressions for the pres- 
sure distribution to such wings is not recommended. Instead, the adaptability 
of modern high speed digital computers suggests numerical ways of evaluating - 


*) Numbers in brackets refer to References, page 13. 


; 
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~ supersonic wings with straight leading and trailing edges and side edges parallel 
to the free stream direction (see Figure 1). 
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the drag in terms of the coefficients of the power Series which approximates the 
optimum warping. 

The process of solving supersonic flutter problems includes the evaluation 
of aerodynamic forces which often are computed by means of aerodynamic 
influence coefficients. These influence coefficients are defined as the pressure or 
velocity potential induced at a point of the lifting surface by its known motion. 
In general it is not possible to give closed form analytical expressions for 
aerodynamic influence coefficients and therefore numerical methods have to be 
applied. 

Recent investigations on the prediction of nonsteady aerodynamic forces on 
elastic wings have shown that the so-called ‘box method’ can easily be adapted 
to programmed digital computation and leads to satisfactory answers for the 
pressure or velocity potential distribution. In the ‘box method’ a grid of boxes 
is overlayed on the disturbed flow region and the downwash within each box 
is considered to be constant. For the evaluation of generalized aerodynamic 
forces the use of velocity potential influence coefficients instead of pressure 
influence coefficients greatly reduces the amount of computational labor since 
closer approximation is observed for the same grid. 


Planform 


Figure 1 
Typical wing planform with diaphragm. 


Accuracy investigations on a digital computer program for the evaluation 


_of oscillatory aerodynamic forces on supersonic wings by means of velocity 


potential influence coefficients lead to the idea of computing the steady state 


_ pressure drag due to the warping of the wings in the same manner. The program 


referred to is part of a digital computer routine for the flutter analysis of 


. 
: 
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In this program the disturbed flow region is overlayed with a ‘Mach box grid’, 
as proposed by L1[16] e. g., a system of rectangular boxes, the diagonals of 
which are parallel to the Mach lines. The velocity potential influence coefficients 
are based on linearized potential flow theory neglecting wing-body interference. 
The fundamental modes are expressed by the first 21 terms of a power series in 
the stream- and spanwise coordinates. 

Setting the frequency equal to zero and identifying a fundamental mode with 
the mean chord line of the warped wing and another one with the corresponding 
slope distribution, the drag due to lift is obtained in the form of a generalized 
aerodynamic force of this flutter program. If the warping is expressed by the 
first 10 terms of a power series in the stream- and spanwise coordinates, the 
drag due to each term can be computed as outlined and, finally, the minimum 
problem can be solved as described earlier. 


The advantage of using such a numerical procedure is that there is practi- | 


cally no restriction with respect to the generality of the wing planform, e. g., the 


leading and trailing edges may have the form of piecewise continuous curves. — 


With some modification, the procedure can also be used to solve the problem 
of the optimum thickness distribution of wings with arbitrary planforms, 
subsonic or supersonic leading edges, and blunt trailing edges. The author has 


succesfully used velocity potential influence coefficients for the evaluation of | 


the thickness drag in numerous cases. Verification of existing literature results 
showed excellent agreement. 


Drag Due to Lift 
Computed from Velocity Potential Influence Coefficients 


According to linearized potential flow theory the velocity potential w(x, y)_ 


at a point (x, y) of the disturbed flow region of a thin wing can be expressed 
in terms of a distribution of sources: 


1 ff wl, n) dé dy 
Daa V(~— g)? = (es (y we nye 


pls, y) = — (1) 


There, 4) plane is the flattened mean chord surface of the wing, w(&, ) the | 
downwash in the point (&, ), and A(x, y) is the disturbed flow region within the — 


forward Mach cone of the point (x, y) (see Figure 1) 


The pressure distribution (x,y) and the velocity potential w(x, y) are 


related by 
De, WOO (2) 


where U is the free stream velocity and o the density of the air. 


ee «-@ " 
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The pressure drag D of wings with zero thickness is defined as 
D = || Ap(«,¥) a(x, y) dx dy , (3) 


Ss 


where S is the wing area and AP the pressure difference between the upper and 
lower surfaces. In this equation 


0 O(A 0 
Ap = bu Bi =U ez lyy— w= Ve =2U9p a) 


and «(x, y) represents the streamwise slope distribution of the mean chord. 
Denoting the mean chord line of the wing by Z(x, y), the streamwise slope 


a(x, vy) becomes 
OZ (x, 
c(t, y) = C2) (5) 


and the downwash is given by 


OZ (x, y) 
Ox ~ (6) 


wis, y) = U a(x, y) = U 


for points on the planform. 
For points in the disturbed flow region off the wing planform, the down- 
wash is determined by the condition 


Ay = yy — 9, = 0. (7) 


Subdividing the disturbed flow region by a ‘Mach box grid’ and assuming 
the downwash within each box to be constant and equal to the downwash in 


‘its center, the velocity potential becomes 


eA ; i ; dé dy —. 8) 
v(x, 9) = — Xm | Ve By 0) : 


In equation (8) the integration region A; is the area of box 7 and w,; the 
downwash in its center. The summation has to be carried out over all the 


boxes the centers of which are within the disturbed flow region and the forward 
Mach cone of the point (x, y). 


-_ 
va 


If the following dimensionless coordinates are introduced: 


E Bn 
= 5 i= ao E= bi’ Mk b,? (9) 


“where b, is the streamwise dimension of a box and f is defined as /M 2_ ], 


a 
é 


the velocity potential difference between the upper and lower surface can be 


> 
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written in the form 


A(x, ye es 2 w(x, jae % ff 7 ile. ut. ) dé dn, y (10) 


— &,)? (v4 —- nae 


This expression has the advantage that the double integral over A; does not 
depend on the Mach number and can be written in terms of the relative and 
nondimensional coordinates 


y=n—yv and p=m—ypB, 


as indicated in Figure 2. The symbols ”, m, v, and uw are all integers; ” and 
denote the number of boxes downstream and m and wu the number of boxes in 
the spanwise direction from an arbitrarily chosen origin box (0, 0). 


Mach line | ates (»,n)or(m,p) 
(x€) or(n,v) 


Figure 2 


‘Mach box grid’ and reference systems. 


> e 


The double integral in (10) is usually defined as the velocity potential | 
influence coefficient and will be abbreviated in the following by 


: 
t 
Rix “al 76 dé, dy 11) 
i 
= $3)? (ya )* ( } 
Integrating (11), we obtain the formulae | 
| 
RY, =—1-0, | 
F+1/2 ro, 
R,=-= / sin-1 (=) dx = 2| See 1] (Sar lbeih O~ Fae 
p-1/2 ne ee be 12? 


Vol. X, 1959 A Numerical Method of Evaluating the Velocity Potential ... vi 


y+1/2 


th Rees. 2 ete T/2 mu 
Res _ = [sin 1( 2a al ) — 5 ax 
7-1/2 
y-1/2 
=—1+ ~ (- 5) bee. ie 2 7541/2 
wt 2 lt lM m7 1 ) 
¥+1/2 
2 ~ 1/2 
Rat = [sin (4") — sin- (4 mal dx 
7-1/2 ‘ (12) 
B-1/2 
S 2 (i As m7 i ao b= ap 
M4 ad i be pr ere 
p=12 
w+1/2 
x a ea | Sin-* »: 1—/V1— p? EST 
= (w+ a lu In( bb ) _ B+1/2 
p—1/2 
(0 <m<>) 
The substitution of (11) into (10) yields 
b 
AY n,m oF pth & Wyn , (13) 
», I 


where the summation has to be carried out over all the boxes within the region 
of disturbed flow and the forward Mach cone. 
Combining equations (3), (4) and (5), the pressure drag can be written in 


the form Us & 
alee SY SE ay, dy (14) 


If the slope distribution 0Z/0x is continuous and can be differentiated, the 
partial integration of (14) with respect to x leads to 


ne rae AE! Sasi Aeanaha (15) 


where (LE) refers to leading edge and (TE) to trailing edge. Along the leading 
edge Ay is zero, therefore (15) reduces to 


D=agtl | [40a 


| (16) 


In carrying out the double integration in (16) it is assumed that the velocity 
potential difference Ay as well as 0?Z/0x? within a box are constant and equal 
to the corresponding value in the center of the box. The Ay values of boxes 


Nhe 
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which are cut by the trailing edge undergo a weighing process. Hence, the 
double integral in (16) can be replaced by a double summation over the entire 


planform: 
02Z _ B A 
I Ap jar di dy =F SAV m (Sa), (17) 
S ? 
Similarily, the eee: along the trailing edge is replaced by a sum 
OZ 

Ay 22 dy ==» SAY, om ; (18) 

‘i ps aS foe Purp, ce eae 


where (0Z |0%) n.. », is the streamwise slope at the point of the trailing edge 
with the y-coordinate m 6,/B, and AVn pps is the corresponding Ay value 
obtained by linear oon With (17) and (18) the drag becomes 

WA 


D=ou{ Eee Gare es) ai: (19) 


From (19) the drag coefficient 


D D 


CoG Sap URS 


can readily be evaluated. 


The Evaluation of the ‘Minimum Drag Warping’ 


a) Basic equations 


(20) 


The unknown warping «(x, y) which minimizes the drag due to lift can be — 


approximated by a series of functions: 


= Dy A; ; (x, ¥) . (21) 


The pressure difference between the upper and lower surface induced by this 
warping is 


Ap(x, y) = Ue oe ee a; Ay, = Xa Ap; , (22) 


where Ay; is the velocity potential difference and Ap, the pressure difference 
caused by «,(x,y). With (21) and (22) the drag coefficient C p becomes a_ 
quadratic form of the unknown coefficients a,, ds, ... , pe 


a= ff [rontes)|[ Sa? AP aay 


1, mae a 


intend 2a ee 
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or 
Ch >, a; a; || «,(x, ¥) BO se dy (23) 
t=1 ve 
j-1 5 


Introducing the abbreviation 


LS pi Apt 


the expression for Cp can be rewritten as 
Cy = a Cy;; a; a; . (25) 
i 
j=1 


For subsonic leading edges, expression (25) has to be supplemented by the leading 
edge suction in terms of the coefficients aj, dy, ... , A- 

Finally, the aerodynamic and geometric constraints are also expressed in 
terms of the coefficients a,, dg, ..., a): 


Pky ty) © fades... ayaa 0's "fz (21, ds, +. 55.4, O02 (26) 


Thus, with A,, A.,..., A, as Lagrangian multipliers the minimum conditions 
are given by 


daz |Co+ 2h] =O PERO RAE. (27) 


_ Since the evaluation of the coefficients a@,, dj, ..., a, from equations (26) and 
_ (27) can be considered as a familiar and straight forward operation, the re- 
mainder of this section deals with the more interesting part of the problem, 
_ the evaluation of the Cp,,’s. 


b) Computation of the Cp,,'s 
In a flutter investigation using a modal analysis the contribution of the 
i-th fundamental mode gq, to the 7-th generalized aerodynamic force Q; is 
- defined as 
= |] Abie, ») pile») ax dy, (28) 
iS 


I 


“where Ap,(x, y) is the pressure distribution induced by the 7-th fundamental 
"mode ;. Setting the frequency equal to zero, Ap; is given by 


dy,/0x) dé dn 
Ap, = U 9 (dy) = UP 0 5, ( =) be rage 


PRN EOS NTE RE: 


10 Max ANLIKER ZAMP 


Substituting (29) in (28) and integrating by parts with respect to x, one obtains 


do; 
Oe | [49 dy — |[ 4y, AL oe (30) 
TE 3 


Comparing (16) and (30), one immediately realizes that for 
0 
9i(%, 9) = ay Vil® Y) (31) 


equation (30) yields the pressure drag of a warped wing planform, the mean 
chord of which is defined by 


Z(x, ¥) = pil®, 9) - (32) 


In the case of pointed tip wings KAINER[6] has obtained satisfactory : 
answers for the optimum warping by approximating it with the first 10 terms 
of a power series in the nondimensional stream- and spanwise coordinates 

x 


phe 
root 


y 

ee pees 

ieee soe (33) 
Encouraged by this, the mean chord is expressed accordingly as a power series 
in the same nondimensional variables and each term is defined as a fundamental — 
mode. Then, the Cp,,’s are determined from the generalized aerodynamic — 
forces by multiplying them by a conversion factor K: ' 
1 


Ke ee (34) 


Comparison with Analytical Results 


As a check of the accuracy achieved in computing drag by means of velocity © 
potential influence coefficients, the Cp,,’s of a delta wing with 60° sweep have — | 
been calculated for 3 different Mach numbers, corresponding to a subsonic, ; 
sonic and supersonic leading edge condition: 


M =1-715; M=2:0; M=2°5. 
The results for M = 2-0 are tabulated in Table 1 in the form 


CSAS tsicey =)M*—1C5, (35) 


eR te 


together with the corresponding analytical results derived by KAINER. The 
fineness of the Mach box grid used in the numerical approach is indicated in the 
maximum number of boxes (m4) along the reference chord, The accuracy 


. 
; 
. 
i 
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the numerical approach with the indicated value for 1m. Can be deserted by 
the maximum deviation A jax of the numerical Cj,, from KAINER'S analytical 


values: 
M = 17155 tore = 20, Ames = Bie 
M = 2.0, Nig Oi aap fe 


M=255, Mman = 18, Agere 2:5%. 


Since KAINER’s analytical work has been extended through a numerical inte- 
gration process{6] to drag computations on pointed-tip wings, the Chis 8 
obtained with the two independent procedures for an arrow wing and a diamond 
wing are also compared (see Figure 3). The maximum deviation in these cases 


are: 
Arrow wing: M = 20,” nan = 39, Amag = LO%: 


diamond wing: M = 2-0) mga, = 20, a4 = e- 


60° 


é Diamond wing 
Arrow wing 


Croat 


Figure 3 
Pointed-tip wings. 


It should be noted that KaingEr’s results for other than delta wings can not — 
be considered as rigorous since they are obtained through a numerical inte- 
gration process. The fact that in the present numerical setup using velocity 
potential influence coefficients 1,,g, can easily be made larger than 30 indicates 
that more accuracy can be achieved than actually is anticipated in using 
linearized potential flow theory. t 

In order to illustrate the adequacy of the approximation of the optimum 
warping by the first ten terms of a power series in the nondimensional stream-_ 
and span-wise coordinates, the near optimum drag for a prescribed lift has _ 
been computed for the wing planform shown in Figure 4. For two Mach num-_ 
bers corresponding to a sonic and supersonic leading edge condition the results. 


are compared in Table 2 with the corresponding values obtained by using the» 
harmonic method [13]. ; 


ee ee 
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Table 2 
Comparison of Harmonic and Numerical Methods 


Cp Cp Cp 

Mach number (gar) flat plate (ear) harmonic (rer) numerical 
Ne = 2:0 0-234 | 0-1835 | 0-1865 
M = 2:5 


0-2315 | 0-2169 | 0-218 


60° | 
15 CRoot 
0:25 Crnot 
Figure 4 


Swept wing with side edges. 


Concluding Remarks 


The values for Cp,, obtained with this numerical method show good 
agreement with the corresponding analytical and semi-analytical values in the 
case of delta and pointed tip wings. There is no indication that the numerical 
method should be less accurate in cases of more general planforms, e. g., plan- » 
forms with leading and trailing edges in the form of piecewise continuous 
‘curves, for which no analytical results are available for verification. 

The possibility of a check of the final answer in the form of the near optimum 
drag against the value obtained from the harmonic method makes the numerical 
procedure outlined in this paper a convincing and practical method of solving 


the optimum drag problem. 
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Zusammenfassung i § 


Es wird ein numerisches Verfahren zur Ermittlung der annahernd optimalen: 
Verwindung von diinnen Uberschallfliigeln beliebiger Form aufgezeigt. Die un- 
bekannte Verwindung wird durch eine Reihe von Funktionen der in der Strémungs- 
und Fliigelrichtung gemessenen Koordinaten approximiert. Mit Hilfe von Einfluss- 
zahlen wird die jedem der Glieder dieser Funktionsreihe entsprechende Geschwindig- 
keitspotentialverteilung berechnet, wobei die Giltigkeit der linearen Potential- 
stromungstheorie vorausgesetzt und die Wechselwirkung zwischen Rumpf und 
Fligel vernachlassigt wird. Dies erlaubt, den durch die Verwindung verursachten 
Widerstand als quadratischen Ausdruck in den unbekannten Koeffizient : 
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lich werden diese Koeffizienten bestimmt, indem der Widerstand zu einem Minimum 
gemacht wird unter Beriicksichtigung von aerodynamischen und geometrischen 
Nebenbedingungen. j 

Das Verfahren lasst sich auf Tragflachen anwenden, deren Vorder- und Hinter- 
kanten die Form von stiickweise stetigen Kurven haben. Ferner lisst es sich aus- 
dehnen auf Kombinationen von Fliigeln, Héhenflossen und Flugzeuge des Enten- 
typs, sofern diese annahernd in einer Ebene liegen. Ausserdem kann es ohne weiteres 
zur Ermittlung der optimalen Fliigelstarkenverteilung in Fliigeln beliebiger Form 
und Fliigelkanten in Unter- oder Uberschallstrémung sowie stumpfen Hinterkanten 
herangezogen werden. 

Fiir Fliigel mit geraden Vorder- und Hinterkanten und Seitenkanten parallel zur 
freien Strémungsrichtung ist ein Programm fiir ein Digitalrechengerat aufgestellt 
worden. In diesem Programm wird die unbekannte Verwindung durch die ersten 
zehn Glieder einer Potenzreihe in den dimensionslosen Koordinaten in Strémungs- 
und Fliigelrichtung angenahert. Analytische Ergebnisse fiir Deltafliigel sind anhand 
dieser numerischen Methode verifiziert worden, wobei gute Ubereinstimmung be- 
obachtet wurde. 


(Received: May 17, 1958.) 


Theory of Oscillation Type Viscometers IV: 
A Thick Disk’) 


By Atronso GILt AzpeiTi1a?) and GORDON FRANK NEWELL), Providence, R.I., USA 


Abstract 


The viscous drag exerted by a fluid on an oscillating disk is determined by 
‘methods similar to those used in Part III. The present paper differs from 
Part III in that the boundary layer thickness is assumed here to be small com- " 
pared with both the thickness and the radius of the disk. An extrapolation of 
the formulae derived here agree sufficiently well, however, with those derived 
for the thin disk (boundary layer thickness large compared with the thickness 
but small compared with the radius) that it is possible to make an interpolation 
between the present formulae and those of Part III that will permit the evalu- 
ation of the drag for arbitrary shape disks provided the boundary layer thick- 
ness is small compared with the radius. 
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of Scientific Research of the Air Research and Development Command, under contract No. AF 18 
(600) 1548. Reproduction in whole or in part is permitted for any purpose of the United States 
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by a grant from the Alfred P. Sloan Foundation. 

2) Brown University, now at the University of Massachusetts, Amherst, Massachusetts. 


8) Brown University. 


TN ee eee ae 


<u 


16 Atronso GILL AZPEITIA and GORDON FRANK NEWELL ZAMP 


1. Introduction 


This series of papers has dealt so far with the oscillation of a cup filled with 
a fluid [1,2]4) and a thin disk surrounded by a fluid [3]. The treatment of the 
thin disk dealt with a situation in which the thickness of the disk was small 
compared with the boundary layer thickness of the fluid which in turn was 
small compared with the radius of the disk. 

Even if one has a disk with a thickness very small compared with its radius, 
the restriction on the boundary layer thickness considerably limits the range 
of fluids for which the results of [3] can be used satisfactorily. We consider 
now another range of conditions, namely that for which the boundary layer 
thickness is small compared with both the thickness of the disk and the radius. 
We refer to this as a ‘thick’ disk to distinguish it from [3] even though the disk 
may in fact be thin from the point of view of geometry alone provided the above © 
conditions are true. A 

A significant fraction of the experimental data that has been taken using 
an oscillating disk falls into this latter range and. the main difficulty that re- 
stricts the use of the instrument as a viscometer is the lack of an adequate 
theoretical treatment of the fluid motion near the edges where the cylindrical 
walls meet the top and bottom surfaces of the disk. Part III also dealt mainly 
with an edge correction but it was based initially upon the analysis of the fluid 
motion near the edge of a disk of zero thickness (near a 180° turn) whereas the 
present paper will be based mainly on an analysis of the fluid motion near right 
angle turns. 

Although we use many of the same methods as described in Part III, the 
details are quite different. The notation and general formulation have been 
described in Part [ and also reviewed briefly in Part III, so they will not be 
repeated again here. 

In essence the problem is to evaluate the integral 


mic .07s ow | 

Di) = 25 [fe = ao, (1) 

A > 

equation (1,16), in which w is a solution of the equation | 
| 

o2w 3 dw o2w q 

af TE Ver aes oe (2) 


4) Numbers in brackets refer to References, page 34, Equation numbers and other references _ 
to earlier parts of this series will be denoted by the prefix I, II or III. For the most part, howey 4 
the present paper is dependent only on Part I which also contains a table of notation. 
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with the boundary conditions 
w=1 on thesurface A of the disk. (2a) 


The geometry and some of the notation are shown in Figure 1. It is con- 
venient to deviate slightly from our previous convention, however, and define 
the coordinate 7 to be measured from the bottom of the disk rather than from 
the plane of symmetry. The coordinates x, 7 measure positions relative to one 
edge of the disk whereas €, 7 are the usual cylindrical coordinates measured from 


Figure 1 


Geometry and notation. 


the center of the bottom surface of the disk. The radius of the disk is &, the 
thickness 2 m) and all lengths are measured in units of the boundary layer 
thickness 6 = (v/w,)'”, as before. Because of the symmetry, it suffices to con- 
sider only the quarter space shown in Figure 1b, with the condition 0w/0y = 0 | 
at the symmetry plane 7 = mp and w finite on the axis § = 0 or x = &. The 
assumption that the boundary layer thickness is small compared with both the 
thickness and the radius of the disk is equivalent to the conditions 


No->1 and Ey >1. (3) 


It is clear that the surfaces of constant w must be nearly vertical (w inde- 
pendent of 7) in a region along the vertical walls of the disk and nearly hori- 
zontal (w independent of &) along the bottom of the disk except in the relatively 
small region within a distance of order 1 from the edge where the surfaces 
must join each other smoothly. The integral in equation (1), however, extends 
over the entire surface of the disk and the contribution to D(s) from the region 
of distortion near the edge will be relatively small compared with the contribu- 
tion from the remaining part of the surface. 

a 

ZAMP X/2 
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One approximation we can make here is to neglect the edge entirely and 
consider a function w, defined in the various regions of Figure 1 by 


eS Ny ay, 
EK, (& Vs) ’ 
eer ie an Rae 


exp (V/s) .. in III | 


with K, denoting the modified Bessel function of order 7. This function does 
satisfy equations (2) and (2a) except for the fact that it is discontinuous across — 
the boundaries A’, A” of region II. Substitution of this into equation (1) gives 
for D(s) the approximation . 
3/2 K, (Vs) 
D,(s) =m #? b + 4 me ay (5) | 
in which ‘ 


man oO ET nes (6) 
An asymptotic expansion of equation (5) for & >> 1 and |s| ~ 1 gives 


3 


Dyls) =m 9 | + Se (1 + 


3 
. te +-\L, Gall 


; 
| I 
z 
an expression which has been used previously by KESTIN and WANG [4] and] 
by Kerstin and Moszynski [5] to evaluate experimental data. 

Although equation (4) represents a very good approximation to w far from 
the edge, it is very poor in and near region IT. To 1 improve this situation, it is é 
convenient to modify equation (2) slightly. Everywhere that w itself has a 
value appreciably different from zero, the first derivative term in equation (2) ; 
1s smaller than the term sw, or possibly one of the other terms of equation (2), : 
by a factor of order &>' or less. We shall, therefore, consider here the simpler 
equation 

oe op Be (7), 


subject to the same boundary conditions, equation (2a). ; 
If we should also neglect the edge effects in a manner analogous to the 
above, we would obtain an even simpler approximation than equation (4) 


namely 

| exp (Vs x) indy 
) crawl he j (8 
exp (V/s) in IIL. 


Wy = 
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This satisfies equation (7) except for the discontinuities on A’, A”, and leads 
to the expression 


Day(s) = m9 (1+ 4-1). (9) 


Equations (8) and (9) represent the leading terms in the asymptotic expansion 
of equations (4) and (5). 

Rather than trying to determine a continuous solution of equation (7) di- 
rectly, we shall seek instead a second solution of equation (7) 


W, = W— Wy, (10) 


which will give a continuous solution when added to w,. The advantage of this 
is that w, is expected to decay exponentially in all directions. We may, there- 
fore, disregard, for the moment at least, the boundary conditions at =p = i 
and x = &, > 1 by asking instead that w, > 0 for 7 or x oo. From the bound- 
ary conditions on w and the properties of w»), we obtain the following conditions 
for w.: 


@(4;%) = 0 on AL, (11a) 

w9(0,, 4) — wa(0_, 7) =—exp(/sm) on A", (41) 

: | 4 is continuous across A” , (11c) 
W(x, 4) = We(n, x) . (11d) 


Conditions (11b), (11c) and (11d) guarantee that w and 0w/0x are continuous 
across A” and A’. 

An evaluation of w, in section 2 will immediately give a first edge correction 
to D,(s), (correct. to the lowest power of &,'), which will be improved in sec- 
tion 3 using variational methods. 


‘| 2. Solution at a Right Angle Turn 

g The technique for solving equations (7) and (11) is based in part on the 
Wiener-Hopf method used previously in section III 2., but the geometry of the 
present problem is much more difficult to handle. Although we are interested 
in the solution of equations (7) and (11) for complex values of s and particu- 
larly for s ~ 1, it is convenient to think of s as a real positive number. We can 
always obtain formulae for complex values of s by suitable analytic continua- 
tion of formulae obtained for real positive s. 
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We are interested in a solution w, defined only for 7 <0 or x SO but it 
is desirable to have w, formally defined everywhere. We, therefore, choose 
w, = 0 for x > Oand 7 > 0. Actually we can simultaneously obtain other use- 
ful solutions of equation (7) which are defined for 7 > 0 and x > 0. For example, 
we could solve the corresponding problem for the fluid motion inside a cup near 
the edge, a problem which was considered previously in section II 3. There are, 
however, easier ways of finding solutions of equation (7) for the inside of a 
right angle turn than the method to be described here and there is not much 
to be gained by considering the two problems together. 

Since w, is bounded and goes to zero exponentially for x or 7 > +00, the 
double Fourier transform 


+00 +00 


A(p, q) = [dx | dn O**° w4(x, 1) (12) 


=O, = 'CO) 


is assumed to exist for all real # and g and have an inversion 


+0cOo +00 


wale, 9) = (22=)-* [dp [dg e~'***# AUB, g) (13) 


—0o 


If we take the double Fourier transform of equation (7) and integrate by 
parts, we see that A(p, q) satisfies the equation 


(p2 + gq? +s) A(p, 9) Sai ean {|e —ip wa] ae a = 2p w|_t 


—~co 


+00 
ae Baal fe uy ee es 
+f Woe at 2 n=0_ mM oa 7 
Using the boundary conditions, equation (11), this can be simplified to ‘ 


(p+ @ +s) Alb, 9) = f(g) + 0) —p(g—iys)'—q(p—iys)*, aay 


in which 


jt) = face [2] as) 


Although A(p, q) and {(p) are both unknown, equation (14) at least ae : 
restrictions on the functional form of A(p, q). Furthermore we see from equa- | 


tion (15) that is defined and i : : | 
Pe f(P) ined and analytic for 0 S Im 4, ie., in the upper ae | 


‘ 
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We now turn our attention to the properties of w. »(x, 0_), which is known to 
vanish for x > 0. 


: +00 0 
h(p) = fers W(x, 0_) dx = fees W(x, O_) dx (16) 


is therefore defined and analytic for Im # < 0, in the lower half p-plane. 
We can obtain a relation between the unknown functions h(p) and f(p) b 
making a partial inversion of A(p, g), namely 


ip) = Jim oh fem" Alp, a) a (17) 


The integral can be evaluated by substituting A(f,q) from equation (14). 
Since /(7) is analytic for Im g = 0, the path of integration can be closed in the 
usual way by a semi-circle in the upper half plane. The only singularities en- 
closed by the contour are simple poles at g = 7 (s + £2)" and g=iVs. We thus 
obtain 


2 2)1/2 p fu TA 2)1/2 i[ys — (s + p2)??] Vs—tp 
(s + p?)"" h(p) = fp) + f(s + py”) + p Rapremener uct 


This equation, as derived, is valid for real values of f, but the analytic con- 
tinuations of f() and h(p) must satisfy equation (18) everywhere. Furthermore 
{(f) must be analytic in the upper half plane and h(f) analytic in the lower half 
plane. This information will permit us to determine the singularities of h(p) in 
the upper half plane and the singularities of (p) in the lower half plane. 

_ The presence of f(7 (s + #?)1”) as well as f(p) in equation (18) will cause 
tather complicated branch points at p = +7 Vs. If we were to investigate the 
properties of /(f) and /(p) in the p-space, we would sooner or later be forced to 
perform an analytic continuation up infinitely many Riemann sheets. We can 
avoid this unpleasant necessity and at the same time simplify the algebra by 
making a multiple-valued transformation of the independent variable from p 
to 6. We let 

p=iVssin 0, (19) 


which gives 
| (s + p2)! = |/s cos 6 = J's sin (F 3 0). (19a) 


-- 


__ This transformation produces a one-to-one mapping of the entire p-plane 


a 


into a vertical strip in the complex 6-space with — 2/2 < Re 6 S 2/2. It also 


Prey 
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maps the p-plane into every strip (2m — 1) /2 <Re 9 < (2m-+1) 2/2 for 
every integer m. 
For convenience in notation we write 


(0) = 4 (és sin8) = 106) nO) =h (is sin0) = h(p) (20) 


so that'equation (18) takes the form 


2\/s cos h{6] = f(0] + f [F 0| + sec + cosecO + tan6 — cot@. (21) 


Since f(p) is analytic in the upper half plane, /[0] is analytic in the image 
region 0 < Re 06 < 27/2, regions I, I’ of Figure 2. In particular, f(p) is analytic 


; 


p-plane 8-plane 
he: Figure 2 
“4 multi- calues mapping. equation (19). The p-plane is on the left, the 0- eee on the right. : 
at p=17Vs and since sin(z/2 + 0) = sin(z/2 — 6), equation (20) implies that 
the analytic continuation of /[@] to the right of regions I, I’ must satisfy the 
condition 


[5 +0] =1[5-9]. hehe 
By similar arguments, we conclude that A[@] is analytic in the strif 


— 2/2 < Re 0 S0, regions II, II’ and that the analytic continuation of (0) 
to the left of these regions must satisly the condition 


nl +6] =4/-F—o]. ) 


Equations (22) and (23) imply that /[6] is analytic at least in regions I, I’, II | 
and III’ whereas h[6] is analytic at least in regions II, II’, IV and IV’. : 

Equation (21) now shows that with A{@] analytic in II, II’, IV and IV’, the 
analytic continuation of f[6] into these regions has only one singularity, | 
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simple, pole on the boundary at 0 =.—z. By similar arguments we conclude 
that h[@] has a second order pole at § = 2/2. We could now proceed to find the 
singularities of /[@] in the next two strips to the right of III, II’ and h[6] in 
the strips to the left of IV, IV’ using equations (22) and (23). Equation (21) 
would then determine the singularities of [6] and /[0] respectively in these 
strips. By iterating this procedure, which is equivalent to the analytic continu- 
ation up successive Riemann sheets of the #-plane, we could determine the 
singularities of /[0] and A[6] in every strip. 

A systematic way of doing this continuation is by means of a difference 
equation. If we solve equation (21) for 4[6] and substitute into equation (23), 
we obtain an equation involving only the function /[6]. By using equation (22) 
also, this equation can be put into the form 


fian+o)+f[*F +6] +/ E + 6] + #(0] = —2 cosecd . (24) 
The most general solution of this a difference equation is 


(0) = A(8) 6) cos” + B(6) sin” + (0 )+6 D(6)| cos (2 6) — cosec#, (25) 


Gin which A(8), B(6), C(8) and D(6) a unspecified periodic functions of period z/2. 


If we allow A(6) to have a singularity at some point 6 = 6), then it will also 
have identical singularities at 6) + 7 2/2 for every integer value of 7 including 
four points in the strip - a < Re 0 <7. Since /[6] is analytic everywhere i in 
this strip, all of these singularities must be cancelled by singularities of B(6), C (6), 
D(6) and/or cosec§. The coefficients of 4(6), B(0), C(0) and D(@) are linearly 
independent functions, however, and except for 9) = 0 or 2/4, the matrix of 


the four coefficients evaluated at each of the four points in question is non- 


singular. We conclude from this that it is impossible to choose A (8), B(@), C(A) 
and D(@) in such a way that their singularities cancel in equation (25) at all 


four points. The only places that these functions may have singularities are at 


the points 0 = j 2/2 and 6 = 2/4 +j/2. The former singularities are necessary 
to cancel the poles of cosec@ in equation (25). The latter may occur only in C (8) 


and D(6) and then only because the coefficient cos(2) vanishes at these points. 


ay,** 
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The most general form for the A(6) etc. consistent with the above descrip- 


tion is 
A(6) = A cot (2 6) + A’(6), 
= (ON 
B(0) = B cot (2 6) + B’(6) oe 
C(6) = C, cot (2 6) + C, tan(2 8) + C’(9) , 


D(6) = Dy, cot (2 6) + Dz tan (2 6) + D’(6) , 


et ee 
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in which A, B, C,, C., D,, and D, are constants and A’(9), B’(9), C’(0) and D’(6) 
are entire periodic functions of period z/2. The additional requirement that the 
poles of the right side of equation (25) must have zero total residue at each of 
the four points 9 = —x/2, 0, +2/2, and m gives us four simultaneous linear 
equations for the four constants A, B, C,, and D,, the solution of which is 

4 4/3 2 


4 
B => — —— CF => and D == ira (26a) 


Le erg 9 ’ 3 


Equation (25) now takes the form 


(10) = = cot (2 0) [2 cos oe % sin ae + (1 — a) cos (2 0) 
20 


— cosec8 + (Cy + 6 Ds) sin(2 0) + A’(6) cos" + BY(6) sin ¢ (27) 
+ [C’(6) + 6 D'(6)] cos (2 8) . 


We finally determine /[@] uniquely by forcing it to vanish as Im 6 > +00. ; 
The necessity of this can be seen by the following arguments. We conclude from — 
equations (15) and (16) along with a few obvious restrictions on the behavior 
of w, near x = 0, 7 = 0, that /(p) > 0 for Im > + coand that p h(p) > — 1. 
for Im p> — oo. This means that f[6]>0 for Im 6 > + co with 
0 < Re 0 < x/2 and because of equation (22) also with 2/2 < Re 0S a. 
Similarly for Im 6 + + oo, sin# h[0] > 0 for — a < Re 0 S 0. These conditions 
along with equation (21) guarantee that {[0] > 0 also for Im 9 > + oo and 
—a< Re #0. Finally equation (24) guarantees the same for all values of 
Re 6. 

It is easy to select the unknown terms of equation (27) in such a way as to | 
cancel the known infinite terms for 9 = +700. By choosing B’(@) = 4/3, one 
can eliminate the singularity due to the first term in equation (27). By choosing 
A'(6) = 4/(3 V3), C, = 2/3 and D, = —4/(3) and finally C’(6) = D'(@) = 0, - 
one successively cancels the singularities of the remaining terms. By using these _ 


values of A’(9) etc. and some simple trigonometric identities, we obtain as a 
possible solution ; 


sat 2 4 40 0 - 
/(0] 3 sin (20) Fz cos| as g) + 1-47 — 3 cos6]. (28) 


/ = ; 

If [0] were any other solution, then /’ [0] — [6] would be of the same form. 

as the unknown part of equation (27), i.e., | 
am 


f'(0] — f{6] = (Cz + 0 Dt) sin (26) + A"(6) cos7,2 + B16) sin 


+ [C"(6) + 6 D"(6)] cos(2 6) , 
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with A"(6) etc. entire and periodic with period 7/2. Furthermore /’(0)] —7(6] 
must vanish for Im 9 > + co, Re 0 arbitrary. 

The coefficients of 4"(0) etc. in equation (29) are not periodic with period 
x/2 and if we replace 6 by 0 + 72/2 for various integer 7, we obtain new expres- 
sions involving A" (6 + 7/2) = A”(8) which must also vanish for Im 9 > + oo. 
One can easily construct in this way six linearly independent expressions in- 
volving the six unknowns Cj, Dj, A"(9) etc. all of which must vanish for 
Im 9 > + 0e. We conclude from this that C| = Dj = 0 and that A"(0), B"(6), 
C"(8) and D"(6) all must vanish individually for Im 9 > + co. The only entire 
periodic function which vanishes for Im 0 + + 0 is the null function. We 
therefore conclude that A”() etc. are all identically zero and /(6] = f’[0] isa 
unique solution. 

Having found /[6], we can now determine /(f) from equation (20), A(f,q) 
from equation (14) and finally w,(x, 7) from equation (13). Fortunately we will 
not need to evaluate w,(x, 7) explicitly because all the things we wish to cal- 
culate can be expressed directly in terms of /[0]. 

Our immediate objective is to find D(s), equation (1), which is given in 
terms of our approximate w by 


on 


No £ 
= pany esmesis| fen PS] fee 2) Pith 


D(s) ~ Dols) + 9 88s > | B [S| do 
4 


(30) 


The integrands in the above integrals decrease exponentially and so for 
~&5, 79 > 1, we replace the upper limits of integration, jy and 5, byoo. We then 
expand (1 — x/&,)3 and notice from equation (15) and (11d) that the integrals 
in question can be expressed in terms of f(f) and its derivatives at p = 0. These 
in turn can be expressed in terms of /[6] and its derivatives at 0 = 0 through 
equation (20) and evaluated explicitly from equation (28). We thus obtain 
another approximation to D(s), namely 


16 4 il 10 ' 
Ds) = Ds) + ms" 32 (2 Nae reais eed 


| The errors in this approximation to D(s) are of two types. One type consists 
of negative exponentials in 4 and & or quantities of comparable magnitude 
that arise from neglecting the boundary conditions at 7 = % and €é =). The 
other errors are a result of replacing equation (3) by equation (7). 


2 eh bea OO, 
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The complete asymptotic expansion for D(s) is expected to be of the form 


Pe are Sa teen hte [ae Stated |e 
(m 53/2) ‘gp (V/s &)? Lae pe, (V/s &)? itn (G2) 
+ exponentials in (&, o) - | | 


Equation (31), in effect, determines the coefficient a, of equation (32) but the 
approximations made above are such that it does not give the correct value of 
a, nor any of the 0,’s. Actually the 0;’s have nothing to do with the edge effect — 
and they are correctly given by equations (5) and (5a) or by replacing ‘D,(s) 
in equation (31) by D,(s). = 

By using a variational method we shall now try to find the correct a, and 
obtain an estimate of some of the exponentially small terms. vee 


3. Variational Method 


The results of the previous section can be improved by using the variational 
form for D(s) derived in section III 4, 


Dis) = 25~ fave (Sr) + Cae) Pom] OM 


To use this effectively we must construct an approximation to w which 
exactly satisfies the boundary conditions of Figure 1. The best solution that — 
we can construct from the solutions at hand is the following: 


1 


(34) 
| 7 ¢ 
This is a solution of equation (7) which satisfies the boundary conditions at 
= No. AS, equal to 1 on A and is continuous across A’ and A”. It is deficient — 
in only two respects; it does not satisfy equation (2) and there is a slight (ex- 
ponentially small in 7) discontinuity in the derivative of w across A”. { 

We can integrate equation (33) by parts to obtain 


(s) T ‘ of w on M até w | oes eee 


| W(x, 4) + We (x, 24 — y) + Wolx, y) \ in Land If, | 
w(x, n) oy: oie 


W9(x, 7) + W(X, 7) in III. 


in which the surface S must now include the two sides of the surfaces Al sa A" 
as well as the surface A. By substituting equation (34) for w, using the boundary 
conditions on wy and w, and the fact that they are solutions of equation (7) 
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we can decompose equation =) into the following terms 


= | \80 | 
D(s am. | — 2’) Wy | ow 
af alt fet é cle hak a “| 


& \3 OW 
+ fa (1 -- a = le r fant [1 — w,(x, )] 
i) 


"o 


X We (x, 2 ny ~ 1) |, + hfs dy (1 — ——)-(w ot Ws) ) + (36 


0 § 4% \4 
are (Wy + wy) + few ” (1 — ae (Wo + Ws) 


—c% -o 


ra) 
x Ox (Wo he Wg) | ? 


with the first five terms coming from the surface integral and the last two from 
the volume integral. 

By using equation (34) for w, we can expect to determine D(s) correctly to 
only one power of €, more than found in the previous section. We will, there- 
fore, discard any terms of equation (36) proportional to &>* or exp(— (— Vs &). 
This leads us to replace the upper limit of integration & in the fourth term by 
coand the (1 — x/&)) in the sixth and seventh terms by 1. We then notice that 
the first four terms give D,(s), equation (31), and that the sixth and seventh 
terms can be integrated with respect to x. The fifth term is exponentially small 
in 7% but since nothing else limits the range of 7, this term is worth keeping. 
We use equation (11d) and the inversion of equation (15) to sEpIace Ow,/0x i in 
this term and so transform equation (36) into 


2ms 


m So 


[anf e7 iP (2% -n) /(p) [1 — ara n)| = ote | 
(37) 


0 


x ' dn [w9(0_, n) + w,(0_, ny? ae [wolSo, n) “ We(Eo, n) 


0 


- The integral over 7 in the second term can be evaluated with the help of equa- 
_ tion (16); the integrand of the third term is 1 because of the boundary condi- 
_ tions; and in the last term w,(&, 7) is exponentially small and ‘can be neglected. 


_ We now obtain’ 


‘too-te : 


6m s 


(3) = D,(s) — 2ms mS [dpe 1) [6 6) “4 ntpyT¥ a [ne + 5] 08) 
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The term proportional to 7) in effect gives us the coefficient b, in equation 
(32), a term which we already knew was correctly given by the third term of 
equation (5a). To extend the results a little further we will now include also the 
fourth term from equation (5a), a term which we know is correct but which we 
can not find from the present calculations using equation (34). The last term 
in equation (38) combines with the last term of D,(s) given in equation (31) to 
produce the coefficient a, in equation (32). 

We are not able to evaluate the remaining integral in equation (38) expli- 
citly but we can reduce it to a more tractable form by changing the integration 
variable from to 6, equation (19), and then substituting for /[0] and /[0] using 
equations (28) and (21). The 9 integration, which is along a vertical line to the 
left of the origin is then deformed to a vertical line through 6 = —z/2 and by 
substituting 6 = 7 6’ — 2/2, we can reduce this to an integral along the real line. 

We thus obtain as our final expression for D(s) 


1 17 
=a 


nye SHEN ee 
0 


D(s) =mso*{1 Sy ts 


13a \3y3 


4 No 3 3 es 
tera 22s eid *a)) G9) 


+ (é, ys)" (ys ns) for 1%, & > 1 


in which 


On) = =p | a0 Oe ee 
0 

In addition to the errors indicated, there is an error due to the discontinuity 
in the derivative of w in equation (34). This discontinuity is of order exp(—2 noV's) 
or for s ~7 of order exp(— V2 4). The variational procedure will give D(s) _ 
correctly at least to order &>* exp (— V2 %) and it is reasonable to assume that 
the error in m~D(s) from this source is no larger than order 1 exp(—2V2 mp). | 
The results of the next section indicate, however, that this error might even — 
be considerably smaller. Since, if need be, we could easily obtain as many terms _ 
as we wish in the series multiplying /& in equation (39), the limitation on 
the use of equation (39) is, at worst, determined by the errors of order (£V/s)-8 
and &>1 exp(—2 V2 ). | 

Present applications of the oscillation type viscometer, however, call for 
a determination of the viscosity from observation of the decrement of oscilla- 
tion in which case one is usually interested only in the imaginary part of D(s). 
Also one wants to know D(s) for s = + (i — A) w and in the present situation — 
one finds that @ ~ 1 — A. For good experimental measurements, it is necessary 
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that 4S 1/10 and preferably that A ~ 1/1007 In calculating the smallest 
terms of equation (39), it suffices therefore to approximate them by their values 
a 


| 0-1 
= 
= 
= 
am 
2 of ‘aie oe 
=0+ 
cf py 
-03 
Figure 3 
The vertical scale applies to curve A. 
Curve B is the same graph magnified 10 times. 
=f}h 


We make use of this in two ways. First we note that the terms 
~ms?20(E,3 5-3?) and ms??? (4 no/£)) O(E) 3 s~ 2”) will be real and independent of s. 
They contribute nothing to Im D(s). The errors in Im D(s) therefore depend upon 
the next terms of the series namely m O(&* s~"”) and m(4 mo/&) O(& * s~*”). 
_ Secondly we note that the last term in equation (39) is at least no larger than 
order &5' exp(— V2 No). We may take s = 7 in this term and consider only the 
imaginary part. Since Q(V/s 9) is multiplied by s, it suffices to evaluate only 
Re OWi No) as a function of the one real parameter 7. The integral in question 
was evaluated by numerical integration for several values of 7 and the results 
are plotted in Figure 3. Since this term of equation (39) is very small if one is 
working in a range where equation (39) is a valid approximation, values of 


“Re OWE No) taken from this graph should be sufficiently accurate. 
ye 
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Equation (39) was derived by assuming yj) and &) >> 1 but from our crude 


estimates of the errors we can now define more precisely this range of No and bo: 
We show in Figure 4 the regions of the (79, &9)-plane in which equation (39) is 


ar | \ 
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Figure 4 
The regions above the broken lines labeled 1% and 0:1% are regions in which the thick disk © 
formulae are estimated to be correct to within 1% and 0-1% respectively. The regions below 1 
the solid curves represent the corresponding regions where the thin disk formulae apply. The — 
straight line represents a fixed geometry with €5/No = 0. 


accurate to 1% and 0-1% assuming that the fractional errors in Im D(s) are. 


approximately &>4 and &>1 exp(—2/2 7). Note that the condition 7) > 1 can 
be interpreted very loosely. 

For any given disk, 7/&) is a constant and the various possible values of 
No and. lie on a line of slope 7/& through the origin. As a practical matter, 
equation (39) is probably best suited to disks for which y/& = 1/4. In such a 
case the error in Im D(s) is always dominated by the error of order &> ‘4 and one 
can obtain accuracies of 0-1% for 5 = 7 and 1% for & = 4. 


4. Reformulation of the Thin Disk 


boundary layer thicknesses one may have a situation in which np) <1< &, 0 : 
in which 1< No <>. In the former case the formulae of the present paper 
apply and in the latter case the formulae of Part III apply. Figure 4 shows the 


If one is dealing with a very thin disk for which 7)/£ <1, then for = 
region of the (7%, &9)-plane in which the thin disk formulae of Part III apply’ 


| 
€ 
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and one sees that a line of constant 19/& << 1 willpass sain both the ‘thin 
disk’ and ‘thick disk’ regions. 

Unfortunately there are practical femttations on how small one can make 
Ho/&o and most of the experimental data do not clearly fall into ranges where 
either the thin disk or thick disk are are valid. Much of these data fall - 
into ‘the intermediate range , = O(1) < &. We will now try to demonstrate 
the possibility of interpolating dae the limiting cases and thereby obtaining 
a procedure for dealing with the intermediate case. 

We begin by writing D(s) in the form 


ee 32 | 2 / 3 7. 3 3 
D(s) = ms ji += an + 28s  4Reat+guet 
4 1 3 5] ’ 
oe Pape ar tage Na ae 
ex (no V3), &2 (me Vs), &s (m0 V8) 
dex ae ae a ays er : 


in which ¢,(7 's) etc. are as yet unspecified functions which are in fact defined 
by this equation. 

_ The first series in (e,s”?)—! is taken from the formula for the thin disk, 
Part III, whereas the series multiplying 4 7/& is taken from equation (39). 
This choice is somewhat artificial but it does have the desirable feature that 
for only moderately large values of &, the first series is correctly given and 
dominates the part proportional to /& which is incorrect for 1 ‘<1. For 
‘moderately large 7, 1 << p< &>, the error in the first series is small while the 
‘series proportional to 7 is described correctly. The burden of finding D(s), 

Bowever, falls on finding the ¢,(7% Vs) etc. and the real motivation for writing 
Ds s) in this form is to make these terms small over the entire range No << So 

while keeping the formula reasonably simple. 
a From equation (39) and equation (58) of Part III, we see that the €;(7 /s) 
Satisfy the following equations: 


J , 
g 16/42 2VE 
se ess (; ee a + O(Vs mo) + ate. ) or nol, a 
Noy S) = y 
[AB L-toe no's) — 1-00] + O¢0a for 1jo<1 | 
f [as O (e721?) | nie) 
Anos) = oe | 1b) 
| oem [— log (o/s) — 141] + O(n?) for <1 | 
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' me toler for 79 >1, 
plata ‘ 
E3 (no /s), E4 ( lo J s) lo (ys no) for 1% a is oe: 


Since these give only relatively small contributions to D(s) and we are inter- 
ested only in Im D(s), it suffices to consider only the quantities Im 7 &(7 Vi), 
Im 74/2 @,(n Vi), Im €3(% Vi) and Im i-/2 @,(7o Vi). The third of these quanti- 


0 0-2 0-4 0°6 0°8 4:0 


Figure 5 


The curve 7%) > 1 corresponds to the thick disk formula, the curve 7) < 1 to the thin disk formula, 
equation (41a). The broken line at the right indicates the magnitude of the error estimated in 
section 3. The broken line at the left represents a possible interpolation between the two curves, 


ties has the limiting values 0 and O(n) <1 and is to be multiplied by & ‘ 
whereas the fourth is at most of order 1 and is to be multiplied by &) *. We neg- 
lect these. The limiting forms of the remaining two quantities are shown in 
Figures 5 and 6 by the solid lines. 

The behavior of the curves for 7) > 1 is rather spectacular, particularly in 
Figure 5. First of all the range of 7) shown in Figure 5 lies mostly outside the 
range where equation (39) is supposed to be valid and if we assume an error 
—exp(— 2 V2 4), which would be consistent with our previous estimates, we 
would obtain the broken line on the right side of Figure 5. The curve No > 1 
extended to ™ <1 comes surprisingly close, however, to the correct curve for 
%9 <1 shown by the other solid curve. In view of this, it seems that the errors 


] 
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‘of order’ exp(— 2 V2 ) are probably only a verysmall multiple (perhaps as 
little as 1/50) of exp(— 2V2 No) One should also note that we have already sub- 
tracted away from the coefficient of (& /s) “tin D(s) an amount 2 + 4 y and 
that the vertical scale of Figure 5 runs only to 0-4. Also shown in Figure 5 by 
a broken line is a possible interpolation between the two curves 1 <1 and 


No > 1 in the region 7) ~ 0-1. From the shape of these curves it seems reason- 
able to assume that this interpolation curve may be correct to within about 


Figure 6 


7 A graph of equation (41b) with the same notation as in Figure 5. 


0-01 (absolute error). The resulting error in D(s) is, however, smaller than the 
error in Im 7 ¢, (m9 Vi) by a factor &51 <1. 

In Figure 6 we note that the range of 7 is smaller than that of Figure 5 and 

the interpolation curve is more questionable. It seems unlikely, however, that 
‘the broken line would be in error by more than about + 0-05 and this would 
contribute to D(s) an error of only (+ 0-05)/&. 
a We speculate that equation (40) along with Figures 5 and 6 or the limiting 
‘relations, equation (41), gives Im D(s) correctly to within fractional errors of 
‘roughly (0-01)/&, (0-05)/&2 and 1/& each of which is less than 0-002 for yeaah 
‘Collectively the errors in Im D(s) should be less than 0-5% for & > 5. By the 
‘Same arguments we would conclude that equation (39) itself is probably correct 
to 0-5% for &, > 5 and yy > 0-2, a range which greatly exceeds the more conser- 
vative estimates of the previous section. 
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5. Conclusions 


The main problem treated here was the evaluation of the function D(s) 
which is, in effect, the viscous drag exerted by the fluid on the oscillating disk. 
The functions D(s) is to be used in the formulae ot Part I, particularly equation 
(1,44), to obtain a relation between the decrement of oscillation and the vis- 
cosity. The actual determination of the viscosity from the measured values of 
the decrement can be done in nearly the same manner as described previously 
for the large cup in Part II and the thin disk of Part III. This evaluation is not 
significantly more difficult with the new formula for D(s) than with less accu- 
rate approximate formulae for D(s) presently in use. Equation (39), however, is 
sufficiently accurate to open the way to the use of the oscillating disk for preci- 
sion absolute measurements of viscosity. 

It is also worth noting that the method of solution for the fluid flow near 
the right angle turn obtained in section 2, can also be applied to the solution 
of certain other problems, for example heat transfer near a right angle turn. 
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Zusammenfassung 


Der Einfluss der Zahigkeitskrafte auf eine in einer Fliissigkeit schwingende 
Scheibe wird nach den gleichen Methoden wie in Teil III behandelt. Jedoch wird 
im Unterschied dazu vorausgesetzt, dass die Grenzschichtdicke klein ist im Ver- 
haltnis zur Dicke und zum Radius der Scheibe. Die abgeleiteten Formeln lassen 
sich mit gutem Ergebnis fiir die diinne Scheibe extrapolieren (hier ist die Grenz- 
schichtdicke jedoch gross, verglichen mit der Scheibendicke, aber klein gegeniiber 
dem Radius), so dass es moéglich ist, zwischen diesen Formeln und denen oe Teils [11 
zu interpolieren. Damit wird die Berechnung des Reibungsmomentes fiir Scheiben 
unterschiedlichster Form erméglicht, vorausgesetzt, dass die Grenzschichtdicke 
klein ist gegeniiber dem Radius der Scheibe. 
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Sur la résolution numérique dans les cas limites 
des €quations différentielles linéaires 4 coefficients sinusoidaux 


Par JEAN Parry, Ziirich!) 


Introduction 


Les équations différentielles linéaires 4 coefficients périodiques font l’objet 
de publications de plus en plus nombreuses, car leur importance est considé- 
rable et la détermination numérique de leurs solutions pose déja des problémes 
difficiles 4 résoudre. Il y a quelque temps?), l’auteur a présenté trois méthodes 
générales pour résoudre les équations linéaires a coefficients sinusoidaux. La 
forme des solutions et le procédé 4 employer pour calculer ces derniéres dépen- 
dent du domaine de convergence envisagé. Ces domaines sont au nombre de 
trois, sont séparés les uns des autres par des courbes de forme trés simple 
(droites ou cercles suivant les cas) et seront désignés au cours de cette étude par 
_ l’expression «domaines fondamentaux». Dans la publication susmentionnée, la 
validité des méthodes de résolution n’a été démontrée qu’a l’intérieur des do- 
maines fondamentaux. I] s’agit maintenant d’étendre ces résultats aux limites 
elles-mémes. 
Aprés un rappel des propriétés fondamentales de ces équations (chapitre 1), 
il sera démontré que les s solutions, ainsi que les s premiéres dérivées de chacune 
d’elles, sont réguliéres aux points singuliers, si une condition trés simple est 
satisfaite. Dans le cas général, on peut montrer que les séries correspondantes 
_ convergent aussi sur les limites des domaines fondamentaux (chapitre 2). La 
- méthode numérique, valable dans le domaine intermédiaire, se révéle aussi 

utilisable, si ce domaine disparait, c’est-a-dire, si les deux limites de conver- 

gence sont superposées, a la condition que les solutions restent réguliéres sur 
ces limites (chapitre 3). Si les deux rangées de points singuliers sont elles- 

-mémes confondues, la méthode peut étre légérement modifiée, afin d’améliorer 
la précision des résultats. L’étude numérique d’une équation particuliére de 
ce type (chapitre 4) permet de comparer les valeurs obtenues a l’aide de ces dif- 
férentes méthodes avec les expressions que donne la résolution analytique de 
cette équation. La correspondance est telle, que la valeur pratique de cette 
étude peut étre considérée comme démontrée. 
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NS A 
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1. Les propriétés fondamentales 


11. Les équations de base 


Les bases analytiques pour l’étude de la convergence sur la limite des do- 
maines fondamentaux se trouvent dans la publication mentionnée dans l’intro- 
duction et qui sera désignée par I. Il est cependant avantageux de rassembler 
ici celles qui seront nécessaires en indiquant éventuellement le numero dans I 
de la formule considérée, afin d’en faciliter le contréle. 

L’objet de cette étude est l’équation différentielle suivante: 


= we ve ne ae (11.1) 
oy (et he + g,e*) Tite (I:11.1) 


avec les constantes quelconques é,, f, et g,. La substitution 


7 ae (11.2). ¥ 
(I:11.2) § 
conduit a l’équation polynomiale équivalente: 
: dku (11.3) 
F,+ E,Z+ G,Z*) Z* =U 
oy Ent Ex eae Ail (1:11.3) 
avec les nouvelles constantes 
E, = DS Diem Em (E,= es) , : 
m=k ‘ 
Fi, = ad alae Pe= , Weg) 
+ Zilonle Ua ar} 
: ‘ 
Gi = pay diem Em (G, a 8s) | 
= ‘ 
L’équation polynomiale a quatre points singuliers: | 
Z=0," Z =o; 
i et higinhy (11.5)% 
pen? © s @ yf = | 
2, 2g, (1 V1 49), 45-5, yl ee (I:11.6) — 
Le paramétre ¢ est ainsi défini: 
Is 8s sa 6). 
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Dans le cas général, les deux singularités particuliéres Z ; et Z, sont déter- 
minées, en plus de (11.5), par l’inégalité 
(11.7) 


\24| < 
(1:11.8) 


Z,| . 


Selon le théoréme (I:2), chacune des singularités est un pdle, si elle ne coin- 
cide pas avec une autre singularité. En plus, si les deux points Z, et Z, sont con- 
fondus, ils correspondent a un péle, si une condition trés simple est satisfaite: 


11.8) 
Foy + 2,72, + 6, Z2=0. ( 
as a ae (£:11.12) 
En introduisant dans ]’équation (11.1) l’expression 
-++00 
- aeons (11.9) 
“ees (1:11.13) 


pour la solution, on obtient un systéme d’équations récurrentes, le systéme 
fondamental: 


se (11.10) 
A(u+ n) D,.+ Blu +n) D, + C(ut+n) Dry = 0 (1:11.18) 
avec les définitions: 
A = #i(n—1)*¢., Bly) = 1* uF e, , 
(1) p> (u )* Sx (u) p> fb & (11.11) 
She (1:11.19) 
Cay = let 1) fh 
k=0 


L’expression (11.9) se simplifie dans le domaine de l’origine (|Z| < |Z,|) et 
dans celui du point infini (|Z| > |Z,|), car D, est nul pour des valeurs néga- 
tives de m dans le premier cas, pour des valeurs positives dans le second. La 
solution s’obtient alors par récurrence directe 4 partir de Dy. Dans le domaine 
intermédiaire, la série est infinie dans les deux directions (—co <  < +00) 
et la méthode directe n’est plus applicable. Les domaines fondamentaux, qui 


_ correspondent aux domaines de convergence des différentes formes de la solu- 
_ tion sont ainsi définis: 


|Z|<|Z,| domaineI, |2|> 74 domaine II , meee 
|Z,|< |Z| << |Z,| domaine III. 


"Les limites de convergence (frontiéres des domaines fondamentaux) sont donc 


ee 8 be LO tak LA ay 


ee as 


les suivantes: 
|2| = |e*| = [44 (11.13) 
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et [2] = lee | Z| j 


Si ces deux limites sont confondues, la validité de la méthode numérique, 
valable a l’intérieur du domaine III, doit étre démontrée, car ce domaine dispa- 
rait. Dans de nombreux problémes, il est cependant nécessaire de trouver une 
solution convergente sur cette limite. 

Dans de nombreuses applications, l’équation différentielle présente cer- 
taines propriétés de symétrie, qui conduisent a des simplifications importantes 
pour la recherche des solutions. Nous considérerons d’une part les équations 
symétriques: 


(11.14) 
pep OF fe eden: 
ca fe= Gl) (I:21.1a’) 
et les équations antisymétriques d’autre part: 
aK (11.15) 
Co, = 90, fe = &x(—1)* ay (i2isag 


Les équations différentielles réelles (définies dans I:21.1b) ne permettent pas 
de simplifications importantes et ne seront pas étudiées spécialement. 

Les solutions des équations avec symétrie doivent jouir de propriétés ana- 
logues a celles de l’équation différentielle elle-méme. On en déduit les relations 
suivantes: 


> : F j 11.16) 
ie ee D n ( 
HO BO eae (1:21.3a) 
et en particulier, 
fe ly? 
H=0, D,= 2D.) be 2 Do De ( 
is os Pin 01.94 7a) 


Les propriétés de symétrie se retrouvent encore dans le systéme fondamental 
pour les équations symétriques: 


A(u+n)=C(—p—n), Biu+n) = B(R—p—n) 
et pour les équations antisymétriques: 
A(u+n)=—C(-p—n), Bwt+n)=—B(-p—n), BOO)=0. 
(21.2a0) 


Ces propriétés de symétrie ont encore d’autres conséquences qui seront, éven- 
. , | 
tuellement, mentionnées en temps opportun. 
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12. La méthode numérique 


La résolution numérique du systéme fondamental (11.10) dans le domaine 
III présente des difficultés supplkémentaires, car les séries avec 1 positif et 
celles avec n négatif ne convergent, dans le cas général, que dans les domaines 
I et II respectivement, selon le théoréme de Porncark (théorémes I: 4 et 5). 
Pour assurer la convergence dans le domaine III, il faut déterminer les valeurs 
caractéristiques de « ou d'une autre grandeur, la grandeur caractéristique. 

La méthode numérique permet de tourner ces difficultés en réalisant d’em- 
blée la convergence des séries. L’équation caractéristique exprime l’égalité des 
coefficients de FourRIER D,, calculés 4 partir des  positifs et ceux calculés a 
partir des n négatifs. 

Le procédé de calcul est le suivant. On introduit de nouveaux paramétres: 


Dey (12.1) 


D. n— 
eS Bip Sil Fae 


Les valeurs de dy et de ¢_y sont définies par approximation: 


dy = 2B (ut N) [lL -Ji- 4p NPD, 


(12.2) 
1 
ey =z B(w—N) [1-/1—4e(u—N)]. 
La fonction g( + 7) est une généralisation de la grandeur ¢: 
A(u+n)C(u+n—1) (12.3) 
PUT?) = Bian) Baan a) (1:13.2) 


Le nombre N doit étre choisi suffisamment grand pour que la précision des. 
résultats corresponde aux nécessités, mais il n’est pas infini. Le signe devant la 
racine dans (12.2) est choisi de maniére a rendre la valeur absolue des deux 
_grandeurs considérées la plus petite possible, afin de garantir la convergence. 
Les paramétres 6, et e_, pour les autres valeurs de n (0 Sn <olViuse 
_déduisent des équations récurrentes suivantes: 


A(utm+i1CWwtn)  , _ Atu—nC(w—n—}) rs 
£0, = Sires 6 Ss Bipen- th Se.4 (1:14.7) 
Léquation caractéristique 12.5) 
i B(u) = 80 + £0 (1:14.8) 


se résoud par approximations successives (voir: I, paragraphe 14). 


YP: ENOL NS 


ERA em, 
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La validité de cette méthode se démontre 4 partir du théoréme suivant 
(I: théoréme 13): 


Théoréme 1: Les trois conditions suivantes sont nécessaires et suffisantes 
pour garantir la validité de la méthode numérique: 


: (12.6) 
[chy — ay | < My ory] (I:15.2) 
(12.7) 
| — a9] < Mz || (1:15.3) 
et 
N 00 (12.8) 
DD Diy |Z") +>, |Dad) = ies 
De Wa aN ZN Sey eae as (1:15.4) 


Dans ces relations, M,, M, et Mz sont des grandeurs bornées qui peuvent étre 
choisies aussi petites que l’on veut. Les valeurs idéales sont désignées par a, et Dy, 
les valeurs approchées, calculées par la méthode numérique, par a, et par D,,. 
Enfin, il faut tenir compte de la définition des ap: 

bn Cut) Duss (12.9) 


on B (ut WB (se ae aD ee (I:15.1) 


Les paragraphes 15 et 16 de la publication I montrent que ces conditions 
sont remplies, si l’on reste a l’intérieur du domaine III. 


13. Application du théoréme de Fuchs 


La méthode directe pour calculer les solutions de (11.1) dans les domaines 
latéraux I et II (I: paragraphe 12) se base sur le théoréme de Fucus (I: théo- 
réme 3, seconde partie). Ce théoréme peut aussi étre utilisé pour obtenir les 
solutions autour des points singuliers Z, et Z., si le point considéré n’est pas une 
singularité essentielle, c’est-a-dire dans tous les cas ou Z, et Z, sont distincts et, 


en outre, si Z, et Z, sont confondus et que l’équation (11. 8) est satisfaite. 
La Ebbetitudan, 


L=1L,42 (13.1) 
dans Il’6quation (11.3) améne l’origine sur le point singulier Z, envisagé. L’équa- 
tion différentielle prend la forme 


; | 
ge , , , , 
2 | n+ Ey 2’ + GZ") (Z, + Ze seis 3, 3.2) 
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Les nouveaux paramétres Ej, Fi et Gi se déduisent des anciens A l’aide-des 
relations 


G.= Gz; B= E+ 2 Giza ? F,= Fy+ E, Zit G, yA (13.3) 

Les relations (11.5) et (13.3) permettent de calculer les valeurs de Ey, Fe etG.: 
G.=G,=8,, F,=0, | 

~ 88 (Z,, Z,) e, 7s (1 a 1 : * 2) . 


Z,, désigne le second point singulier. Si donc les deux points singuliers sont 
confondus, £% est nul. 


Selon le théoréme (1:7), une solution au moins peut s’écrire: 


| (13.4) 


“= Bae. (13.5) 


m=0 
On en déduit le nouveau systéme fondamental 


s4+1 


2 QO, (u+ m) Py, =O. (13.6) 


Les nouvelles grandeurs P,, et Q; (uw + m) jouent dans cette étude un réle ana- 
logue a celui des D,, et des A (uw + n), B (u + n) et C (wu +n) dans la méthode 
directe habituelle. Ces relations seront trés utiles au chapitre 2, ot la valeur de 
la solution dans le domaine du point singulier sera étudiée en détail. 


2. La convergence des séries 
21. Les solutions analytiques autour de Z, 


La convergence des séries sur les limites de convergence peut étre étudiée a 
l’aide de deux méthodes bien différentes. Une étude analytique de toutes les 
solutions dans le domaine du point singulier permet de déterminer, si elles sont 
toutes réguliéres dans ce domaine. I] est alors vraisemblable que les séries, dont 


le domaine de convergence est limité par ce point, convergent aussi en ce point. 
‘Une étude de la convergence elle-méme des séries est aussi possible dans bien 


des cas. Elle doit faire appel 4 un critére de convergence trés précis, car les 
critéres classiques ne sont pas suffisants. Ces deux méthodes seront utilisées 


dans ce travail aussi loin que possible et conduiront a des résultats identiques. 


L’étude du domaine autour de Z, et celle du domaine autour de Z, sont 


identiques, si l’on désigne par Z,, le point singulier envisagé et par Z;, Vautre 
point singulier. Pour simplifier l’écriture, nous ne considérerons que le point Z, 


pour Z,, Z, remplacant Z;,,. 


ce ne \ he 
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Les solutions valables autour de Z, sont définies par (13.5) et (13.6) 
L’équation caractéristique est la suivante: 


Q(u) =9, (21.1) 


et détermine les valeurs caractéristiques de qui seules sont intéressantes 1c1. 
= Fe. is b , 
Les grandeurs Q;(m) se rapportent aux termes suivantes de léquation diffé- 
rentielle: 
au 
aZ's- = 5) 


Pos = 
Ou) + EZ’ Zia + Fa Ze 


2 75 WA HED rgs-1_ Gu 4 (21.2) 
O,(p) = GZ? Z BES 222) 2 + ba! Za eee 


s—2 
+ Fo ie abbas? etc. 


Be = Dh eZee ee Soe 


soe 1 


On en déduit le théoréme 2: 


Théoréme 2: Les grandeurs Q,,(u) peuvent s’écrire, si k est plus petit que 
(s — 1), sous la forme d’un produit de (s — k) facteurs: 


Qu(H) = Mel) (u — k)-(u—k—1)---(u—s4+2). (21.3) 


Le premier de ces facteurs est de degré (k + 1) en wu et dépend del’ équation différen- 
trelle envisagée. Les autres sont définis par I équation (21.3) elle-méme. 


Ce théoréme permet de déterminer immédiatement (s — 1) valeurs caracté- 

ristiques de wv: 
y= k-1Ss—2. (21.4) 

Dans le seul cas, ot Qy serait identiquement nul, l’équation caractéristique ne 
serait plus (21.1) et le nombre des valeurs caractéristiques entiéres diminuerait. 
Comme ce cas particulier sera étudié en détail plus loin, il n’y a pas de raison de 
l’approfondir ici. 

Le théoreme 2 permet de démontrer en plus la propriété suivante: 


Théoréme 3: Les (s — 1) solutions se rapportant aux (s — 1) valeurs carac- 
téristiques entiéeres de wu définies par le théoréme 2 sont toutes de premier rang. 


Démonstration: On déduit du théoréme 2 les relations suivantes: 


Q.(m)=0 pouwr0<Sm—k<s—2—k. (21.5) 


Vol. X, 1959 Sur la résolution numérique dans les cas limites .., 43 


Les solutions sont définies par la relation (13.5), ou l’on pose: 


'm =O Pour m <0. (21.6) 


Les premiéres relations de récurrence (13.6) sont toutes identiquement satis- 
faites, quels que soient les P,,, pour 0 < m < s — 2. Larelation 


On GO, (§ —1) Pj = 0 (21.7) 


est la premiére qui permet de déterminer un coefficient P,,(P,—,). Ainsi, les 
(s — 1) premiers coefficients de la série représentant la solution, peuvent étre 
choisis arbitrairement et donnent lieu a (s — 1) solutions linéairement indépen- 
dantes de la forme (13.5), c’est-a-dire de premier rang. 

La derniére solution nécessaire dans la région du point singulier Z,, se rap- 
porte a la valeur caractéristique déterminée par l’équation 


ols) = 0. (21.8) 


Seule cette solution ou l’une de ses s premiéres dérivées, peut présenter un pdle 
en Z,. Il nous faut donc calculer cette s-iéme valeur caractéristique. 


22. Les coefficients dy», de la transformation linéaire 


L’équation caractéristique (21.8) peut s’écrire en tenant compte de (21.2): 
gots) = Z, Ei (tp—s +1) + Fy =0. (22.1) 


Comme nous l’avons remarqué, cette équation peut aussi étre identiquement - 
-satisfaite. Il faut, dans ce cas, étudier g,(). Pour cela, il est nécessaire de cal- 
culer les coefficients dj, de la transformation linéaire (11.4) avec k = s — 1 et 
k=s — 2. 
De la substitution (11.2) on peut déduire: 


du : du 
4 eel aw baal | 
du du du 
si Bis = jo 
9 ie ’ 
nd dx? me d ‘3 i t (22.2) 
au U - my U 
dx? f8 ze = 38 28 ge 9 oar 
d*u d‘u Bu > au du 
qe gee 2 age | eget + mag SE 


EON ee ee eT RAN cate ak Sie Bhs 
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Il en résulte les relations suivantes: 
i, x = 1s , 


gk +1—s 


Fn = 4 (deat & din) = (BR + 1) iT (22.3) 


; k4 + 3 RP 58 + RY op oe 
dy nig = (deans + & de, +1) =( or ID ea 3 


La derniére relation fut calculée de la facon suivante. On peut s’assurer que 
dy, x49 doit étre de la forme: 


Dnpag a (aR +b Bch? od kh ake. (22.4) 
Le terme constant est nul, car on a: | 
do.=0. (22.5) | 


Les valeurs de dy, pour k plus petit que 5 permettent de déterminer les — 
coefficients a, b, c et d. . 

Les coefficients d,, pour m > k + 2 peuvent étre calculés de la méme 
maniére, mais nous n’en aurons pas besoin. 


23. La valeur de mu, 


La derniére solution de ]’équation différentielle, celle qui correspond a la 
valeur caractéristique ,, reste bornée ainsi que ses s premiéres dérivées dans | 
le domaine du point Zj, si la partie réelle de mu, est égale ou supérieure a s. On © 
peut supposer que les séries valables dans les domaines I et III convergent, 
dans ce cas, aussi sur la limite de convergence. 

La valeur de yw, est déterminée a partir de (22.1): 


Se ee dad 


= Ft i 
bs aS are 1 a in Be . (23.1) 
Les équations (13.3) et (13.4) conduisent a la relation: | 
Fig tieBt2a te Gane : 
=Ss 1 Sek pogo oS) Pia ise 1 
i & Za Z) Zy oo 


La transformation linéaire (11.2) et les valeurs des dy, calculés au paragraphe 
précedent, permettent d’écrire: 


v Pee is Dee Ly =P Gy Ls 


J = 23.3) 
= 4 (fea +24 loca + 23 5-1) + (fo Ze, + Zig ieee ( 


Vol. X, 1959 Sur la résolution numérique dans les cas limites ... 45 


La deuxiéme partie du membre de droite est nulie, en raison de l’équation, de 
définition de Z,. 


Pour terminer, la valeur de mu, peut_se calculer dans le cas général a l’aide 
de la formule suivante: 


s=s—l Ts fst 8s—1 =) a 
g=s—1+ ee art =r (23.4) 
avec les définitions 
2 1 es ao ‘ 
B= — zt (1—/1-4@), ma (1+V1—49). @3.5) 


Les résultats d’une étude autour du point Z, sont les mémes, sauf que le signe 
devant la racine est changé. 


Théoréme 4: Toutes les solutions de l’équation différentielle (11.1) ainsi 
que les s premiéres dérivées de chacune d’elles sont réguliéres dans le domaine des 
singularités Z, et Z, respectivement, si la relation suivante est satisfaite: 


te. 4 t 8s 4 bans ‘jes.| 1 R (Af eA) il 23.6 
Re( es 2 8s 2f, tV1—4¢ ae 28, 2fs a ( ; 


Le signe devant la racine est positif pour le domaine autour de Z,, négatif pour 
celui autour de Z4. 


| 24. Le cas particulier Z, = Z, 


Le point singulier Z, = Z, n’est une singularité essentielle que si la relation 
(11.8) n’est pas satisfaite, c’est-a-dire si g9() n’est pas identiquement nul. 
De (23.3) on déduit: 


fees aE LZ, Cs_1 ets ZA tres = 0 (24.1) 
Ravec : 
& Difes wre 
Se eae © ean aoe Fe (24.2) 


On en tire la condition préliminaire 


g, et pee §s—1 3} 
e, Fag fo OT eee) 


‘Les valeurs caractéristiques de « sont déterminées par la grandeur ¢,(/): 
(4) =(u—s) (ut 1—s) 4G. + (u—s) (u+1—s)s 2, E, 


; (24.4) 
+ (w+1—s)Z, Ey + (url —s)(s—1) Fiat Fee 


x 


* 


ONT uN eT 
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Cette fonction est quadratique en mw et donne les deux valeurs caractéristiques 
qui manquent, car le théoréme 2 n’en détermine que (s — 2). 
La transformation linéaire (11.4) et les relations (22.3) conduisent aux ex- 


pressions 


Gee os 
Ey = @(2;— 23) = 0, , 
By = — tsa — 24 24 Bsa > (24.5) 4 
yas 0) 5 | 
ie —(fece + Zale ae 20 eee 
On obtient enfin: ] 
2 i Be iene GTC ; 
qx(4) = (mu — s)? ree (tS) ae +1 ve i: Te ) | pi | 
+i ss tans Siem Sa aa 


Pour déterminer si les solutions correspondant a yu, et 4,-, sont réguliéres ainsi ; 
que leurs s premiéres dérivées au point singulier, il faut tenir compte du fait 
que le plus petit exposant de Z’ est égal a (u — 1). 


Théoréme 5: Si s = 2, les deux solutions et les deux premuéres dérivées de 
chacune d’elles sont régulieres dans le domaine du point singulier Z, = Z, siVona 
My entier et positif et ; 

Re ff, > 3 (24.7) 


Oe ee ee 


ou st la. premiere relation est remplacée par 
Re a, 930 (24.8) 


Théoréme 6: Si s est plus grand que 2, la condition nécéssaire et suffisante 
pour que les s solutions et les s premiéres dérivées de chacune d’ elles soient réguliéres 
au point Z, = Z,, s écrit: 


Re iy Ss 4 Reig se (24.9) 


On peut enfin déduire de toutes ces relations le théoréme 7 pout s =i 


Théoréme 7: La solution et sa premiere dérivée d’une équation différentiell 
de premier ordre sont réguliéres dans la condition suivante: 


Re(u—1)21 (24.10) 
avec | 


&1 ey 


(u—1) = 25 (2 — +). (24.11 
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Démonstration du théoréme 7: De (24.4) on tire pour s = 1: 
u(t) = (u — 1) (ZFC! + Z, El) + Z, Bf. (24.12) 


Les relations (24.5) conduisent alors a la valeur caractéristique (24.11). 

La condition particuliére (24.7) du théoréme 5 provient du fait que la solu- 
tion correspondant a y, est certainement du premier rang, car il n’y a aucune 
valeur caractéristique en dehors de yw, et vy. Pour des valeurs supérieures de s, 
il est quelquefois nécéssaire d’introduire des solutions de second rang. La 
s-iéme dérivée d’une telle solution contient un terme de la forme: 

ds 


Boe (2 In 2’) = (uw — 1) (w— 2) ++ (w— 5) Zn Zp, (24.13) 


Ce terme est infini au point Z’ = 0 si la partie réelle de wu n’est pas plus grande 
que (s + 1), ce qui conduit au théoréme 6. 


25. Les propriétés particuliéres des équations différentielles avec symeétrie 


Les équations différentielles avec symétrie jouissent de propriétés parti- 
euliéres en rapport avec les problemes envisagés. 


Théoréme 8: Les solutions des équations différentielles avec symétrie ont 
le méme caractére (régulier ou singulier) autour des deux points singulters. 


Démonstration du théoréme 8: Ce théoréme se déduit directement des pro- 
priétés de symétrie, car les deux points singuliers sont symétriques l’un par 
‘rapport a l'autre. Il peut aussi étre démontré par l’analyse. De (11.14) ou de 
(11.15) on tire les relations: 

Say ee ee eee (25.1) 
; ‘ é, te 8s 
La condition (23.6) se simplifie: 


e a 8s-1 

y Re (22h) S 1. DN 
a E ( §s ( ) 
‘Elle est donc la méme pour les deux points singuliers. 

a 


_ Théoréme 9: Pour les équations différentielles avec symétrie, la condition 
ie do(tt) = 0 dans les cas Z, = Z, est toujours satisfaite. 
va 


 Démonstration du théoréme 9: De (24.2) on déduit: 


2,=2,=- cea s (25.3) 


5 


: 


Ss 
ta relation (24.1) est alors toujours satisfaite, si l’on tient compte de (25.1). 
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L’équation caractéristique (24.6) peut aussi se simplifier: 


qi(u) ==e (u ec) s)}? Ae (u s) (1 24 fs-3.) 24 §s-1 a8 2 eso Z £52 a 0 ; 


8s 8s §s es 5 . 
(25.4) 
Les valeurs caractéristiques se calculent a l’aide de la formule: 
1 a A Calon 4 oe au: ass ny 2 bs2 ‘ 25.5 
Lapa a a ie a 2 f, ms (25.5) 


On en tire une condition suffisante pour que les solutions et leurs s premieres 
dérivées restent réguliéres : 


Théoréme 10: Une équation différentielle d’ ordre supérieur a 2, et jourssant 
de propriétés de symétrie, a des solutions réguliéres ainst que les s premieres dérvvées” 
de chacune d’elles, si toutes les fractions de la relation (25.5) sont réelles et si la” 
condition supplémentaire suivante est satisfaite: 


3 a §s—1 1 C59 &s—2 : 
2 ie ( 8s ) = 2 % 2 é, 2 8s ; (25.08 
Ces conditions sont suffisantes, mais non nécessaires. La premiere de ces conditions 

est réalisée dans la plupart des applications pratiques. , 


26. L’étude directe des séries 
Les résultats du paragraphe 23 peuvent étre confirmés par l'étude directe 
des solutions normales valables dans les domaines I, II et III. Il est ainsi pos- 
sible d’établir une condition garantissant la convergence des séries correspon- 
dantes, non seulement a l’intérieur du domaine, mais sur la limite elle-méme. 
Si l’on demande que la s-iéme dérivée de la solution converge aussi sur cette 
limite, la condition obtenue est identique a (23.6). } 
Les critéres classiques de convergence ne s’appliquent pas a la limite-méme 


du domaine de convergence. I] faut donc se baser sur le théoréme suivant (voir I: 7 
théoréme auxiliaire 2) : 


Théoréme auxiliaire 1: La série S,(u), définie ci-dessous, converge soul 
toutes les valeurs positives de n, car sa valeur est comprise entre les deux limites 
suivantes: 

1 a 1 ene 1 (26.1 
ee yy 
nM : > h+p nate? "<9. 75g 7m 

On peut en déduire la condition pour le rapport D,_,/D,,: 


n* 


ID. p t 
np Za Ua (ie ye eee 


n—->co Dy 
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avec l’inégalité 


Repos. 1: (26.3) 
et la définition 


1 


Te aro (26.4) 


> s * ta Ne LS 5 LS , , a , 
L’opération de la dérivation a pour effet de diminuer l’exposant ¢ d’une unité. 
Pour garantir la convergence de la s-iéme dérivée, il faut remplir la condition 


Ret>s+1. (26.5) 


Nous considérons le systéme récurrent a la limite. Nous obtenons: 
7B) 
Sas =(1+#y) Z, (26.6) 
et | 
= iol A a Sos SA eee || 
+ Z, (1+ ty) [e. (ut n)* —te,4 (w+) 


(26.7) 
+ [fs (t+ m+ 1)*—tfa (ut nt ij t]=0. | 


_Nous pouvons négliger les termes quadratiques en y aprés avoir divisé toute 
V’équation par (4 + 1)s: 


Zi (b+ 2ty)[g,(L—sy)—tgeay]+ 24 (L+¢y) [ee—teay] 
+{f,l+sy)—tfa1y)=0. 


Cette équation contient un terme constant qui est nul par définition, un terme 
linéaire qui détermine la valeur de ¢ et un terme quadratique qui peut étre 
_négligé: 


t(2 g,Z3 + ¢,Z,) —s (Zig,—fe) — 4 (25 sa + 21 eoa + fea) =O. (26.9). 


| (26.8) 


_A cause de (11.5) on peut écrire: 

2¢,2%+2,2Z,= Z3g,—f,= —(4.¢, + 213) (26.10) 
Z et enfin ieee 
‘ 2g, 22+ ¢,Z,=Z,¢,1—49. (26.11) 


- Pour finir, on obtient la valeur de ¢: 


i U §s—1 a pant Ae 1 ( a §s-1 Vfoy a es-1 : 2612 
ss ea 27) it Se rae, 
Ta condition (26.5) est donc identique a la condition (23.6), de sorte que le 


AY 


+ 


théoréme suivant est démontré. 


_ ZAMP X/4 


TNT NS 


so \ a 
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Théoréme 11: Les deux propriétés suivantes sont réalisées stmultanément 
si la condition: 


Reise Ee te Re (46 : fees. fs a 


Ss 


est satisfatte: 

a) Les points singuliers Z, et Z, respectivement sont des poles pour aucune des 
solutions de lV’ équation différentielle, ainsi que pour aucune des s premieres 
dérivées de chacune @ elles. 

b) Les séries qui représentent les solutions et leurs s premieres dérivées dans les 
domaines I et III, resp. II et III convergent aussi sur la limite de ces do- 
maines. 

Cette condition n’est pas valable si les points Z, et Z, sont confondus. Elle prend 

une forme simplifiée 

Re (i f=) Sal (26.14) 


\ Ss 


st  équation différentielle est symétrique ou antisymetrique. 


On peut supposer que les propriétés a) et b) sont aussi réalisées en méme 


temps, méme si Z, est égal 4 Z,. Les conditions correspondantes furent étudiées 


au paragraphe 24. 


3. La méthode numérique dans le cas | Z,|=|Z2| 
31. Les conditions fondamentales 


La validité de la méthode numérique fut démontrée (dans I, paragraphes 
15 et 16) pour l’intérieur du domaine IIT: 


[21] < |2| < [2]. (31.1) 


Il est possible d’étendre ce domaine de validité en deux phases. Tout d’abord 
on montre que les séries valables dans le domaine III correspondent aux solu- 
tions de l’équation différentielle, aussi sur la limite de convergence, si les dites 
séries y convergent. Ensuite, on étudie la validité de la méthode numérique 
quand les valeurs absolues de Z, et de Z, sont égales. Le domaine III selon 
(11.12) a entiérement disparu. I] ne reste qu’une courbe: 


\Z,| = |Z| =|Z,| domaine IIT’ (31.2) 


sur laquelle la méthode numérique est applicable. 

Le théoréme (I: théoréme 15), qui détermine la validité de la méthode 
numérique dans le domaine III, ne peut cependant pas s’appliquer au domaine 
III’. I faut donc se baser sur les définitions (12.9) et sur le théoréme 1. 


/ 
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La premiére condition (12.6) se démontre dans le cas général (|Z,| < |Z,]) 
de la facon suivante. De nouveaux parametres », sont ainsi définis: 


cad , (31.3) 
Ly = +n)+n. 
01 (u ) F Ny (1:15.24) 
La relation de récurrence (12.4) peut s’écrire: 
hoy er) (31.4) 
uigtat Get eae 
et conduit 4 une équation récurrente pour les | 7/,|: 
, ad, 5 
| Mn—a| 2 l7n| g |O1 (u + n) = 2s. (u Cun a 1)| roe 
avec la définition: 
Qa(utn) | (31.6) 
il ir 
CA eee gern a (1:15.26) 


La démonstration envisagée n’était pas valable si 7’ est égal a 1. La fonc- 
tion 0, (u + ”) peut cependant étre développée en série de MacLAuRIN: 


Bee = ONO Te er aye (31.7) 


Si 7 est assez grand, on peut écrire: 


2k || 


wae Te (31.8) 


laut —a(ete—Di<{ 


: ou o” correspond au premier coefficient 0’, o”, etc., qui soit différent de zéro. . 
De (31.5) et de (31.8) on déduit 


Inp-al <2 [| 3 : sen aay 
a—l im, (utn+m)rti ju tn— 1|* 
> 
mal VEC 
Pee bp (31.10) 


‘si l’on tient compte du théoréme auxiliaire 1. 
Théoréme 12: Les paramétres nn diminuent donc towjours au moins 
“comme (1/1) quand n tend vers V'infin. La condition (12. 6) se trouve réalisée dans 


a les cas, méme si seul le domaine III’ substste. 


4 
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32. Les autres conditions du théoréme 1 


La seconde condition du théoréme 1 est toujours réalisée en méme temps — 
que la premiére, de sorte qu’il n’est pas nécessaire de l’étudier. La troisiéme 
condition ne fut démontrée que si le point Z est a l’intérieur du domaine III. Ul 
faut donc étudier spécialement le cas, ol. ce point se trouve sur la limite de 


convergence. 
La série dans l’équation (12.8) se subdivise en trois parties. 


n=N”+1 
ee, AUeZ 
n=N-+1 


La valeur de N” est choisie de telle sorte que la relation suivante est satisfaite: 


Si= Sf = Dal lz, (1:15.34) | 


32.2) 
6a,| Sdn N" <n=<N ( é 
[ban] S[Seuvl ( ) ia 
avec la définition : ' 
6 ane hy by (32.3) & 

i on (e522) 
Il est enfin possible de poser: t 
” Ay 
Dyn = Dyn rn (32 )e 


(1:15.36) 


La somme partielle 3’, tend toujours vers 0 quand N tend vers oo, car elle” 
se compose d’un nombre fini de termes qui ont tous la valeur limite 0. La somme 
partielle 3’, correspond au reste d’une série qui définit la solution exacte. Ce 
reste tend vers 0 aussitét que la dite série converge. 

La somme partielle 3’, doit étre examinée de facon trés attentive. De (32.2) 3 
de (32.4) et de la définition (12.9) pour «, on peut tirer: 


iene PNT ER. 9 ot My 


[D,—D,| 21D, Fier bag oe 32 : 
| 
Cette grandeur peut étre transformée a l’aide du théoréme auxiliaire 2: 


Théoréme auxiliaire 2: Dans tous les cas, l’inégalité suivante est satisfaite: 


(1+ [d ay|)""%" —1 <n |b ay| exp (m |b ay) . (32.6) 
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Démonstration: La formule bindmiale de NEWron a la forme: 


1+ Leal —1=(n—N’") | dey, +-(n — N") (n — N* —1) Pe 


(32.7) 


La comparaison avec la série représentant la fonction exponentielle conduit a 
la relation suivante sin > N": 


1+ |day|)"*" —1< (n—N") |b ay| exp {(n — N”) |6 Gilp.  hSaae) 


Les nombres » et N" sont tous les deux positifs ou nuls et le premier est plus 
grand que le second. La relation (32.8) est donc valable. Le remplacement de 
(n — N") par n augmente le membre de droite de (32.8) et conduit au membre 
de droite de (32.6). 

Les termes de »’, se composent de deux facteurs, dont le second peut étre 
transformé a l'aide de (32.6). Le premier peut aussi étre remplacé par une ex- 
pression simple, quand on tient compte de la convergence de la série: 


M 
/ We Z| = yk? > (32.9) 


k’ étant plus grand que 1. La relation (31.9) conduit a l’expression 


N 


oy |Pa— Dal [241 < 3 a (32.10) 


=N"+1 


car la fonction exponentielle de (32.6) est bornée (k = 1). M et M’ sont deux 
grandeurs qui restent bornées. On peut dans chaque cas les considérer comme 
_ des constantes. 
Le théoréme auxiliaire 1 permet de calculer la somme dans (32.10): 


MM’ L , 
y : =e [ee 2 : : N%k/-2 = NE (k a 2) , 
Z. IM’ N 
gee ae eee _ 2a 
Ee = Nk In N" - 4) : ( ) 
. M’ N2-k’ ; 
poe ONE: ae 


Spans les premiers cas, la limite est nulle quand N tend vers oo, car le facteur 

_ N* dans le dénominateur est plus important que InN. Dans i troisiéme cas, 

_Pexposant de N est (2 — k’ — k). Il est toujours négatif, car k et k’ sont plus 
grand que 1. La troisiéme condition est donc remplie, si la série converge et sila 
elation (31.9) est satisfaite avec k = 1. 


4 


$ 
? 
y 
1 
Md 
4 
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33. La premiere condition dans le cas |Z,| = |Z3| 


Le cas |Z,| =|Z,| se différencie du cas général, seulement par le fait que 
les deux grandeurs 0,(00) et 9,(00) = 1 — @,(00) sont conjuguées complexe et 
ont la méme valeur absolue. Cependant, ce n’est pas toujours le cas pour 
0,(u +n) et (4 + n), si (w+ n) reste fini. 


Théoréme 13: Les trois conditions du théoreme 71 sont aussi réalisées, st 
|Z,| est égal A|Z,| et st, par contre, Z, n'est pas égal a Zy, ala condition que les 
séries convergent en ces deux points. | 


Démonstration: La série de MACLAURIN pour @, (u + 7) a toujours la forme , 
(31.7) quand Z, ¢ Z,, car la racine ne donne lieu 4 aucun membre avec |’ex-— 
posant (1/2). La relation (31.9) fut obtenue avec l’hypothése 7’ = 1. Elle est | 
donc aussi valable si Z, et Z, ont la méme valeur absolue. Ce cas ne se diffé- 
rencie pas de facon essentielle du cas général. 

I] serait cependant intéressant de trouver un procédé spécial, afin d’amé-_ 
liorer la précision des résultats. Soit f 


(33.1) 


ey es = 
Pp (u Ar n) a (u ae n)™ (Po ) ([:13.11) 7 
et f 
= De ; 
or (ut+n) =—-[1- p+ apeeaE | @=/1-4) (33.2) 
: 
La valeur de «y sera ainsi définie: 
, af 
i [1 at BS eae (33.3) 


Le nouveau paramétre y doit étre choisi de fagon que la relation de récurrence _ 
n’en modifie pas la valeur: 


7 (1-8 + ea Po t+ %1/ (u +N) 


FFN=T)-Gmtte-roe | OA 


On en déduit: 


i 1 — 
O-P)+y(Gowor— pew) —4mt he. GBI 


Les termes constants sont égaux selon la définition de 


aa ob- 
fcnplidedateraiten ier LIES 


_ 29; 
Vee (33.6 
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car tous les termes d’ordres supérieurs sont négligés. La grandeur fb peut 
prendre n’importe quelle valeur a l'exception de 0. La relation (33.6) correspond 
cependant exactement a la formule-approchée (31.3). La méthode habituelle 
conduit donc aux résultats les plus précis, aussit6t que 6 n’est pas nul, 


Théoréme 14: Dans tous les cas, oi les points singuliers Z, et Z, sont séparés 
(6 & 0), la relation suivante est satisfaite: 
M" 
Ne 


| Ooty, —. 


(33.7) 


M" est un nombre réel, positif et borné. 


Cette inégalité est en général plus sévére que la relation (31.9). La méthode 
numérique est donc valable si les séries obtenues convergent. 


34. La premiére condition avec Z, = Z, 


Les séries valables dans le domaine III’ ne convergent, si Z; = Z,, que si la 
condition préliminaire (24.3) est satisfaite. 


Théoréme 15: La condition préliminaire (24.3) conduit a la relation: 


p=0. (34.1) 


Démonstration: Les définitions (11.11) pour A (w+), B(u+n) et 
_C (w+ n) et celle de w (u + 7) nous donne la relation suivante: 


&s-1 §s—2 | 
aa site 2 
‘< ees ee, A Se 
eee aa as REARS etn 
” at i 
ss Bi ua | (34.2) 

t Hen Y ma fs—2 

3 w+n (pe 9NP oS 
y x - 

ee fae be C52 
+n (u + n)? 


Vira os 


Le membre de droite peut se développer en série de MacLaurin: 


p(utn)=9[1 wig (fe + ges) eA. (34.3) 


wv 


i 


{aA t, 


La valeur de q, selon la définition (33.1) est donc nulle. 
Le théoréme 15 a une conséquence trés importante. La relation (31.7) est 
aussi valable quand Z, = Z, si les séries convergent, car l’expression sous la 


m 


AN 


SAVE Nee 
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racine dans la définition de 0, (u +) est du deuxiéme ordre en (“+ n). La 
grandeur & dans la relation (31.9) est au moins égale a 1 et la méthode numé- 


rique est applicable. 


Le procédé de calcul peut étre amélioré, et l’écriture simplifiée, en raison des 


relations suivantes, qui sont valables dans le cas étudié ici: 


1 il 
p=—q alo) = 5- 


Soit les paramétres suivants: 


1 il 
bn =z V1—49(u+), De i i 


La méthode normale donne: 


qv = Py 
et 
— Int 2Pn 
In-1 all ane a In 
On en déduit: 
Ulsrnt = Cis = 


(34.4) § 


(34.8) 


: mt _ ‘ s iy ae 
Cette solution n’est pas idéale, car les paramétres g, devraient étre en premiére 


approximation inversément proportionnels a (w+ 7). Il est donc avantageux 


de modifier la relation (34.6): 
In = ¥ Pw 
et de définir y de facon a satisfaire ]’équation 
Qv-1= Y Py-1- 
La relation de récurrence (34.7) est toujours valable: 


y bu + 2 bi 
tee Pe 


Y Py-1= 
c’est-a-dire: 


2? by Pyrat Y (Py—1 Di 2 py =0, 


‘oat ; ( zz - a -. i ( - fue J+ ae ' 


Cette relation peut étre simplifiée si N est trés grand: 


1 1 N N-1 1 Pn 


Pre oe Paes V 02 V% V9 Sy Dra 


(34.9) 


(34.10) 


(34.11) 
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On en déduit pour terminer: 


/ o* V 2 “ 1 , / ae 1 
ty = Yt ae t+ 46 — I. (34.14) 


Les valeurs correspondantes des paramétres q, seront désignés par g), pour les 
différencier des valeurs obtenues par la méthode numérique normale. Les er- 
reurs commises sur g,, sont du second ordre en (4 + N)-1, de sorte que le 
théoréme 13 peut étre encore étendu au cas Z, = Z,: 


Théoréme 16: La méthode numérique est applicable a Vintérieur du do- 
maine III sans restriction, elle conduit aussi a la solution sur les limites de con- 
vergences, st la série correspondante converge sur cette limite. La méthode numé- 
rique est encore valable dans le domaine ITI’ si la solution est représentée par une 
série convergente dans ce domaine. II est toujours possible de déterminer la valeur 
approchée de xy de fagon a satisfaire la relation 


M* 


Rs ectteen nent (34.15) 


Les critéres de convergence ont été exposés dans le chapitre précédent. 


4. Un exemple numérique 
41. La solution analytique 


étude d’un procédé de calcul ne saurait étre complete, s’il n’est pas dé- 
montré pratiquement a l’aide d’un exemple numérique. II est alors avantageux 
de choisir dans ce but une équation différentielle, dont laou les solutions peuvent 
étre déterminées avec toute la précision nécessaire. Une équation simple, ne. 
comportant qu’un ou deux paramétres arbitraires, présente aussi de grands 
avantages. 
L’équation différentielle suivante remplit ces deux conditions: 


(1 + 2@cosx) S* 4 2p9sineu=0. (41.1) 


Nous considérerons p et # comme réels, quoique cela ne soit pas nécessaire. 
C'est une équation du type antisymétrique qui a une solution symétrique dans 
le domaine IIT: : 


a Se A 


Mur = 9 (41.2) 


selon (I, paragraphe 21). Cette solution est unique, puisque l’équation est du 
premier ordre. 


=) 


> ie oe" hia! 


Pee 
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La solution analytique de (41.1) est la suivante: 


au _ _2p$sinx K (1+ 20cosx)?+. (41.3) | 


u = K (1+2¢%cosx)? P 


Les points singuliers sont déterminés par la relation: 


1 
COS%9 = —Fa- (41.4) 
Ces valeurs de x» sont par paires conjuguées complexe si | #| est inférieur a 0,5, — 
elles sont symétriques par rapport a l’origine, dés que || dépasse 0,5, elles © 
sont enfin superposées (Z, = Z,) quand || = 0,5. f 

Le cas # = 0,5 est particuliérement intéressant. De (41.3), on peut déduire: — 


2p 


(1 + cosx)? = 2 (cos3) (41.5) 


et par dérivation: 
du 


_ Ke 26 # a idee $ 
M  —K2 psin4 7 (cos 5 y (41.6) 


La valeur de du/dx pour x = x) dépend de la valeur de p. Dans le cas général, i 
on obtient: 


(a ee deer (ae =oo (p<1). (41.7). 
Dans le cas # = 0,5, il vient par contre: 
ae) sty (6>3), Eb =10e (e< =). b=05. (418) 


Cette différence dans la valeur limite de # au sujet de la convergence des séries 
s’est déja manifestée au chapitre 2, ot le cas Z, = Z, a di étre étudié sé- 
parément. 


j 
: 

Enfin, il faut remarquer que la solution reste toujours réelle pour x réel, : 
sauf si # est plus grand que 0,5. La parenthése de (41.3) change alors de signe 


, 


et la solution peut devenir complexe. 


42. Les méthodes générales de calcul 
Le systéme fondamental correspondant a (41.1) s’écrit: 
(42.1) 


B(utn—1—))D,.+ (ut) D,+9(ut+n+144)D 


n+l 0 . 
On en déduit immédiatement les valeurs caractéristiques : 


f=—1—p domainel; w=1+ domaine II ; ##=0 domaine III. 
(42.2) 
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La valeur caractéristique pour le domaine III se déduit de la propriété de 
symétrie de l’équation différentielle. 


La condition de convergence se déduit dans le cas général du théoréme 4: 


Re @y >" ‘10g 233} 


Cette valeur limite correspond exactement A celle que donne la solution ana- 
lytique et que précise les relations (41.7). 
Dans le cas spécial # = 0,5, on déduit du théoréme 7: 


ee py Ly pS OSs (42.4) 


la valeur déterminée par la solution analytique (voir 41.8). 
Les équations spéciales de la méthode numérique sont les suivantes: 


1) — g2 @—1—?P) (m+) pp" 
o(n) = eo (1 a | (42.5) 
et 
De. os 
8, = — O(n 6 + 1) 3 = # An ag Les (42.6) 


En raison de la symétrie, il n’est pas nécessaire de calculer les ¢,. 
On peut enfin remarquer que la solution présente une forme polynomiale 
selon I: théoréme 29, quand # est un entier positif: 


A(n=p+1)=0, Cn=-f—1)=0, (42.7) 
_car la valeur de m définie par la relation (I:22.1) est la suivante: 
M=—2p— 21. (42.8) 


Dans ce cas, le deveoppemcnt de Four!eR de la solution analytique (41.3) ne’ 
-comporte, en effet, qu’un nombre fini de termes. 


43. Les résultats numériques pour |Z,| < |Z,| 


Les résultats généraux obtenus par la théorie des méthodes numérique et 
“directe correspondent exactement a ceux d’une étude de la solution analytique. 
g ‘Une étude numérique avec les valeurs 


G= 025, Z,——0,267949, Z,— —3,732051 (43.1) 


PEEL ONT 


i) 
+ 


§=05, Z,=Z,==1 (43.2) 


ous permettra de confirmer cette correspondance. 


eiieaiees il HN 


60 JEAN PATRY ZAMP 
Le cas (43.1) ne présente pas grand intérét, sauf pour les valeurs de x: 
Jee Zale Ou ee ae ale (43.3) 


Par suite de la symétrie, les résultats sont, d’ailleurs, les mémes pour ces deux 
valeurs. 
Considérons les nouvelles grandeurs: 


D = |Ze\eD 2. (43.4) 
La solution u(x) peut s’écrire: 
+00 , 
te (aE a Se Dine END = tn eae (43.5) 


Les tableaux 1 et 2 contiennent les valeurs de D/, calculées par la méthode 
directe, les tableaux 3 et 4, celles que donne la méthode numérique. Le nombre 
N de la méthode numérique est égal a 32, mais les valeurs de D3, et de Dj2 sont © 
considérées comme insuffisamment exactes. Les valeurs choisies pour p se 
situent de part et d’autre de la limite p = 1, afin que l’influence de cette limite 
soit visible. Ces mémes tableaux numériques contiennent aussi les valeurs de 
(uw +n) Dj, et de (u +n)? Dj, qui permettent d’évaluer la convergence des 5 
séries. 

Le tableau 5 contient les valeurs numériques de D,,, calculées par la méthode ~ 
numérique. Elles conduisent directement a la solution avec x réel. Dans le 
tableau 6 sont, enfin, rassemblées les valeurs de wu et wu’ pour différentes valeurs 
de x. Les valeurs analytiques, tirées de (41.3) avec K = 1 sont désignées par 
Uy,, les valeurs calculées et normalisées par u,. Pour calculer ces derniéres, les” 
séries ont été prolongées au dela de n = 30 a l'aide d’une formule approchée. : 

La correspondance entre uy, et u, est pratiquement idéale, sauf pour les” 
valeurs de u'(Z,). La limite de convergence n’est pas respectée, car elle se situe 
a p = 1,055 dans les résultats de la méthode directe, 4 p = 0,945 dans ceux de la 
méthode numérique. Ces seules différences proviennent d’une inexactitude dans" 
le prolongement de la série. I] faut remarquer que l’exposant k, utilisé pour le 
prolongement, n’est connu qu’a 5% environ. Comme cet exposant est trés_ 
proche de 1, la valeur de ( — 1) peut présenter une erreur supérieure a 100% 


7 
; 
44. Les résultats numériques pour Z, = Z, . 


Le cas Z, égal a Z, est encore plus intéressant que le précédent, car la 
méthode numérique n’est pas applicable selon les théorémes généraux (cf. I 

paragraphes 15 et 16). Il a été démontré dans le chapitre 3 qu'il est avantageux 
d'introduire une valeur de «,, différente de celle que l’on utilise ordinairement. 
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Les tableaux numériques contiennent des résultats obtenus a partir des deux 
formules, afin de permettre une comparaison. 


La limite de convergence se situe ap = 0,5. Choisissons, pour cette raison, 
les valeurs suivantes de p: 


p= 05 et 2-06. (44.1) 


La relation (34.14) et la définition (42.5) donnent la valeur: 


PFI Los p 
/ ia We a ic 


P 1 
eat N- 
Théoréme 17: La relation (44.2) pour qy est aussi valable pour qj, quel que 


sott n. Elle est donc rigoureusement exacte. 


Démonstration du théorveme 17: Dans le cas considéré, en obtient a l’aide de 
la relation de récurrence (34.7): 


P 4 P(p+1)_ 
; 2n ' 2n(n—1) p 
qa = 1+ p/n ye ee, 1), nae) 


ce qu'il fallait démontrer. 
L’étude analytique peut étre encore prolongée dans le cas p = 0,5: 


Théoréme 18: Les coefficients Dj, de la solution exacte (valeur approchée 
Qy = valeur exacte) de (41.1) avec p = 0,5 correspondent 4 I’ équation: 


! (— AN ! 
eye eae 0" (44.4) | 
Démonstration du théoreéme 18: De (44.3) ainsi que des définitions pour 6,, 


a, et gn, on peut déduire: 


6, =1- a, = (1 —)/2.. (44.5) 
ia relation de récurrence (42.6) permet de calculer la valeur de 6,: 
ok pbs Fs 
Bans le cas p = 0,5. On obtient, enfin, le premier rapport DADs: 
: LOE OS cee ee Se Le (44.7) 
Do 1,5 3 es a | 
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Le cas général conduit a la relation: 


aes nr 2 P a Oe Fer a Ne ee =e ae 
ls alae Pia ls Lea) 2 n+i+p mo eee 


n 


En multipliant le numérateur et le dénominateur de (44.8) par 4 (x + 1/2), on 
obtient : 


Pats cones: (44.9) 


Diy 4 (n+ 1)?—1° 


La relation (44.4) se déduit de (44.7) et de (44.9) par induction totale. 
On peut encore remarquer au sujet de théoréme 18, que les coefficients in- 
diqués par ce théoréme, sont les coefficients de FourTER de la fonction 


| (1 + cosx)"? = y2 | cos (x/2) | (44.10) 


a un facteur prés. La méthode numérique modifiée conduit donc, dans ce cas — 


extréme, a des résultats d’une précision analytique. 


Les tableaux 7 et 8 contiennent les valeurs calculées pour f = 0,5, les ta- ; 


bleaux 9 et 10, celles qui correspondent 4 # = 0,6. Pour permettre une comparai- 
son entre la méthode numérique normale et la méthode améliorée, les valeurs de 
gn et de n D,, sont aussi indiquées. Les différences (g, — gn) et (n D, — nD») di- 
minuent quand  diminue, malgré que g, et  D, augmentent dans les mémes 
conditions. Les erreurs relatives sont rapidement trés petites. 

Quelques valeurs de u(x) et de w’(x) sont réunies dans le tableau 11. Les 
valeurs de u,, et u, correspondent de nouveau aux valeurs analytiques avec 
K = 1, et aux valeurs calculées et normalisées. La correspondance peut étre 
considérée comme idéale. 


45. Quelques remarques sur le cas |Z,| = |Z,| 


Aucun calcul numérique n’a été effectué avec une valeur de # plus grande 
que 0,5. Il est cependant possible d’examiner quelques difficultés qui pourraient 
se présenter. Tout d’abord, les valeurs de «, ne restent pas réelles. Il y a deux 


ee he ee 


| 


solutions possibles qui sont conjuguées complexe. Ce résultat semble en contra- 


diction avec la définition (41.3) de la solution analytique. Dans le cas # > 0,5, 


l'expression entre parenthéses change de signe. Comme # n’est pas un nombre 


entier, la fonction w n’est pas enti¢rement réelle et ses coefficients de FOURIER 


ne le sont pas non plus. Les deux valeurs possibles de a%  conduisent a deux 
possibilités pour la représentation de (— 1)?. 


Une difficulté se présente aussi si 6 est trés petit. La valeur de N doit étre_ 


trés grande, sinon la précision des résultats laisse A désirer. Cette difficulté se 
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présente toujours, si les deux points singuliers sont trés rapprochés sans étre 
confondus. Ce domaine fait la transition entre les conditions normales (23.6) 
pour la convergence des séries et les conditions exceptionnelles du paragraphe 
24. Cette transition s’effectue de fagon progressive dans le sens que les condi- 
tions exceptionnelles semblent réalisées pour les faibles valeurs de n, alors que 
les conditions normales apparaissent pour des valeurs de n d’autant plus gran- 
des, que la distance Z, — Z, est plus petite. 
La derniére difficulté qui pourrait se présenter, est liée a la relation: 


€é-=0, p(w~)=oco. (45.1) 


La méthode numérique est alors inutilisable. Seule une méthode directe dans 
le domaine du point Z, peut conduire au résultat cherché. Comme I’on a: 


Oe oe (45.2) 


tout le domaine III’ se trouve a l’intérieur et sur la limite du domaine de con- 
vergence de la série ainsi obtenue. 

Il est ainsi démontré que les méthodes directes et numériques, exposées en 
détail dans la publication I, sont applicables dans tous les cas, ot les solutions 
et leurs premiéres dérivées (jusqu’a l’ordre s y compris) sont convergentes. Ce 
résultat ne peut pas s’obtenir par l'étude directe des fractions continues. 


< 


Zusammenfassung 


_ Die linearen Differentialgleichungen mit sinusférmigen Koeffizienten weisen im 
allgemeinen drei Konvergenzgebiete auf, in welchen die Lésungen als quasi-perio- 
dische Fourier-Reihen dargestellt werden kénnen. Die Theorie dieser Gleichungen 
wird auf die Falle erweitert, in welchen die Lésungen auf die Konvergenzgrenzen 
gesucht werden. Es zeigt sich, dass die Losungen auf den Grenzlinien dann konver- 
sieren, wenn besondere Bedingungen erfiillt sind, und dass die numerische Methode 
. das mittlere Gebiet immer anwendbar ist, wenn die gesuchte Lésung konver- 
giert, auch in dem Fall, wo die beiden Grenzlinien zusammenfallen und wo deshalb 
e.: mittlere Gebiet keine Flache mehr aufweist. Diese analytischen Ergebnisse 
erlautern die numerische Auswertung eines Beispieles. 


(Regu: le 3 Février 1958.) 
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Tableau 1 


Les coefficients de Fourier Dj, par la méthode divecte 


ZAMP 


17 


aAPWwWNe 


p = 0,90 p = 0,95 
Dy | (u+n) Dy, | (+n)? D,, Dy (w+) Dn | (w+ nm)? Dp 
1,00000 — 0,90000 0,81000 1,00000 — 0,95000 0,90250. | 
0,96462 0,096462 0,009646 | 1,01821 0,05091 0,00255 | 
0,012923 0,014215 0,015637 | 0,040924 0,042970 0,04512 | 
0,013389 0,028117 0,059046 | 0,0065793 0,013488 0,02765 | 
—0,0075077 | —0,023274 | —0,072149 |—0,0036874 | —0,011247 | —0,03430 | 
0,0047605 0,019518 0,080024 | 0,0023069 0,009343 0,03784 | 
— 32889 16773 — 85542 — 15763 — 7960 —4020 | 
24120 14712 89749 11455 6930 4193 | 
— 18472 13115 — 93117 — 8705 — 6137 — 4327 Ii 
14618 11841 95912 6842 5508 4434 | 
—0,0011868 | —0,010800 | —0,098280 |—0,0005322 | —0,004997 | —0,04522 
9836 9934 10033 4552 4575 4588 
— 8291 — 9203 = 10215 sts — 4219 — 4662 
7088 8576 10377 3250 3916 4719 
6132 — 8033 — 10523 — 2800 — 3654 — 4768 
0,0005363 0,007558 0,10657 0,0002438 0,003425 0,04812 
~4727 = 7438, — 10778 S243 3225 = 4854 
4201 6764 10890 1898 3046 4889 
— 3760 — 6430 — 10995 — 1693 — 2887 — 4922 
3386 6129 11093 1520 2744 4953 
—0,0003066 | —0,005856 | —0,11185 |—0,0001372 | —0,002614 | —0,04980 | 
2790 5608 11272 1245 2496 5004 
2550 = 5381 — 11354 = 135 — 2389 — 5029 | 
2340 5171 11428 1039 2291 5052 | 
2155 — 4978 — 11499 = 955 = 2201 — 5073 
0,0001992 0,004801 0,11570 0,0000881 0,002119 0,05096 | 
— 1846 — 4633 ~ 11629 =s15 — 2042 —5115 
1716 4479 11690 756 1970 5132 # 
— 1599 — 4333 11742 — 704 — 1904 — 5150 # 
1495 4201 11804 657 1843 
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Tableau 2 
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Les coefficients de Fourier Dj par la méthode directe 


p= 1,05 p = 1,10 

n 
198 | (uu + ) Dy, | (u + n)? Dy aah (u +n) Dy, (u + n)? Dj, 
1 1,00000 —1,05000 1,10250 1,00000 —1,1000 1,21000 
2 1,12539 — 0,05627 0,002814 1,17898 —0,11790 0,01179 
3 0,105542 0,10026 0,09525 0,142158 0,12794 OFT S15: 
4 |—0,0061946 | —0,012079 | —0,023554 |—0,0117736 | —0,022370 — 0,04250 
By 0,0034261 0,010107 0,029816 0,0065043 0,018862 0,05470 
6 |—0,0020865 | —0,008242 —0,032556 |—0,0039081 | —0,015242 — 0,05944 
| 13943 6902 34165 25826 12655 6201 
8 — 9947 — 5918 — 35212 — 18254 —10770 — 6354 
9 7441 ‘ay Bt 35938 13548 9348 6450 
10 — 5770 —4587 — 36467 — 10434 — 8243 — 6512 
11 0,0004601 0,004118 0,03686 0,0008271 0,007361 0,06551 
12 — 3753 — 3734 = ThD — 6709 — 6642 — 6576 
‘al | 3118 3414 3738 5546 6045 6588 
14 — 2630 — 3143 — 3756 — 4658 — 5543 — 6596 
15 2248 2911 3770 3965 Silas} 6598 
16 |—0,0001943 | —0,002710 | —0,03780 |—0,0003414 | —0,004745 — 0,06596 
17 1696 2536 3791 2969 4424 6592 
18 — 1493 — 2381 — 3798 — 2604 — 4140 — 6583 
} 19 1324 2244 3804 2302 3890 6574 
} 20 — 1182 — 2122 — 3809 — 2049 — 3668 — 6566 
t 21 0,0001061 0,002011 0,03811 0,0001835 0,003468 0,06555 
22 —958 — 1917 = sel2 — 1652 — 3287 — 6541 
23 869 1821 3815 1495 3125 6531 
24 — 792 — 1738 — 3815 — 1359 — 2976 = (Sai117 
25 725 1664 3819 1241 2842 6508 
| 26 |—0,0000666 | —0,001595 | —0,03820 |—0,0001137 | —0,002717 — 0,064.94 
27 614 1532 3822 1046 2605 6486 
28 — 568 — 1474 = 3825 — 965 — 2499 — 6472 
29 527 1420 3827 893 2402 6461 
30 — 490 — 1370 — 3829 — 829 — 2313 — 6453 
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Tableau 3 
Les coefficients de Fourier Dj, par la méthode numérique 
p = 0,90 pb = 0,95 
Ma / 
IDSs (u+ n) Dy (u + n)? Dy, Dy (uu + 2) Dy (u + n)? Dj, 
—3 0,000006 —0,000018 0,000054 0,000003 — 0,000009 0,000027 
—2 |—0,000214 0,000428 | —0,000856 |—0,000112 0,000224 | —0,000448 
—1 0,060867 — 0,060867 0,060867 0,063947 — 0,063947 0,063947 
0 1,000000 0,000000 0,000000 1,000000 0,000000 0,000000 
1 0,84777 0,84777 0,84777 0,89065 0,89065 0,89065 
2 |— 0,041506 —0,083012 | —0,16602 |—0,021767 — 0,043534 | —0,087068 
2 0,015058 0,045174 0,13552 0,0075323 0,022597 0,067791 
4 |—0,0078558 | —0,031423 | —0,12569 |—0,0038338 | —0,015335 — 0,061341 
5 0,0048499 0,024250 0,12125 0,0023277 0,011639 0,058193 
6 |—0,0033038 | —0,019823 | —0,11894 |—0,0015659 | —0,009395 — 0,056372 
7 24015 16811 11767 11268 7888 55213 
8 — 18278 — 14622 — 11698 — 8505 — 6804 — 54432 
9 14398 12958 11662 6652 5987 53881 
10 — 11648 — 11648 — 11648 — 5348 — 5348 — 53479 
11 0,0009627 0,010590 0,11649 0,0004395 0,004835 0,053180 
12 — 8096 — 9715 — 11658 — 3677 — 4412 — 52949 | 
is 6908 8980 11675 3123 4060 52779 
14 — 5967 — 8354 — 11695 — 2686 — 3760 — 52046 | 
15 5208 7812 11718 2335 3503 52538 | 
16 |—0,0004587 —0,007339 | —0,11743 — 0,0002049 | —0,003278 —0,052454 q 
17 4073 6924 daly TAlL 1813 3082 52396 | 
18 — 3642 — 6556 — 11800 — 1616 — 2909 — 52358 | 
19 SPATE 6226 11830 1450 2155 52345 | 
20 — 2965 — 5930 — 11860 — 1308 — 2616 — 52320 } 
Pq 0,0002696 0,005662 0,11889 0,0001187 0,002493 0,052347 
22 — 2462 — 5416 —11916 — 1081 — 2378 — 52320 
23 2258 5193 11945 989 2275 52318 | 
24 — 2079 — 4990 — 11975 — 908 — 2179 — 52301 | 
25 1921 4803 12006 837 2093 52313 @ 
26 |—0,0001780 | —0,004628 | —0,12033 — 0,0000774 | —0,002012 | —0,052322 | 
27 1654 4466 12058 718 1939 52342 | 
28 — 1542 — 4318 — 12089 — 668 — 1870 — 52371 | 
29 1441 4179 12119 623 1807 52394 | 
30 — 1350 — 4050 12150 — 582 —1746 


— 52380 | 
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Tableau 4 
Les coefficients de Fourier Dj pay la méthode numérique 
p = 1,05 > = 1,10 
Dy, (M+ n) Dy | (w+ nD, Di, (u+n) Dy, | (w+ n)? D/ 
— 0,000003 0,000009 | —0,000027 |—0,000005 0,000015 | —0,000045 
0,000122 — 0,000244 0,000488 0,000255 — 0,000510 0,001020 
0,069988 — 0,069988 0,069988 0,072949 — 0,072949 0,072949 
1,000000 0,000000 0,000000 1,000000 0,000000 0,000000 
0,97481 0,97481 0,97481 1,01601 1,01601 1,01601 
0,023744 0,047488 0,09498 0,049409 0,098818 0,19764 
— 0,0074198 | —0,022259 | —0,06678 |—0,014613 — 0,043839 | —0,13152 
0,0035882 0,014353 0,05741 0,0068814 0,027526 0,11010 
—0,0021057 | —0,010528 | —0,05264 |—0,0039681 | —0,019841 | —0,09920 
0,0013810 0,008286 0,04972 0,0025686 0,015412 0,09247 
— 9752 — 6826 —4778 — 17924 — 12547 — 8783 
7237 5790 4632 13187 10550 8440 
— 5578 — 5020 — 4518 — 10090 —9081 — 8173 
4428 4428 4428 7958 7958 7958 
—0,0003598 | —0,003958 | —0,04354 |—0,0006430 | —0,007073 — 07780 
2980 3576 4291 5299 6358 7629 
— 2508 — 3260 — 4238 — 4438 — 5769 — 7500 
2139 2995 4193 3769 5277 7387 
— 1845 — 2768 — 4152 — 3239 —4859 — 7288 
0,0001608 0,002573 0,04117 0,0002812 0,004499 0,07199 
— 1413 — 2402 — 4083 — 2463 — 4187 —7118 
a2 Sb 2252 4054 et HD: 3915 7047 
—1116 — 2120 — 4028 — 1934 — 3675 — 6982 
1001 2002 4004 1730 3460 6920 
—0,0000903 | —0,001896 | —0,03982 |—0,0001556 | —0,003268 | —0,06862 
819 1802 3964 1407 3095 6810 
— 746 —1716 — 3947 — 1278 — 2939 — 6761 
682 1637 3929 1166 2798 6716 
— 626 — 1565 — 3912 — 1068 — 2670 — 6675 
0,0000576 0,001498 0,03895 0,0000981 0,002551 0,06632 
— 532 — 1436 — 3877 — 904 — 2441 — 6590 
493 1380 3864 836 2341 6554 
— 458 — 1328 — 3851 —775 — 2248 — 6518 
427 1281 3843 721 2162 6485 
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Tableau 5 
Les coefficients de Fourier Dy, 


p = 0,90 0,95 1,00 1,05 1,10 
0 1,000000 1,000000 1,000000 1,000000 1,000000 
is 0,227160 0,238650 0,250000 0,261190 0,272240 
2 — 0,002980 —0,001563 0,000000 0,001695 0,003547 
3 0,000290 0,000145 0,000000 —0,000142 —0,000281 
4 —0,000041 — 0,000020 0,000000 0,000018 0,000035 
E 0,000007 0,000003 0,000000 —0,000003 — 0,000005 
Tableau 6 ‘ 
Les valeurs de u(z) et u’(x) pour & = 0,25 Z 
Pp 0,90 0,95 1,05 1,10 
u47,(0) 1,440397 1,469896 1,530720 1,562070 
u,(0) 1,440405 1,469898 1,530723 1,562066 
1 4y,(7) 0,535887 0,517632 0,482968 0,466517 
U,(I) 0,535894 0,517637 0,482972 0,466511 
2 4p, (70/2) 1,000000 1,000000 1,000000 1,000000 
i ut, (70/2) 1,000000 1,000000 1,000000 1,000000 
5 u', (70/2) 0,450000 0,475000 0,525000 0,550000 
E u!(s/2) 0,450005 0,474994 0,524999 0,550001 
q 
S Un (Z1) 0,000000 0,000000 0,000000 0,000000 
g U(Z4) 0,000037 0,000059 0,000012 0,000032 
= un(—Z1) 1,86607 1,93187 2,07053 2,14354 
= u,(—Z,) 1,86608 1,93187 2,07053 214950 
uty (23) fo) oo 0,00000 0,00000 
wo(Z4) oo) Soph 0,4283 0,4294 
u',(—Z)) 0,72723 0,79470 0,94139 1,02100 
uo(—Z)) 0,72727 0,79472 0,94139 1,02094 
® » te(Z4) 0,000204 —0,000078 0,000103 —0,000129 
fies tt (— Za) 1,86607 1,93187 2,07053 2,14354 
$e (2) 0° 00 00 0,6676 
=u u}(—Z;) 0,72713 0,79470 0,94138 1,02088 
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Tableau 7 
Les résultats pour & = 0,5 et p = 0,5 
In Qn Dy n Dy n Dj 
0,250033 0,250000 0;333333 0,333319 05333533 
0,125082 0,125000 — 0,0666667 — 0,1333096 — 0,133333 
0,0834608 0,0833333 0,0285714 0,0856803 0,0857142 
0,0626722 0,0625000 — 0,0158730 — 0,0634476 — 0,0634920 
0,0502166 0,0500000 0,0101010 0,0504500 0,0505051 
0,0419275 0,0416667 — 0,00699300 ) —0,0418921 — 0,0419580 
360193 357143 512820 358206 358974 
315992 312500 — 392157 — 312849 — 313726 
281712 277778 309598 277653 278638 
254378 250000 — 250627 — 249532 — 250627 
0,0232095 0,0227273 0,00207040 0,0226540 0,0227744 
213601 208333 — 173914 = 207383 — 208697 
198022 192309 148149 191170 192594 
184734 178573 — 127715 — 177268 — 178801 
173279 166667 111236 165212 166854 
0,0163314 0,0156250 — 0,00097753 | —0,0154653 — 0,0156405 
154577 147060 86581 145328 147188 
146862 138889 — 77221 — 137029 = 138998 
140009 131579 69301 129593 131672 
133890 125000 — 62540 — 122892 — 125080 
0,0128400 0,0119048 0,00056722 0,0116819 0,0119116 
123453 113636 — 51680 — 111289 — 113696 
118980 108696 47282 106232 108749 
114921 104167 — 43422 — 101587 — 104213 
111228 100000 40016 97305 100040 
0,0107858 0,00961538 | —0,00036996 | — 0,0093345 — 0,0096190 
104776 925926 34305 89672 92624 
101952 892857 — 31898 — 86251 — 89314 
993595 862068 29739 83062 86232 
969763 833333 — 27785 — 80076 — 83335 
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Tableau 8 
Les vésultats pour 3 = 0,5 et p = 0,5 

n | Of, In D;, n Dy Tas 
31 0,00947824 0,00806451 0,00026021 0,007728 0,008067 
3? 927606 781250 — 24420 — 7465 _ 7814 
33 908958 757576 22962 WPAN I 7577 
34 891747 735294 — 21631 — 6983 — 7355 
35 875855 714286 20412 6762 7144 
36 0,00861178 0,00694445 — 0,00019294 | —0,006552 — 0,006946 
37 847624 675676 18265 6354 6758 
38 835110 657895 — 17316 — 6165 — 6580 
39 823562 641026 16439 5985 6411 
40 812913 0,00625000 — 15627 — 5814 — 6251 
41 0,00803104 0,00609756 0,00014874 0,005650 0,006098 
42 794082 595238 — 14174 — 5494 — 5953 
43 785798 581395 IB5SZZ 5344 5814 
44 778208 568181 — 12914 —5201 — 5682 
45 TAM2HS 555556 12346 5064 5556 
46 0,00764957 0,00543478 —0,00011815 | —0,004932 — 0,005435 
47 759227 531915 11318 4806 5319 
48 754055 520833 — 10851 — 4684 — 5208 
49 749398 510204 10413 4567 5102 
50 745262 500000 — 10001 — 4454 — 5000 
5L 0,00741610 0,00490196 0,00009613 0,004345 0,004903 
52 738422 480769 — 9247 — 4240 — 4808 
53 735680 471698 8901 4138 4718 
54 733368 462963 — 8574 — 4040 — 4630 
55 731470 454546 8265 3945 4546 
56 0,00729974 0,00446429 — 0,00007972 | —0,003852 — 0,004464 
‘Wf 728867 438597 7695 3763 4386 © 
58 728138 431035 — 7432 — 3677 — 4311 
59 727778 423729 7182 3593 4237 
60 727778 416667 — 6945 — 3511 — 4167 


q 
: 
: 
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Tableau 9 
Les résultats pour 9 = 0,5 et p = 0,6 
qn In ve n D, nv Dy 
0,300013 0,300000 0,375000 0,374994 0,375000 
0,150043 0,150000 — 0,0576922 =O, LSS 75 — 0,115384 
0,100073 0,100000 0,0224359 0,0672939 0,0673077 
0,0751043 0,0750000 —0,0117057 — 0,0468048 — 0,0468228 
0,0601374 0,0600000 0,00710703 0,0355128 0,0355352 
0,0501718 0,0500000 — 0,00473802 ; —0,0284014 — 0,0284281 
430645 428571 336649 235342 235654 
377439 375000 — 250529 — 200067 — 200423 
336148 333333 193116 173404 173804 
303199 300000 = L53035 — 152591 — 153036 
0,0276319 0,0272727 0,00124011 0,0135924 0,0136412 
253992 250000 — 102358 = 122297 — 122830 
235169 230769 85800 110963 111540 
219101 214286 — 72871 = 101398 — 102019 
205237 200000 62594 93225 93891 
0,0193167 0,0187500 — 0,00054298 | —0,0086166 — 0,0086877 
182573 176471 47511 80014 80769 
BASAL 166667 — 41891 — 74605 — 75404 
164888 157895 37189 69816 70659 
157449 150000 — S327 — 65546 — 66434 
0,0150768 0,0142857 0,00029834 0,0061719 0,0062651 
144743 136364 — 26930 — 58269 = 59246 
139289 130435 24420 55145 56166 
134336 125000 — 22236 — 52301 — 53366 
129824 120000 20325 49702 50812 
0,0125703 0,0115385 — 0,00018644 | — 0,0047320 — 0,0048474 
121931 Lee at 17158 45128 46327 
118471 107143 — 15838 — 43103 — 44346 
115291 103448 14661 41229) 42517 
112365 100000 — 13607 » — 39489 — 40821 


', & = Oh bd yee p Date 


A 


Toy Lk tae 


I 
tS 


n 


In 


0,0109668 
107180 
104883 
102761 
100799 


0,00989854 
973082 
957579 
943255 
930031 


0,00917835 
906603 
896276 
886802 
878134 


0,00870229 
863048 
856556 
850722 
845516 


0,00840913 
836888 
833421 
830493 
828086 


0,00826185 
824776 
823847 
823387 
823387 
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Tableau 10 
Les résultats pour } = 0,5 et p = 0,6 


ZAMP 


Ga Dy n DF 
0,00967741 0,00012660 0,0037870 
937499 — 11806 — 36358 
909090 11033 34944 
882352 — 10331 — 33616 
857142 9693 32312 
0,00833333 —0,00009110 ; —0,0031198 
810811 8577 30088 
789474 — 8088 — 29043 
769231 7639 28057 
750000 — 7225 — 27120 
0,00731707 0,00006843 0,0026232 
714285 — 6490 — 25389 
697674 6163 24587 
681818 — 5859 — 23822 
666666 5576 23090 
0,00652173 —0,00005313 | —0,0022393 
638297 5067 PAW IPS} 
624999 — 4838 ~- 21086 
612244 4623 20472 
599999 — 4422 — 19885 
0,00588234 0,00004233 0,0019324 
576922 — 4056 — 18782 
566037 3890 18264 
555555 — 3733 —17761 
545454 3585 17276 
0,00535714 — 0,00003446 | —0,0016811 
526316 3314 16359 
517242 — 3190 — 15927 
508475 3072 15505 
500000 — 2960 — 15096 
Tableau 11 


Les valeurs de u(x) pour & = 0,5 


p = 0,50 


1,41421166 1,41421356 
— 0,0000010 0,00000000 
1,00000000 1,00000000 
0,5000635 0,50000000 


1,51571609 
0,00000076 
1,00000000 
0,600157 


/ 
n Dj, 


0,0039246 
— 37779 
36409 

— 35125 
33926 


— 0,0032796 
31735 

— 30734 
225,92 

— 28900 


0,0028056 
— 27258 
26501 

— 25780 
25092 


— 0,0024440 
23815 

= 23222 
22653 

— 22110 


0,0021588 
— 21091 
20617 

— 20158 
19718 


— 0,0019298 
18890 

— 18502 
18125 

— 17760 


1,51571657 
0,00000000 
1,00000000 
0,60000000 
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On the Viscous Core of a Line Vortex II 
By Nicuoras Rott, Ithaca, New York, USA?) 


In a previous note [1]*), present-author discussed some aspects of BURGERS’ 
solution [2] for a vortex in a viscous fluid, which is aligned with an axisym- 
metric stagnation-point inflow. In the present note, the temperature distri- 
bution in this flow field will be calculated under certain limiting conditions. 

It is assumed that a steady axisymmetric state is reached, where all func- 
tions of state depend only on the radius ry. The energy equation is, under these 
conditions (see, e. g. [3]), 

quT F=nO+ — = (kr Z), (1) 
where o¢ is the density, 7 the temperature, S the entropy, wu the radial velocity, 
k the thermal conductivity, « the viscosity and ® the dissipation function, 

which has to be calculated for the flow pattern under consideration. Let the 
radial, circumferential and axial velocity components be 


b= by, Vuh), w=Zas, (2) 


or, the flow consists of an axisymmetric stagnation-point flow with the ‘inflow- 
gradient’ a and a circulatory motion with a velocity v(v) not yet specified. For 
the velocity field (2) the dissipation is found to be 


pO=n (a — >) + na. (3) 


The second term is due to the ‘basic’ stagnation point flow; it will be neglected 
in the limit of ‘high Reynolds number’ circulatory motion which will be con- 
sidered presently. Let v be given by BurcErs’ solution 


Ree te eee”) (4) 


Z PREIS 
; (v = w/o). The total dissipation per unit height can be calculated neglecting the 


1) Cornell University. 
2) Numbers in brackets refer to References, page 81. 
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second term in (3): 


dv v 5 
D=u/(Z-->) 2504 GF (5) 
0 
The result is af 
pan Osco 
D=+—* (6) 


independent of ju, as was noticed by Burcers. Let the value (6) be compared 
to the neglected term, which is in principle infinitely large; but let the radius 4 
be calculated for which the neglected term equals the value D given by equa- 
tion (6): 


This leads to the radius 


where the ‘viscous core’ radius is defined as 


oa 2. 6 


For high values of I’,,/y, i.e., for high Reynolds numbers, it takes a very big 
radius (compared to 7*) to give a ‘basic flow’ dissipation comparable to D.. 
On the other hand, the bulk of D is contributed by the dissipation in a region 
with a radius which is not too big compared to 7*, e.g., 0-98 D is found within 
10 y*, Thus, the proposed simplification implies that ‘infinity’ means a not too 
big multiple of 7*. 

Replacing S in equation (1) for a perfect gas by 


T dS =c,aT ~~ dp (9) 


the pressure will be computed, consistently with the previously made restric- 
tion to the vicinity of the core for high Reynolds numbers, from 


dp v2 ; 
Peet aa (10) 


(neglecting pressure changes due to the basic stagnation point flow) so that 
equation (1) becomes 


OTS: UtN" GO ie ee Aid te 
ou (cy dy — 7) =0(G-%) Sacer (aoa (11) 
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This relation will be used for the calculation of the temperature distribution. 
An alternate form can be obtained by combining equation (11) with the Navier- 
Stokes equation for the circumferential velocity component v(r). The latter 
equation is, without any simplification for high Reynolds numbers or low com- 
pressibility, 


v du\ _ d*v 1 dv ( du (dv uv 
ci | rt 4 i Mays Ha oe ae dy a 5, i (12) 
Multiplying (12) by v and adding it to (11) leads to the following energy 
equation: 
dT dv\ 1 
dy ) ¥ 


aT du 
ou (¢, sorte (A tt De i v), (13) 


ay 
The problem of the vortex core at high Reynolds numbers in a viscous com- 
pressible gas involves the simultaneous solution of equations (12), (13) and the 
continuity equation. At this point, a further simplifying assumption will be 
made, namely, compressibility effects will be considered small. The proposed 
solution uses the velocity distribution (2) and BuRGERs’ solution (4), which 
fulfill continuity and equation (12) for incompressible flow, i.e., w= —ar, 
og = const, and 4 = const. The pressure distribution under this assumption 
was already given in the author’s previous note; here, the temperature distri- 
bution will be calculated from equation (11) or (13), using again w = —ar and 
constant values for 0, “, and k, together with v given by (4). From the values 
of # and T found in these steps, the state of the gas is known and an iterative 
recalculation of the velocity field is possible. This plan leads to an expansion 
of the solution in powers of a ‘Mach number’ which will be defined below. 
Actually, in this note no iterative steps will be carried out; the only quantity 
to be calculated is the temperature, following the plan outlined above. The, 
‘result represents the limiting temperature distribution for high Reynolds 
number Re and low Mach number M. 
It is intuitively clear (and can be proved a fosteriort) that the ratio of a 
typical circumferential velocity, €. g., Vax, and the sound velocity c is an appro- 
priate definition of M. The value of the maximum velocity is 


Umax = 0°64 (14) 


ie : 


The condition of high Re and low M leads to the following inequalities (numeri- 
“cal factors are omitted) : 


‘2 y\1/2 
pete erte= (22) 


AN aN TE 
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The possible range for J’,, between these ‘Jimits’ is still very large, as the two 


limits lead to the inequality 
C2 


for the velocity gradient a. As v is of the order cl, where / is the mean free path, 
this condition is fulfilled by a wide margin for inflow-gradients of practical 
interest. 

Ina first example, let the temperature be calculated in the limit of vanishing 
thermal conductivity k, or infinite Prandtl number. In this case, the calculation 
of the temperature from equation (11) leads to a simple quadrature; v has to 
be introduced from equation (4), and for «= —ay the result is found in the 
following form: put 


a y y2 
x= Seed = ye (15) 
and 
IPSS IPS 
eee es eee oo. 
H* = 16 a2» 8 m2 v*2° (16) 


The latter quantity is of the order of magnitude of the kinetic ‘head’ associated 
with the maximum circulatory velocity; actually from equation (14), 


2 


VU ae Es 
H* = 2-44 St | (16a) 


With these notations, the temperature distribution for k = 0 is found to be 


cy T = c, T,, — H* |— (1—e-*) tiie emit (17) 


1 
x 
The result is plotted in Figure 1 together with the ‘kinetic’ or ‘stagnation en- 
thalpy’ which equals c,7 augmented by 


v2 H* 
ae =, ares (1 a Cae ; (18) 


It is seen that in the center, the (total) enthalpy is reduced as compared to 
the state at infinity; from (17), it is found that for x = 0, 


Colo = tpl, —H* , (19) 


a reduction more than twice as big as the kinetic head associated with v,... 
This remarkable result is the consequence of two opposing effects which ae 
termine 7° for k = 0: the dissipation acts as a steady heat source for the flow 
but the gas is also cooled as it drops into a ‘pressure pot’ while moving towaed 
the core; the net effect, for k = 0, is an appreciable ccoling in the center. 
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The situation is further clarified by considering the case of an incom- 
pressible fluid, i.e., a liquid, for which the results differ essentially from those 
obtained above for the ‘low Mach nuinber’ gas flow. The reason for this dif- 
ference is that equation (9) is not-appropriate for a liquid with an equation of 


0-5; 


7 _ - 
| 
i a ¢; (T-Teo) 


-05 


= 


Figure 1 


Temperature distributions near the viscous core of a vortex for zero heat conductivity (infinite 

Prandtl number). The solid curves represent the temperature and stagnation temperature (or 

_ enthalpy) for a gas. The broken curve (top) shows the temperature field in a liquid. The dash-dotted 
curve gives the temperature for an inviscid-isentropic (potential) vortex flow of a gas. 


state quite different from the gas law. The proper equation to replace (9) is 
Las 6,08)", (20) 


where c, is the specific heat of the liquid at constant volume, and volume 
changes are neglected. Thus, for a liquid, the energy equation for k = 0 is 
simply 
: dv Pe 
— arc, ass =+(5— ) ‘ (21) 


The result of the integration is now 


2 
¥ P= Gt | SR ay fare ie) bee Cl e-*)*| (22) 
: where the function 43 i 
- 2 x 
F(x) = [ (e-* — e-®*) (23) 


Peay Se eeca sy ee a 
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already appeared in the author’s previous note for the pressure calculation ; it 


is easily found that 
F(0)== log 2. (24) 


The function (22) is also plotted in Figure 1, and now, as is to be expected, the 
temperature is always higher than 7,,. Compared to the case of the gas, the 
heat sources remained the same, but the effect of cooling by expansion is 
absent. The temperature for x = 0 is now 


ty Tp = €; Ty, + (2:log 2 — 1) H* = c, T,,, 40-386 H*, (25) 


Also shown in Figure 1, for comparison, is the temperature T,,, for an in- 
viscid-isentropic gas flow; within the limits of the approximations for low Mach 
number flow, 


3, H* 
sagt =o Te (26) 


Cod pe = Cp ee 


Next, the thermal conductivity will be re-introduced. Clearly its effect will 
be a reduction of the temperature-differences found for k = 0, 1.e., less heating 
in liquids and less cooling in gases. In what follows, only the case of the gas 
will be treated. 

Introducing the Prandtl number 


aia os 
Gres Gk 
and using temporarily the variable 
dT ; 
Dae one ie v2 (27) 
equation (11) becomes, for constant »y and u = —ar, 
oar aaa dv v\2 dv 
oT Y +5 =- (at o)+2@o-noS (28) 


Putting 
We Z(r) ea 
the determination of Z is reduced to a quadrature, which can be carried ou 


conveniently with the arrangement of the inhomogeneous term indicated i 
(28), where v from equation (4) has to be introduced. A second quadrature | 
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yields 7, and the final result, after some calculation, is 
Rey ova, r= A* oe zs la —=Za) F(x) — re —e~*) | 


ax : ; 1x 


+ H*'o (20-1) | e-°* = f (1 —e-*) ere 


J * cae 

x 0 
The result fails for o = 2, which, however, does not represent an interesting case 
and can be calculated separately if needed. An interesting simplification is 


Figure 2 


Temperature and stagnation temperature for a gas with Prandtl number 1/2. 


found for the value o = 1/2, when the term containing the double integral in 
(29).vanishes, and the explicit result is, using the function F defined by (23), 


bd] 


A Poe ls SS: E ic ae (3) (30) 


a C 


This relation is plotted in Figure 2. The cooling effect in the middle is reduced, 
compared to the case o = oo; using the limit (24), the result is now 
Ss eM ee eee & log 2) H* =c, T,, —0-462 H*. (31) 


For Prandtl number 1/2, the ‘non-recovery’ at the core still approximately 
equals the maximum kinetic head! Writing 


Kg dgasClecie He) Ad ® ys) 
“a ‘non-recovery factor’ is defined, which is a function of the Prandtl number. 
‘Besides for o = co and o = 1/2, can be calculated easily for o = 1, The inner 
integral in the last term of (29) can be evaluated as a series: 


: : 

s ‘ dx x? at 

4 IG) =) (et 22) = 2 (sort aap + ~} 
a i 

b= 

a 


ea 
x 


ica bain’ 
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and the double integral becomes 


eee ae —— 


so that 


This value together with (1/2) = 0-462 brackets the result for all gases within — 
narrow limits. 

The final result is that the steady line vortex aligned with a stagnation 
point represents for all gases (with 2/3 So < l)a ‘refrigerator’. It is charac- 
teristic for this phenomenon — as a glance at the stagnation-enthalpy curves ; 
shows — that a small amount of cold gas is produced at the expense of a large” f 
amount of ‘lukewarm’ gas. Clearly the stagnation-enthalpy has to be considered $ 
for this final appraisal, and its maximum excess over c, J, is far smaller than ~ 
its greatest defect below its value at infinity. To say that the defect is formed : 
at the expense of the excess is literally true: the flux of the quantity 


; ; 

n(1) = 2log 2— ~ = 0-564. "t 
; 

: 

; 

2 


re 


Glaeser tint de 


; 


in the axial direction, across a plane z = const, can be shown to be zero for 
all Prandtl numbers. This flux, as w = 2 az = const on the control surface, 
equals 


CO 


2nw fc T+ . Weve, ) y dy 


0 


1 
= wW ( T + oe Bs ie} 7 


dT d 
—nw{ (cy +1 5) ew 


The first term vanishes at both limits, and the integral can be shown to vanish 
by using the energy equation (13) with y = const and uw = —ar. | 

At this point it is unavoidable to think of the case of the Ranque-Hilsch 
tube. A regular solution of the Navier-Stokes equations has been found yield- 
ing refrigeration effects of the order of the maximum kinetic head; nevertheless, | 
the applicability of the results to the Ranque-Hilsch phenomenon is very doubt- 
ful. First, the results so far do not explain the ‘hot’ end of the vortex tube; 
only modest values of heating were found. Second, the results are restricted to 
the limiting cases of high Reynolds number and low Mach number, and the 
latter restriction is very serious for the discussion of most of the vortex tube 
effects. The case of the ‘high M’ vortex is soluble in principle, but not analyti 
cally. A solution of these cases with machine calculation has been initiated in 
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dependently by Donapson [4] (who also gave several analytical limits) and 
is already under way. But even if these restrictions are removed, it is not clear 
whether calculations based on a strongly simplified flow geometry can explain 
the vortex tube effect. In particular, the strong experimental evidence (with 
theoretical support) must be accounted for, which has been gathered at the 
‘Institut fiir Aerodynamik’ (ETH, Zurich (5-7), in favor of an explanation 
of the Ranque-Hilsch effect based on unsteady flow phenomena. It is even 
possible that both kinds of effects must be considered for a satisfactory 
explanation. 

A further point can be clarified by considering the amount of refrigerated 
gas found in the solution above: it is of the order 


waurt?—2azv. 


This is a small amount depending on the kinematic viscosity but independent 
of the inflow gradient a. There seems to be no way to relate this result to ex- 
periences with the vortex tube. 

Thenew result provided by this theory is that recovery rules from boundary- 
layer theory cannot be applied to the viscous core of a vortex. The main dif- 
ference lies in the réle played by the pressure in the two cases. 
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4 Zusammenfassung 

- Die Temperaturverteilung wird berechnet im zahen Kern einer Wirbelstré- 

mung, welcher eine axisymmetrische Staupunktstromung iiberlagert ist. Im Grenz- 

fall kleiner Zahigkeit und geringer Zusammendriickbarkeit werden geschlossene 

Lésungen fiir verschiedene Prandtlsche Zahlen angegeben. Tm Wirbelkern ergibt 

sich ein bemerkenswerter Kiihlungseffekt; somit wird gezeigt, dass die Erfah- 

rungen der Grenzschicht-Theorie nicht auf den Fall des zihen Wirbelkernes an- 


wendbar sind, 
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Estimation of the Critical Viscous Sub-Layer in Shock Wave 
Boundary Layer Interaction . 


7 


By Ayir Kumar Roy, Bangalore, South India?) 


1. Introduction 
; 
One of the main complications caused by the boundary layer in a super- 
sonic flow arises from the fact that a portion of the layer is subsonic. Any 
disturbance to the flow at any point on a body will effect the flow upstream 
through the agency of the subsonic part of the boundary layer which forms | 
a channel for these disturbances to propagate along the body. Although the 
effects of the disturbance die down asymptotically because of the viscous | 
action in the boundary layer, in practice, distances as a measure of upstream 
and downstream effects of the disturbance are defined well enough. 
Experiments in the shock wave-boundary layer interactions were first 
made in Switzerland by ACKERET, FELDMANN and Rorr [2]}?) and in the USA 
by DonaLpson [1] and LrrpmMann [3]. Later on, more detailed study of the 
interaction problems with laminar and turbulent boundary layers on flat plate 
were made by LIEPMANN, RosHko, and DHAWAN [4] and others. In the present 
paper, experimental datas have been taken from the work of LIEPMANN, 
RosHkKo, and DHAWAN [4}. ; 
On the theoretical side, an exact mathematical theory is given by LicHT- 
HILL [5], regarding the upstream influence through shock or expansion wave- 
boundary layer interaction, for small pressure rises. In his analysis, the bound- 
ary layer has been divided into viscous sublayer and outer inviscid layer. In 
the inner sublayer, disturbances are effected by viscosity and the Mach 
number remains small enough for compressibility to be neglected. aa 
mathematical analysis for small perturbations as applied in hydrodynamica 
stability problems has been made in this case of flow without separation and 
the logarithmic decrement of the upstream influence under zero-heat transfer 
condition is estimated which confirms experimental results. 
The present paper is concerned with the estimation of the viscous sub-layer 
a fraction of the total boundary-layer thickness, which as in many other 
boundary layer problems [6] plays a significant role in such phenomenons and 


1) Department of Aeronautics, Indian Institute of Science, 
*) Numbers in brackets refer to References, page 89. 
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for that purpose, similarity considerations, which serve a very useful tool, 
have been applied. 


Vi 


cs\ioa\ Orie iS clad 


2. Notations 


velocity of sound; 


distance of the upstream influence measured from the point of incidence of 


the disturbance on the flat plate; 
constant in the temperature—viscosity relation 


y | EP ie + 4 (assumed 1:5 for the numerical evaluation in the pre- 
Tx 


RS 
SUTHERLAND’s constant; 
unknown function given by equation (1); 
half the length of the plate; 
distance of any point on the plate from the leading edge; 
eh — d); 
density ; 
static pressure ; 
gas constant 
= 8-31 x 10’ erg/°K-mol; 
absolute temperature; 
= 283° K (assumed in the present case for numerical evaluation) ; 
= 293° K (assumed in the present case for numerical evaluation) ; 
= Ap|bo; 
one atmosphere in the present case; 
velocity along the plate within the boundary layer; 
viscosity coefficient ; 
= plo; 
boundary layer thickness; 
thickness of the critical viscous sublayer 
= (% — Xp) VF as given in equation (9) ; 
ratio of the specific heats of the air (assumed 1-4) ; 


sent case under zero-heat transfer condition) ; 


tangent of the angle which a straight line makes with the positive direction 


of the abscissa. 


Subscripts 
Reservoir conditions; 
conditions at the wall; 
conditions at the outer edge of the boundary layer. 
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3. Theoretical Considerations 


Any disturbance caused to the boundary layer (laminar or turbulent) lead- 
ing to a positive (in the case of shock wave) or negative pressure gradient (in 
the case of expansion wave) causes the boundary layer to thicken or thin 
down [5]. In either cases, boundary layer begins to curve slightly upstream 
and the curvature of the boundary layer itself induces an additional pressure 
gradient in the same sense. This in turn, produces further thickening or thinning 
respectively of the boundary layer. The process repeats itself, pressure gradient " 
decaying gradually as we follow upstream. Theoretically, the decay is an 
asymptotic one although in practice, upstream influences are defined well 
enough. 

If Ap be the error induced in the pressure distribution along the flat plate | 
at any point x because of the disturbance incident on it, then we expect the 


error in the pressure reading 


: , 
Ap=zowF, (1) 


i.e., proportional to the dynamic pressure @ w?/2 and also a function of the 
Reynolds number 


Re = (x — x) — - 


(i.e., a characteristic Reynolds number formed with the distance * — %9) only, 
as we neglect the effect of compressibility. The final result will show whether we 
shall be justified to do so. Since the velocity in the boundary layer tends to 
zero at the wall it follows that the error tends to zero as the length (x — x9) 
sufficiently small. The distance y from the wall which influences the propaga 
tion of the disturbance must be determined by experiment and from it, th 
effective velocity w for the given boundary layer profile u(y). Let us conside 
the case when the velocity increases linearly in this region. This approximatior 
is certainly correct when the distance (x — x») is small enough. 
We thus put 


e. 


du 
a a const . (3 


If we denote the velocity at a distance y = (x — x,) from the wall by 
Ne aes Sa (4) 


we expect the similarity of the processes when we have the same value of th 
Reynolds number 

Re = Al—*%) _ (§= 4%)? de 

v v aim ( 
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The error in the pressure reading is then 
1 = 1 
Ap == 0 UZ F=—=9 (« — x)? F, (6) 


where F is the proportionality factor depending on Reynolds number i.e., 


F-F Baber ae |. (7) 


Furthermore, the function F will be dependent upon the dimensionless para- 
meter d/l where d is the distance of the wedge from the plate. Again, we can 
rewrite equation (6) in the form 


A =50[& >| , (8) 
where the characteristic length 
6, = (x — %) yr. (9) 


This characteristic length 6, obviously denotes the thickness of the critical 
region of the boundary layer (namely, the critical viscous sub-layer of L1cHT- 
HILL) which is the determining factor in the propagation of disturbances up- 
stream. 


4. Numerical Evaluation 


The tangential stress at the wall is given by 


to= [us (10) 


In the laminar case, we have 
Pe ee Mey jae 
2 2:918 V1 % 
v,% C- y—1 
8 = 5-836 |/ = = [1 + 0-417 M3 7 | (see Low [7]). (11) 
1 


du ta | 1 ieee 
eile 2-918 My Vo m% | 


Substituting this value of [dw/dy], in equation (6), we find, 


4 Ap kx 2 (Polo) (@o/@s) 1 1 i M2} (2-918)2 
_- 5 a [(#— mo) Ge] ~ “te — a)? uj C? | "6 i 
a edd ” dy 

a 


86 Ajit KumAR Roy 


under zero heat transfer condition when 


My Ty 


C 
Mw Ty [1 ea = mi] 


and putting k = Ap/p, and @ = 4. 
Further, we rewrite 


F = Jee Cy , 
where 
fatal Rx 
‘ee G — 4%)? 
and 
2 (2-918)? R To(go/e:) i [; , 1 
Gre aoe ES = Mi}. 
Again, 
(~ — x)? C32 yg ; 
eee = Re SGo- 
Vx (2-918) [1 a Be mi] p32 : 
where 
Re’ = (% a Xo) 
Ez 
and 
‘ai C82 us? 1 
Cs e-OTe.< 


(1 ST mM?) yg 
6 
In the turbulent case, 


2)", 6 = 0-38 Re-t? » 


os 


Ty = 0-037 0, u2 Re? ( 
and 


Rea . (see LIGHTHILL [5]) . 
1 


Hence, under zero-heat transfer condition 


Ree aaah 2 R T(oo/e3) C? v8 M227 2 
ue (4 — x)? (0-037)? 41876 [1 * é| tag 
where 
R xe 
FF’ => 
(¥ — %)? 
and 


— 2 RT4(0/e1) C? v3 M2)2-7 
oF (0-037) 2 1818 : [1 35 ‘| : 
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Further 


if we put 


mand 
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4% —%,)* 0-037 ugls 
Re = (x %o) ah een en ~~ Re’ G 
yils ¥ 271.8 ; 2 
v C E + ate 35 9is 
_ (% — %)? 
Redon 
pes 0-037 ugs 
? M2\1-35 aa : 
1g Se odie 


5. Result 
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(16) 


The adjoining figure shows the graphs between log F’ and log Re’ in the 
different cases. With the help of the conversion formulae (13), (14), (15) and 
(16), corresponding relations between log F and log Re have been obtained. 


log Re’ 


12-10 -08 -06 -O4 -02 0 Gay ces 06 08 


10 


12 


e@ laminar and O turbulent (reference Figure 19); x turbulent (reference Figure 21); 


. 
— 


as OS 


log F’ plotted against log Re’. 


~w laminar and VY turbulent (reference Figure 22). 


Reference Figures 19 and 22 [4] are concerned with laminar and turbulent 
boundary layers, which are subjected to an incident shock generated by a 
wedge having the semi-aperture angles equal to 45° and 3°, placed in free 
stream of Mach numbers M, = 1-44 and 1-4, the Reynolds number (based on 
the distance from the leading edge) of the flow being 0-9 x 108 in both the cases 
and the distance d of the wedge from the plate being 2:5 and 2 cm, respectively, 
1 = 6 cm in both the cases. 

Reference Figure 21 [4] is adopted from the measurements of FacE and 
SARGENT. Here a shock generated from a concave corner is incident on a turbu- 


ca) Walia th 
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lent boundary layer on a flat plate, the Mach number of the free stream being 
M, = 1-47 and Reynolds number of the flow 6 x 108, 7 = 3 cm. 

The graphs are straight lines and their inclinations do not remain constant 
throughout but change which is obvious because of the changing character of 
the boundary layer. Boundary layers, which are either laminar or turbulent 
initially, may be turbulent after disturbance; there may be separation of the 
flow from the boundary either completely or a separated bubble followed by 
the reattachment of the flow to the boundary. 

In the laminar case (reference Figure 22), ‘m’ of the straight line is at first 
1/5 and then it gradually changes to zero value. The value of 6, in the case 
m = 1/5, is found to 


(17) 


Og = (% — Xo) es | : 


y313 x 10° 


According to the reference Figure 19, in the laminar case, ‘m’ of the straight 
line is —1/12, and afterwards it changes in the same way as in the above case 
(reference Figure 22). With this value of m = —1/12 the formula for 6, in this 


case becomes 
Re-124 
. ='\5 — 4 if 
Mo) lar x ra ve 


In the turbulent case, ‘m’ of the straight lines changes in the same way, namely 
at first m= 1/5 and then it changes to m= —1/5, in all the three cases ex- 
cepting the case of reference Figure 21, where m = 1/8. 

The corresponding formulae for 6, in the respective cases are given by 


8, = (x 9) | a es | | 


y/1-145 « 108 


and (reference Figure 21) 


2 


Re-110 
(x — Xa) ——— 
244 4x “a 

Reito 


Be 
yews] | 
|] 


and (reference Figure 22) 


, 


= (x — 


Reto 


742 x mh | 


Rene | 
14-74 x oar 


15-42 x 102 Ti 


and (reference Figure 19) . 


Xo) 
(x — Xp) 


“ll 
oe 
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These different formulae for 6, have been evaluated for corresponding 
values of x and x, in the respective cases (both laminar and turbulent). It has 
been found that the values of 6, lie between the limits 0-01 mm to 0:04 mm, in 
all the cases. Further, when expressed in terms of the undisturbed boundary 
thickness, the ratio of this critical viscous sublayer to the total boundary- 
layer thickness is about 1/10 in the laminar case, and 1/100 in the turbulent 
case. 

These results support LIGHTHILL’s conjecture that this critical viscous sub- 
layer is very small and is well inside the so-called laminar sublayer in the turbu- 
lent case. 
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Zusammenfassung 


Wenn eine kompressible Grenzschicht (laminar oder turbulent) einer Stérung 
unterworfen wird, sei es durch eigene Instabilitat oder von aussen aufgepragt, so 
macht sich der Einfluss dieser St6rung sowohl stromaufwarts wie stromabwarts 
in der Strémungsgrenzschicht bemerkbar. LIGHTHILL nimmt in seiner mathema- 
tischen Behandlung des Mechanismus der Fortpflanzung dieser St6rung strom- 
aufwarts und stromabwarts an, dass die laminare Unterschicht der Grenzschicht 
diinn genug ist, um Kompressibilitatseffekte darin vernachlassigen zu konnen. 

In vorliegender Veréffentlichung wurde die Dicke dieser kritischen laminaren 
Unterschicht abgeschatzt, ausgehend von Ahnlichkeitsbetrachtungen und unter 
Verwendung vorhandener Messresultate. Es wurde festgestellt, dass das Verhaltnis 


zwischen der kritischen, zihen Unterschicht zur totalen Grenzschichtdicke (un- 


gestért) im laminaren Fall von der Gréssenordnung 1:10 und im turbulenten Fall 


etwa 1:100 ist. 


(Received: June 14, 1958.) 
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to a Pair of Vortices Moving Symmetrically 


By Miss Durca Roy, Calcutta, India?) 


Introduction 


The resistance experienced by a circular cylinder placed in a uniform flow 
with a pair of symmetrical vortices is well known. 

In actual problems, when a cylinder is introduced into a uniform flow, the 
vortices are generated symmetrically, and at great distance from the cylinder they 
move down, thus satisfying the phenomenon of the Karman vortex street. : 

In this paper we propose to calculate the varying resistance as a symmetric 
pair of vortices move down. The motion is no longer steady and the consequent 
effects have been traced. 


Resistance of a Circular Cylinder Due 


Description of the Present Work 


Complex potential when there are two vortex filaments of vorticity + and 
—x at z) and z) and a cylinder of a radius c is introduced with its centre at origin 
into a uniform flow, u along + ve # axis is found with the help of the circle 
theorem and is given by 


int es c2 oe ia ae 
Wa oe [Hos (¢ — %) — log(z — %) — log (]- 40) + log (<- z,)| 


2 
1) (< + =] " 
z 
Denoting the velocity components of the vortex at A(z)) by uw, v9, we try to 


obtain uo, Vp as the components of the velocity of the fluid at A. 
From (1) 


(1) 


1% 
Wa Iq 108 4) + W, , 


where, ignoring a constant, 
1 we 


moo ay [la 7+ ele) ota 


and the complex velocity of the vortex at A is the value —dw,/dz when z = z, and 
is given by ‘ : 


dw, | 


—= U 1 VU, = Zz, 
0 
dz = 9 9 


Pe Le | eRe Zo % 2 (2) 
EG APO EE Foe elt *(1- =) 


1) Presidency College. 
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taking complex 


: ed 1 Zo Zz co 
Zo = “ = ae Seopa 
8 te noe tae aeale ih 
The pressure is obtained from BERNOULLI’s equation as 
p 1 os op 
: + | iii const , 


where p = p,, q® is the square q of the velocity, and @ the velocity potential. 
To evaluate the constant, we assume 


p = x (const), -q=u and ies 


ot 
Hence 
1 om 
=a 2 ~ 2 
pHa zou og? + ve 
where 


= (=) (SS) ae oe See 
oz } \ 02 ot ot ot 


with the help of (4). The total resistance is given by 


x-iv=Jie/ (Qe) % ie | SP (se) > (Blasius theorem) , 
a Yr 


ot \ Oz 


at 


(4) 


(5) 


where X and Y are the components of total resistance in the direction of the 
axes and the integration is performed over the boundary + of the cylinder which 


is given by z Z = c?, where c is the radius of the cylinder. 
From (4) we get 


X—-i¥=(X,-iY) + (%,-1¥%,), 
where 
i Don [i / Ow \2 
X,-i%i=zie/ (=) Oz 
fa 
x? 9 6? = et 
=——g5, "ot *0) Gz — 68) (af — 08) ( — 08) 
y2 
ie 2 70 — 40 
12% OQOUC aa 
and 


$ ; dm /% Pees 
X,-i%=—iof SP (se) ntvoct (4-H). 
Y 


Substituting (7) and (8) in (6), we get 


ri 2 2 
fy _j;y-_ 2° 


Zia, 


(0+ %) Gz 


ee oes VET: A baa bY Nae 


(6) 
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Substituting for z) and z, from (2) and simplifying, we finally get 


xX -iY=-= 2° (+ %) 


ala 


aaj ' 2328 


y 
ie 
Ale) 
c D x 
“) 62) 


Picture of the cross section of the circular cylinder and the position of the vortex filaments on the 
(x, y) plane. U uniform stream along x axis; CDE cross section of the circular cylinder on the 
(x, y) plane; A, B position of the vortex filaments on the (x, y} plane; O origin. 


Let z, = R e’” then from (9) in polar forms 


? 4 sin? » (R? + c?) 
Be aR Aap cee 


Remarks 


Y = 0,1. e., the total force exerted on the cylinder is in the » direction. 
When the vortices are at a great distance from the origin, 


l.e., R>1thenX=0. 
When they lie on the boundary of the cylinder i.e., R = c, then 


2x? @ cos 
ae 


ee 


If the positions of the vortices are known at any time, the total resistance X 


that will be experienced by the circular cylinder can be found out with the help 
of (9) or (10). 
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Zusammenfassung 


Es wird versucht, den veranderlichen Widerstand, welchen ein unter Bildung 
zweier Wirbel in eine gleichformige Strémung gelegter, runder Zylinder erfahrt, 
- ’ 

einer Berechnung zu unterziehen. 


(Received: June 14, 1958.) 


Comments on H. Bergh and B. van den Berg’s paper: 
‘On the Visualization of Laminar Boundary Layer Oscillations 
and the Transition to Turbulent Flow’ ') 


By Wicmot H. Wess, Princeton, N. J., USA?) 


The authors have succeeded in obtaining unusually clear photographs of an 
elusive subject, a boundary layer in transition. A careful examination of their 
published photographs reveals the remarkable similarity between the transition 
vortex structure for an airfoil as observed by the authors and that for pipe flow 
as observed earlier by the undersigned®). The authors may also wish to compare 
with the work of Hama‘) in which the loop vortex structure was observed to form 
from a line vortex shed by a semi-circular rod spanning the breadth of a flat plate 
at its surface, and located near the leading edge. In the work in pipes, it was 
found that the loop structure was the inevitable link between disturbance and 
transition no matter how the disturbance was introduced, providing only that the 
disturbance did, in fact, produce transition. It is interesting to note that the 
authors have observed a similar transition process by impinging sound waves on an 
incipiently unstable boundary layer. The loop vortices appear, with the airfoil at 
zero incidence, in Figures 4e and 5a. For the case of an angle of attack, with a 
thicker boundary layer, the loop vortices become more extended and thus are 
more easily seen. They are observed to occur in all of the angle of incidence 
pictures as the authors have noted. : 

The possibility of carrying out hot-wire studies concurrently with the visuali- 
zation studies may enable the authors to correlate transition traces such as those 
obtained by ScHUBAUER and KLEBANOFF*) with their own observed flow patterns 
and thereby provide a much needed comparison between the results of these two 
means of transition exploration. 


(Received: August 11, 1958.) 


5) G. B. ScHUBAUER and P. S. KLEBANOFF, Contributions on the Mechanics of Boundary Layer 


Transition, NACA TN 3489 (Sept. 1955). 
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Herbsttagung der Schweizerischen Physikalischen Gesellschaft 
vom 13. bis 14. September 1958 in Glarus 


Berichte tiber die Sitzungen fiir angewandte Physik und Mathematik 


Statische und dynamische Eigenschaften supraleitender Zinnspulen. 
Von F, von Batimoos, Ziirich*). 

Anlass zu den folgenden Untersuchungen gab die Tatsache, dass Elemente der 
Schaltungstechnik, wie Verstarker*), Gleichrichter, Schalter’) usw., bei Helium- 
temperatur aus supraleitenden Spulen aufgebaut werden. 

Es wurde bereits an anderer Stelle*) darauf hingewiesen, dass diese Elemente 
gegeniiber der konventionellen Réhren- und Transistortechnik Vorteile bieten. Eine 
genauere Untersuchung von supraleitenden Spulen schien daher gerechtfertigt. 


sn Probe 1 


-100 -50 0 50 100 Oe 
Figur 1 


k/R;\, in Funktion von H, mit H, als Parameter. 


Figur 1 zeigt die Widerstandsiibergangskurven fiir eine Zinnspule der folgenden 
Abmessungen: Spulenlange 25 mm, Spulendurchmesser 10 mm, Drahtdurchmesser 
0,5 mm. Wenn kein ausseres Magnetfeld H, angelegt ist, beginnt der Zwischen- 
zustand, wenn das Spulenfeld H, den Wert des kritischen Feldes H, fiir Zinn bei 
dieser Temperatur erreicht hat. Damit ist erstmals ein elektrischer Widerstand 
messbar. Durch Anlegen eines dusseren Feldes lasst sich das Feld im Spuleninnern 


1) Institut fiir kalorische Apparate und Kaltetechnik der ETH, 
a J. L. Orsen, Ein supraleitender Gleichrichter und Verstirker, Z. angew. Math Phys. 9a, 298 
1958). me 
ron D, A. Buck, The Cryotron — A Superconductive Computer Component, Proc. I. R. E. 44, 482 


4) P. GRASSMANN, Was bedeutet die Tieftemperatur-Physik fiir di ik? i 
eee tig ra y fiir die Technik?, Schweiz. Bauztg. 
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verstarken oder schwachen. Im ersten Falle wird der Einsatzpunkt des Zwischen- 
zustandes um den Betrag des dusseren Feldes zu kleineren Spulenfeldern verschoben 

Ist das aussere Feld dem Spulenfeld entgegengeschaltet, dann wird das Feld im 
Spuleninnern geschwacht, auf der Aussenseite jedoch verstarkt. Das kritische Feld 
wird daher zuerst auf der Aussenseite der Spule erreicht. Dies erklart den Verlauf 
der Ubergangskurven bei Feldschwachung, 


Figur 2 


a Stromverlauf. Der Strom I, verursacht das kritische Spulenfeld H,. U ist die Spannung iiber 
der Spule, wenn der Ubergang supraleitend-normalleitend 6 sprungartig, c gemiass den statischen 
Ubergangskurven erfolgt. Im Oszillogramm d ist die Induktivitat der Spule beriicksichtigt: 

Uy = Lal |dt. 


Das Wechselstromyverhalten der Spulen wurde mit Hilfe von sagezahnformigen 
Strémen untersucht. Figur 2 zeigt die zu erwartenden Oszillogramme. Die Induk- 
tivitat ist fiir viele Anwendungen st6rend. Daher wird das Spuleninnere mit einem ° 
Supraleiter ausgefiillt, der infolge des Meissner Effektes alle magnetischen Feld- 
linien ausstésst und so die Induktivitat der Spule sehr stark verkleinert (die Per- 
meabilitat des Spuleninnern ist dann Null). 

Figur 3 stellt ein aufgenommenes Oszillogramm dar. Der Spannungsverlauf weicht 
gegeniiber dem vorausgesagten ab, indem bei den Ubergangen Normalleitung- 

‘Supraleitung und umgekehrt zusatzliche Spannungsimpulse auftreten. Diese sind 
darauf zuriickzufiihren, dass bei Zerstérung der Supraleitung die Feldlinien in die 
Spulendrahte eindringen. Die dadurch entstandene Flussénderung induziert eine 
Spannung in den Windungen (an den Stellen 1 und 2 der Figur 3). Beim umgekehr- 
ten Vorgang werden die Feldlinien wieder aus den Drahten ausgestossen, wobei 
wieder eine Spannung induziert wird (an den Stellen 3 und 4). Die Hohe dieser Stor- 
impulse ist der Frequenz proportional. Sie nimmt zu mit dem Durchmesser der 
Spulendrahte. Dadurch ist die grésste zu verarbeitende Frequenz gegeben. 

Bei einem Drahtdurchmesser von 0,1 mm liegt die Grenzfrequenz bei etwa 
1000 Hz. Durch Aufdampfen wendelférmiger, diinner Schichten lasst sich diese 
‘Grenze wesentlich nach oben schieben. 
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Herrn Prof. Dr. P. GRassMANN danke ich fiir seine wertvolle Unterstiitzung bei 
der Ausfiihrung dieser Arbeit. 
Diese Untersuchungen wurden dur 


finanziell unterstitzt. 


ch den Arbeitsbeschaffungsfonds des Bundes 


| 


Figur 3 
Beobachtete Strom- und Spannungsoszillogramme, Frequenz 200 Hz, Stromamplitude 1,5 A (das 
entspricht einem H, von 37,6 Oe); H, = 0; Amplitude der Stérimpulse 2,15 mV; 
Uy betragt 0,8 mV. : 


Aufnahme der Hystereseschleife von diinnen magnetischen Schichten. 
Von K. E. DRANGEID, Ziirich®). 


Einleituns 


Bei der Herstellung von diinnen magnetischen Schichten sind Instrumente fir 
die Messung der magnetischen Eigenschaften erforderlich, wobei in erster Linie die 
Hystereseschleife von Interesse ist. Es handelt sich um Schichten in der Dicke von 
100 A und mehr. Die magnetische Induktion ist in der Gréssenordnung von 10! G. 
Fiir kleine Schichtstiickchen haben wir daher sehr kleine Fliisse, die gemessen wer- 
den miissen. Ein grosses Schichtstiick wiirde die Erzeugung eines ausgedehnten 
homogenen magnetischen Erregerfeldes verlangen, was ebenfalls Schwierigkeiten 
bereiten wiirde. Man kénnte natiirlich eine MeBspule mit vielen Windungen ver- 
wenden, um das Signal zu vergréssern. Die Spule wiirde jedoch bald zu gross sein, 
und es ware zum Beispiel nicht méglich, wahrend des Aufdampfens zu messen. Die 
Messung wahrend des Aufdampfens bedeutet, dass die Anordnung im Hochvakuum- 
behalter untergebracht werden muss, wo nebst der Verdampferquelle noch zahl- 
reiche andere Vorrichtungen Platz finden miissen. Wiinschenswert ware eine MeB- 
spule mit einer einzigen Windung. Die gewohnlichen Messanordnungen (ballistische 
Galvanometer, Flussmeter) sind nicht empfindlich genug. Um eine Spannungs- 
verstarkung zu erreichen, gibt es zwei Méglichkeiten: Verstarkung der Signal- 
spannung mittels eines Rohren- oder Transistorenverstarkers, oder Erhéhung der 
Spannung mittels eines Transformators. Die letztere Methode ist die von uns an 
gewendete und hier beschriebene. 


5) International Business Machines Corporation, Forschungslaboratorium. 
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Das Integriernetzwerk 


Verschiedene Griinde haben die Wahl eines Transformators als Spannungs- 
verstarker bestimmt. Die Stéreinfliisse k6nnen mit einfachen Mitteln auf ein Mini- 
mum beschrankt werden, und die Verwendung eines passiven Elementes biirgt fiir 
grosse Stabilitat. Durch das Vermeiden von Verstarkerréhren bleibt man, vor 
Rauschen und vor Schwankungen in V erstarkung und Speisespannung verschont. 


l R, 
U, R, Uy Uy C uy 
a b 
U,<u, 
a ie 
fe uy _ Ur Re 
o PR c pe 


Figur 1 
Das Prinzip des Integriernetzwerkes. 


Ich méchte jetzt das Prinzip der Integration erklaren. Die Schaltungen a und b 
in Figur 1 sind einfache Integrationsglieder. Die Bedingung fiir gute Integration ist 


U2 


Man kann aber zwei Integrationsglieder kombinieren, und jetzt ist die Bedingung 
fiir die richtige Integration 
R,R,C | 1 

dees? 
Die Ausgangsspannung ist in der letzten Gleichung ausgedriickt, und man sieht, 
dass die Integration so wirkt, als ob nur ein Integrierglied vorhanden ware; ent- 
weder L und R, oder R, und C. Wird aber U, in der Berechnung beriicksichtigt, 
zeigt sich die Kombinationsschaltung giinstiger. Fiir dasselbe Ausgangssignal ist 
der Fehler bei der Integration viel kleiner, wenn die Kombinationsschaltung be- 
niitzt wird. Wenn wir einen Transformator beniitzen, um die Signalspannung zu 
erhéhen, und diesen mit einer Induktivitaét, 2 Widerstanden und einer Kapazitat 
belasten, wird das Netzwerk, unter Vernachlassigung der Hauptinduktivitat, durch 
das Schema von Figur 1 beschrieben. Streuinduktivitat und Wicklungswiderstand 
sind in L und R, inbegriffen. Durch diese Art von Integration muss man keine 
-grossen Anforderungen an den Transformator stellen. Das Integrierglied ist in der 
folgenden Art sehr leicht abzustimmen (Figur 2). 
An die Eingangsklemmen werden kurze positive und negative Impulse angelegt. 
Man sieht auf dem Bild, welche Wirkung C auf die Ausgangsspannung hat. Mit dem 
richtigen Wert von C ist das Ausgangssignal eine Rechteckspannung. 


On < Oe 
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Diese Methode haben wir schon in einem Messgerat verwendet, das die Hyste- 
reseschleife von gewéhnlichen Magnetkernen aufnimmt. 


U; 4__,_t_,_1_- t 
f—\__ 4} _f— »t Czuklein 


Up +t CzgroB 


| | : | -t  Crichtig 


Figur 2 


Die Ausgangsspannung fiir verschiedene Werte von C. 


“VOOOT™ 


on hod 


10 V= 


Figur 3 
Anordnung zur Aufnahme der Hystereseschleife von Magnetkernen. 


Als Leiter fiir den Magnetisierungsstrom dient die Nadel (Figur 3), auf die de 
Magnetkern aufgesteckt wurde. Die induzierte Spannung wird derselben Nadel ent 
nommen und dem Integriernetzwerk zugefiihrt. Die Frequenz des Wechselstrom: 
ist 16,5 kHz, und die Ablenkspannung fiir den KO liefert ein Stromwandlerim Trei 

erkreis. In diesem Geraét wird ein gewohnlicher Topfkern als Transformator ver 
wendet. Es werden zwei zusatzliche Wicklungen auf dem Stromwandler angebracht 
die eine, um Spannungen zu kompensieren, die von den ohmschen Verlusten in de 
Treiberleitung herriihren, und die andere, um Spannungen zu kompensieren, die i 
der Verdrahtung induziert werden. 
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Die Amplitude der Kompensationsspannungen kann mittels Potentiometer 
varuert werden. Die von den ohmschen Verlusten in der Stromleitung herriihrende 
Spannung ist in Phase mit dem Strom. Man muss also die Korrekturspannung vor 
der Integration abarthen, Die in der Verdrahtung induzierte Spannung ist gegen- 
uber dem Strom um 90° verschoben: Man muss also diese Spannung nach der Inte- 
gration subtrahieren. Wenn kein Kern eingesteckt ist, muss das Bild auf dem KO 
ein horizontaler Strich sein. 

Bei der Aufnahme der Hystereseschleife von diinnen magnetischen Schichten 
kommt eine 4hnliche Anordnung zur Verwendung (Figur 4). 


Figur 4 


Anordnung zur Aufnahme der Hystereseschleife von diinnen magnetischen Schichten. 


Die Schicht hegt auf einer Erregerspule. Nur das aussere Feld wird beniitzt, 
und die Abmessungen der Spule wurden der Grésse des Messobjektes angepasst. 
Das Feld ist iiber der ganzen Schicht homogen. Die MeBspule umfasst das Mess- 
objekt, und eine Kompensationsschleife ist unter der Erregerspule angebracht. Bei 
richtiger Kompensation darf ohne Schicht keine Spannung induziert werden. Die 
‘Korrektur kann auch durch die gleichen Massnahmen, wie eben gezeigt wurde 
(Figur 3), erfolgen. In diesem Gerat sind ausser der Kapazitaét und einem Wider- 
stand (Figur 4) die Streuinduktivitat und der Wicklungswiderstand massgebend fiir 
die Integration. Der Grund hiefiir ist, dass zur Erzielung eines méglichst grossen 
Ubersetzungsverhiltnisses sehr diinner Draht verwendet werden muss. Der Wick- 
lungswiderstand geniigt dann fiir die Integration. Die Beniitzung des dusseren 
Feldes der Erregerspule erméglicht das Messen wahrend des Aufdampfens. 

Die Figuren 5 und 6 zeigen die Hystereseschleife einer Schicht von 1000 A Dicke, 
1 cm Breite und 6000 G Sattigungsinduktion. Figur 5 zeigt eine Aufnahme, wenn 
‘die Schicht in der magnetischen Vorzugsrichtung (easy direction) magnetisiert 
wird. Die MeBspule hat hier 4 Windungen. Wie man sieht, ist das Bild sehr schartf 


i gee Fe ear oe ave ise gil ss hota a 
amma Tels Se ad Pai SA es 2A el 
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Figur 5 Figur 6 
“Magnetisierungskurve einer Schicht, die in Magnetisierungskurve einer Schicht, die quer 
_ der Vorzugsrichtung magnetisiert wird. zur Vorzugsrichtung magnetisiert wird. 
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und weist keine Stérspannungen oder Schwingungen auf. Hier hat das Integrier- | 
glied eine Zeitkonstante von angenahert der 5fachen Periodendauer. 

Figur 6 zeigt eine Schicht, die quer zur Vorzugsrichtung (hard direction) magne- . 
tisiert wird. Der Ausgang des Integriergliedes ist ohne Zwischenverstarker direkt | 
an den Eingangsverstarker des KO’s angeschlossen. Die Empfindlichkeit der KO’s 
ist 5 mV/cm. 1 cm entspricht also 3,5-10-1° Vs oder 3,5:10-2 Mx. Wenn man sehr 
diinne Schichten messen will, besteht natiirlich die Méglichkeit, gréssere Signal- 
spannung durch das Anschliessen eines gegengekoppelten Verstarkers zu gewinnen. 

Stéreffekte im Transformator (Barkkausen-Spriinge, Nichtlinearitéaten usw.) ~ 
legen eine obere Grenze der Empfindlichkeit des Gerates fest, die wir noch nicht 


untersucht haben. 
Kritik der Schaltung 


Zum Schluss soll auf einige Schwierigkeiten hingewiesen werden. Da man fiir 
sehr diinne Schichten das Ubersetzungsverhdltnis so hoch wahlen muss, dass fiir die 
Integration zum wesentlichen Teil die Streuinduktivitat und der Wicklungs- 
widerstand massgebend sind, muss auf die Kleinhaltung der Wicklungskapazitat 
geachtet werden, da sonst stdrende Einschwingvorgange auftreten. 

Ein wichtiger Punkt ist, dass der Erregerstrom méglichst sinusf6rmig sein sollte. 
Die harmonischen Schwingungen werden sonst mit dem Ausgangssignal tiberlagert, 
und da das Ausgangssignal sehr klein ist, machen sie sich sehr stark bemerkbar und 
sind fast nicht wegzubringen. 


Obere Grenze der Methode 


Ein Nachteil dieser Methode ist, dass man keine Leistungsverstarkung zur Ver- 
fiigung hat. Die gebrauchte Leistung muss dem magnetischen Feld entnommen 
werden. Da der Eingang des KO’s eine gewisse Eingangsimpedanz hat und da die 
Schicht sehr schnell schaltet, muss man aufpassen, dass die Anstiegszeit nicht zu 
lang wird. Durch diese Forderung ist die Ausgangsimpedanz der Schaltung dann 
gegeben. Der Strom in der Me8spule muss auch unterhalb eines gewissen Wertes 
gehalten werden, weil sonst das Magnetfeld zu stark beeinflusst wird. Diese Ein- 
schrankungen setzen das maximale Ubersetzungsverhiiltnis fest. 


Bemerkung : Es soll noch bemerkt werden, dass bei der Wahl vom Kernmaterial 
fur den Integrationstransformator die Grésse der Anfangspermeabilitat von grosser 
Bedeutung ist. Da die Signale ausserst schwach sind, rufen sie Magnetisierungsfelder 
von weniger als 50 mOe hervor, und bei vielen Materialien bleibt die relative Per- 
meabilitat in der Gréssenordnung von 1. In diesem Falle kann der Transformator 
ebensogut durch einen Lufttransformator ersetzt werden. 


Sur une équation classique de la dynamique des réacteurs atomiques. 
Par J. Patry, Ziirich®), 


Introduction 


8) Réacteur SA. 
*) P. Scumip, RAG-Bericht Nr. 5. 
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de facgon approchée de la variation de réactivité. On en déduit quelques phénoménes 
peu connus, qui ont été constatés lors d’une étude numérique pour le Diorit§) 


La formule de Nordheim généralisée 


P. ScHMID’”) a établi l’équation suivante, qui permet de déterminer la variation 
de puissance f(t) en fonction de la variation o(t) de la réactivité. 


t 
P(t) / p(t’) 
—=1+ Gir) | t’) — o) ~~ — (0(0) — ]«. 1 
i (¢— #) | (e’) ae a =) (1) 
0 
La fonction de GREEN G(t — ?#’) est définie a partir de sa transformée de LAPLACE: 
con 1 8 a; 8 
G(s) = ——____—_ — = j}'_# » Gt) = Sra, e~ 7 (2) 
\ B; s tao” es 0 
tT—oO+ So 
t=1 B s+ A; 


Les données numériques sur les neutrons retardés et sur la «durée de vie moyenne» 
des neutrons sont représentés selon l’usage par f;, B, A; et t; o est un paramétre 
arbitraire que l’on peut supposer égal 4 —2, mais qui ne joue pas de réle essentiel 
dans cette étude. Son influence fondamental est que tous les y,; sont plus grands © 
que zéro si o est négatif. 

La formule de NorDHEIM s’obtient en substituant (¢ — 7’) par ¢” et en dérivant 
par rapport au temps. On introduit ensuite la forme exponentielle pour la puissance: 


Ph Pelt. (3) 


Aprés une intégration, il s’en déduit la relation 


Ee ; - yc ee (4) 
eo ame ees cal 


a cause de la relation (2) entre les f;, B et 4; d’une part, les a; et y; d’autre part. . 
La formule de NorpHEIM généralisé s’écrit donc 


8 
= B; tila Set At aes 5 


ll faut cependant remarquer, qu’elle n’est valable que si T n’est pas trop petit 
algébriquement et si le terme 9’T est petit par rapport aux autres. Dans les cas 
contraires, elle ne donne qu’une valeur approchée de la période. 


|e 


Bréve étude de la nouvelle formule 
L’équation (5) peut s’écrire 


@= B(T)+ 


\W 
Noa 


== : 
8) J. Parry, RAG-Bericht Nr. 8. 
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La fonction B(T) représente l’influence des neutrons retardés et peut étre co ae 
comme une courbe trés générale. La figure permet de déterminer @ es en a 
tion des autres paramétres. La valeur de ¢ est la distance entre l’une des courbes 
supérieures représentant 


O1 =] BD) > (courbes continues) 


et l’une des courbes inférieures représentant 


0, = I (droites pointillées). 


Diagramme pour la détermination de la période. 


Il est facile de déduire de cette figure les conclusions suivantes: 


figure, est connue, car elle est déja indiquée par la formule classique (voir vecteurs 
Oe A). 

b) Si @” est positif, il est toujours possible d’obtenir une période positive, méme 
avec une réactivité négative (vecteur B). Il est connu que la puissance d’un 
réacteur augmente dés que les barres de contréle se déplacent vers l’extérieur, 
méme si la réactivité est négative. 

c) Si 9’ est positif, la période sera raccourcie. Cet effet est d’autant plus intense 
que la période est plus longue (voir vecteurs O et C, O, et C,). 

d) Si 9’ devient de plus en plus négatif, les valeurs de T données par la figure aug 
mentent, puis deviennent imaginaires. Les conditions de validité de l’équation 
ne sont plus remplies, car T devient en réalité négatif. La puissance diminue 
malgré que la réactivité reste positive (voir vecteurs O et D) 


Les effets c) et d) ont été observés par l’auteur lors d’une étude numérique sur 
le comportement du Diorit dans un cas spécial (la réactivité externe augmente 
linéairement avec le temps, l’influence thermique empéchant un accident) 8). 
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Note bibliographique 


A la fin de mon article Contribution a la méthode des équations aux diffévences 

[Z. angew. Math. Phys. 9a, 179 (1958)], j’ai discuté briévement le fréquence fonda- 

mentale A, d’une plaque encastrée de vété 1; il faut ajouter que S. TomoTIKA [ Phil. 

Mag. (7° série) 27, 745 (1936)] avait déja, en appliquant la théorie de A. WEINSTEIN 

[Mém. Sci. math., 88 (1937)], obtenu en collaboration avec celui-ci les évaluations 
excellentes 

1295 53.4; = 1302 
JosErH HERScH 


International Symposium on Stress-Wave Propagation in Materials 
June 30-July 3, 1959, at the Pennsylvania State University 


An International Symposium on Stress-Wave Propagation in Materials is to 
be held at the Pennsylvania State University on June 30, July 1, 2, and part of 
July 3, 1959, to be sponsored by the Office of Ordnance Research, US Army. 
Specialists abroad from England, Germany, Israel, India and USSR have been 
invited, in addition to specialists in the US, embracing the fields of high speed 
photography of stress waves, pulse-propagation experiments, scabbing in plates, 
theoretical methods, viscoelastic materials, and related subjects. Joint sessions and 
events will be planned with the eastern section of the American Seismological Society, 
which is meeting at Penn State at the same time. Further information can be ob- 
tained from Norman Davips, Professor of Engineering Mechanics, The Pennsyl- 
vania State University, University Park, Pennsylvania, who is in charge of the 
program, N. Davips 


Symposium on Millimeter Waves, March 31 and April 1, 2, 1959 
in New York City 


Millimeter waves will be the subject of the 9th international symposium of 
the Polytechnic Institute of Brooklyn, Microwave Research Institute, to be held 
in New York City on March 31, April 1, 2, 1959, under the co-sponsorship of the 
Air Force Office of Scientific Research, US Army Signal Research and Develop- . 
ment Laboratory, Office of Naval Research, and the Institute of Radio Engineers. 

The symposium is intended to highlight the present state of research in, and 
applications of, millimeter-wave technology. Accordingly, the program will be 
devoted to invited and contributed papers treating the generation, transmission, 

control, measurement, and detection of millimeter-wave energy. In addition, source 
material and significant advances in basic supporting fields will be summarized 
‘in tutorial papers chosen from appropriate fields in physics and engineering. 
~All correspondence should be addressed to: Professor HERBERT J. CARLIN, 
‘Microwave Research Institute, 55 Johnson Street, Brooklyn 1, New York. 


‘Tenth International Congress of Applied Mechanics 1960 at Stresa 


The Tenth International Congress of Applied Mechanics will be held in the 
‘Congress Building at Stresa (Italy) from Wednesday, 31 August, through Wednes- 
lay, 7 September 1960. 


‘ 
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Apart from a number of invited general lectures, the technical sessions of the 
congress will be held in two sections, viz.: 


Section 1: Fluid Dynamics (hydrodynamics and aerodynamics). ex 
Section 2: Mechanics of Solids (rigid body dynamics, vibrations, elasticity, 
plasticity and theory of structures). 


It should be noted that thermodynamics and computational methods as 
such are not included, although specific applications of computational methods 
to pertinent problems of one of the two sections mentioned above are acceptable 
subjects for papers to be read at the Tenth Congress. 

Previous congresses have demonstrated the desirability of an adequate 
period of time for the presentation and discussion of individual papers. In order 
to allow a period of 45 minutes for each paper (30 minutes for presentation 
and 15 minutes for discussion) a Program Committee will make a selection from 
papers submitted for presentation. Abstracts of papers should be submitted im 
four copies to the Secretary of the International Comittee (Prof. Mekelweg 2, Delft, 
Netherlands) before 7 January 7960. Preferably they should not exceed two 
type-written pages (doubly-spaced) and in no case should they exceed four pages. 
In order to facilitate the work of the Program Committee it is recommended that 
abstracts be in two of the official congress languages (English, French, German 
and Italian). Authors are urged to make their abstracts as clear as possible, 
since selection of papers has to be based upon them. Decisions of the Program 
Committee ave final, and it will be understood that it is impossible to enter into 
correspondence about them with authors of papers. They will be informed 
promptly of the decision on each paper. 

Day-to-day organization of the congress is effected by the Italian Organizing 
Committee (President: Prof. G. COLONNETTI; Secretary: Dr. F. Roitia, Consi- 
glio Nazionale delle Ricerche, Ufficio velazioni internazionalt, Piazza delle Scienza 7, 
Roma). All correspondence (apart from submission of papers) should be addressed 
to the Italian Organizing Committee. Information on accomodation, also regis- 
tration forms, will be obtainable from Dr. Rota on and after 1 September 1959. 

The Executive Committee of the International Committee for the Congresses 
of Applied Mechanics: C. B. BirzENo, President; RicHarD V. SOUTHWELL; 
W. T. KoiTeEr, Secretary (Prof. Mekelweg 2, Delft). 


Internationale Konferenz fiir numerische Informationsverarbeitung, 
15.-20. Juni 1959 in Paris 


__ Die UNESCO fiihrt vom 15. bis 20. Juni 1959 eine Tagung in Paris durch 
tiber das Thema der numerischen Informationsverarbeitung. Es sollen speziell 
folgende Gegenstande behandelt werden: 


1, Mathematische Methoden der numerischen Analysis. 
2. Logische Struktur der Rechenautomaten. 

. Formelsprachen fiir die Rechenautomaten. 

. Automatische Ubersetzung von Sprachen. 


; Allgemeine Probleme der Speicherung und Verarbeitung von Informationen. 
. Allgemeine Struktur von Geraten fiir die Automation. 


Auk WwW 


Wahrend der Tagung kann eine Ausstellung von Geriten besichtigt werden. 
E, STIEFE 
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Errata 


Einordnung : der Gleichungen der inkompressiblen, reibenden und 
schweren Flissigkeit. Von P. GrassMANN (ZAMP 98, Fasc.5/6, 307 [1958]}) 


Auf Seite 310, rechte Spalte, 5. Zeite von oben, muss es heissen: «Bernoullische 
Gleichung: Eu + 1/Fr = (¥) (wi/w® — 1). (Dabei ist wy diejenige Geschwindigkeit, 
die sich fiir die Héhe h = 0 und Ap = 0 ergibt).» 


a i 
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Kerbspannungslehre. Von H. Neuser, 2. Auflage (Springer-Verlag, Berlin 
#958). 226 S., 162 Abb.; DM 36.-. 

Das Buch ist die umgearbeitete und erweiterte zweite Auflage der Neuberschen 
Kerbspannungslehve von 1937. Es enthalt eine systematische theoretische Be- 
handlung der elastischen Spannungskonzentrationen in gekerbten Staiben unter 
den vier Beanspruchungen. Das Problem ist nicht nur theoretisch interessant, 
sondern im Hinblick auf die Haufigkeit derartiger Konstruktionselemente fiir den 
Praktiker von grosser Bedeutung. Die Behandlung ist klar, die Gliederung iiber- 
sichtlich und umfassend; die Mehrzahl der gegebenen Lésungen, die bis zu leicht 
verwertbaren Diagrammen ausgefiihrt sind, ist das Werk des Verfassers und seiner 
Mitarbeiter. 

Es seien, ohne das Buch als Ganzes in Frage zu stellen, einige kritische Bemer- 
kungen erlaubt: Die Erklarung (S. 164ff.) der bei Spitzkerben beobachteten 
theoretisch-experimentellen Diskrepanz aus der Struktur des Materials (und die 
Korrektur der theoretischen Ergebnisse durch Mittelbildung tiber Elemente mit 
unbekannten Dimensionen) scheint dem Referenten unbefriedigend, da ja 
schliesslich auch in Kristallen Abweichungen vom konstanten Spannungszustand 
und sogar Spannungspole méglich sind. Hier spielt doch wohl eher die Tatsache 
mit, dass auch Spitzkerben praktisch auf ihrem Grund etwas abgerundet sind, 
sowie die Plastifizierung des Materials. Gerade das plastische Verhalten wird aber 
nur gestreift (S. 208), und die Ergebnisse, die in dieser Hinsicht vor allem in den 
angelsachsischen Landern erzielt worden sind, bleiben unerwahnt. Schliesslich 
diirfte die Weiterverfolgung des im Anhang (S. 223 ff.) diskutierten nichtlinearen 
Elastizitatsgesetzes zu Schwierigkeiten fiihren, da dieser Ansatz das unterschied- 
liche Ansprechen der in Frage stehenden zahen Werkstoffe auf hydrostatische 
Spannungszustande einerseits und Spannungsdeviatoren andererseits nicht in 
Rechnung stellt. H. ZIEGLER 


Thermodynamics of One-Component Systems. Von W. N. Lacey und 
B. H. Sace (Academic Press, New York 1957). 376 S., 91 Fig.; $8.00. 
Dieses aus langer Lehrerfahrung erwachsene Buch stellt sich die Aufgabe, dem 
mit der Thermodynamik schon einigermassen vertrauten Ingenieur einen Weg 
zu ihrer praktischen Anwendung zu weisen. Dem stellen sich erfahrungsgemass 
folgende Schwierigkeiten entgegen: 
1. Der Ubergang vom geschlossenen zum offenen Prozess. 
2. Die Verwendung eines technischen Mafsystems. 
3. Die Einbeziehung der Reibung. 
4, Die rechnerische Erfassung des in Wirklichkeit immer sehr komplexen Vorgangs. 


ERR ASS ayy ret tO 
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Dem ersten Ziel entsprechend gliedert sich das Buch in 2 Hauptteile. Im 
ersten sind die allgemeinen Prinzipien, die beiden Hauptsatze, Zustandsgleichun 
gen, Mehrphasensysteme und noch kurz irreversible Prozesse und solche mit 
variabler Masse behandelt. Der zweite Teil ist den technisch so wichtigen statio- 
naren offenen Prozessen gewidmet. Darin wird, ausgehend von den Erhaltungs- 
satzen, die Theorie der stationar durchstromten Bilanzgebiete entwickelt und die 
gewonnenen allgemeinen Beziehungen auf die Stromung durch Diisen, durch 
Turbinen und Kolbenmaschinen, auf die Kaltetechnik und Gasverfliissigung 
angewendet. 

In Erfiillung der zweiten Aufgabe kampfen sich die Autoren wacker durch all 
jene Schwierigkeiten hindurch, die jeder auf sich ladt, der mit einem Ma8system 
mit dem Gewicht als Grundgrésse rechnet. Die Studenten, die zu dem hier durch- 
gehend verwendeten englisch-amerikanischen Ma®system abgerichtet werden, 
sind nicht zu beneiden! 

Der 3. und 4. Aufgabe dienen zahlreiche gut ausgewahlte und ausfiihrlich er- 
lauterte Ubungsbeispiele, aus deren Bearbeitung der Leser sicher grossen Nutzen 
ziehen wird. 

Lobend hervorzuheben ist die genaue Festlegung des Giiltigkeitsbereichs der 
Gleichungen. So bedeutet etwa ein in eckige Klammern hinter die Gleichung 
gesetztes P oder p, dass diese nur fiir P = const bzw. nur fiir das ideale Gas 
(perfect gas) gilt. Ein Anhang bringt ein Zustandsdiagramm des Methans (auch 
im Text wird meist auf Kohlenwasserstoffe und nicht auf den sonst so beliebten 
Wasserdampf zuriickgegriffen), eine tibersichtliche Darstellung der Ableitungen 
der wichtigsten Formeln, einige Tabellen und nochmals einige Ubungsaufgaben. 

Das didaktisch geschickt abgefasste Buch kann dem Ingenieur — weniger dem 
Chemie-ingenieur, denn es fehlt die Behandlung der Mehrkomponentensysteme 
gerne empfohlen werden, besonders wenn er sich in amerikanische Betrachtungs- 
weisen und Mafisysteme einarbeiten will. P. GRASSMANN 


Problémes de plasticité théorique. Von W. PRAGER (Dunod, Paris 1958) 
124 S., 54 Abb.; firs. 650.-. 

Es handelt sich um die Ubersetzung ins Franzésische des Werkes Probleme der 
Plastizitatstheorie, das hier [ZAMP 7, 471 (1956)] ausfiihrlich gewiirdigt worden 
ist. Gegeniiber der deutschen Ausgabe sind nur einige Stellen etwas verdeutlich 
und einige Literaturangaben beigefiigt worden, H. ZIEGLER 


t 


Torsionstheorie. Von C. WEBER und W. Ginruer (Friedr. Vieweg & Sohn, 
Braunschweig 1958). 307 S., 160 Abb.; DM 38.-. 
Das Buch behandelt die elastische Torsion prismatischer bzw. zylindrischer 
Stabe. C. WEBER hat dieses relativ beschrankte, aber dankbare und in mehr als 
einer Hinsicht interessante Teilgebiet der Elastizitatstheorie seit vielen Jahren 
machtig gefoérdert und darf heute wohl unbestritten als bester Kenner der Ma 


den Ingenieur von Interesse, sondern auch methodisch anregend und damit ancl 
fiir den Mathematiker sehr reizvoll ist. 


Der Beitrag des Verlegers zu diesem Buch ist dusserst bescheiden. 
H. ZIEGLER 
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Differential Equations Applied in Science and Engineering. Von 
HAROLD Waytanp (D. Van Nostrand Co., Princeton, N. J., 1957). 353 S.,66.ABb. 

Das vorliegende Buch von H. WayLanp stellt fiir Studierende der Physik 
und der Ingenieurwissenschaften eine yorziigliche Einfiihrung in das Gebiet der 
gewohnlichen und partiellen Differentialgleichungen dar. Der Inhalt gliedert sich 
in die folgenden Abschnitte: Vektorrechnung, gewéhnliche Differentialgleichun- 
gen, Losung durch Potenzreihen, spezielle Funktionen (BrssEL, LEGENDRE), 
Loésung durch Orthogonalfunktionen, Randwertprobleme, Hinweis auf Integral- 
transformationen. Jedes Kapitel enthalt eine knappe Zusammenfassung sowie 
Hinweise auf weiter ausholende Literatur; rund 330 Aufgaben dienen der Festi- 
gung und der Erweiterung des Stoffes. 

Entsprechend dem Charakter eines einfiihrenden Werkes stehen die linearen 
Probleme im Vordergrund, und im Falle der partiellen Differentialgleichungen 
sind es diejenigen, die sich durch einen Produktansatz lésen lassen. Der Verfasser 
legt grosses Gewicht auf die genaue Erfassung der physikalischen Probleme (die 
vorwiegend der klassischen mathematischen Physik und den Ingenieurwissen- 
schaften entstammen) und ihre Ubertragung in die Sprache der Mathematik. 
Anhangsweise folgt noch — um den Studierenden zu einem tiefern Eindringen an- 
zuleiten — eine Auswahl von besondern physikalischen Problemen, die in ausfiihr- 
licheren Werken oder in Fachzeitschriften behandelt sind. 

Es handelt sich um ein wohldurchdachtes und gut verstandlich geschriebenes 
Lehrbuch, das in ausgezeichneter Weise in die mathematische Behandlung 
physikalischer Probleme einfiihrt. E, RotH-DESMEULES 


. Subtabulation. A Companion Booklet to Interpolation and Allied 
Tables. Prepared by H. M. Nautical Almanac Office. Published by order of The 
Lords Commissioners of the Admiralty (Her Majesty’s Stationary Office, London 
2958). 54 S.; 7s. 6d. 

Dieses Biichlein ist eine Erganzung zu den Interpolationstabellen, die friher 
von denselben Verfassern publiziert wurden (Z. angew. Math. Phys. 97, 96). Es 
enthalt aile Hilfsmittel, um eine Funktion, die durch aquidistante Stiitzwerte 
‘gegeben ist, mit einem feineren Schritt zu tabellieren. Fir den Einzelrechner kann 
die Methode der «precalculated second differences» bestens empfohlen werden; die 
einfachen beigegebenen Hilfstabellen erméglichen eine sehr bequeme Subtabulation. 
Die im zweiten Teil geschilderten Methoden der «bridging differences» wurden ur- . 
spriinglich von L. J. Comrie entwickelt und eignen sich fiir die Verwendung von 
Buchungs- und Lochkartenmaschinen. Sie haben mit dem Aufkommen der elek- 
tronischen Automaten etwas an Bedeutung verloren, da auf einer programm- 
gesteuerten Maschine ein Interpolationsprogramm auch fiir die Subtabulation ein- 
gesetzt werden wird. E. STIEFEL 


_ H@here Technische Mechanik. Von IsrvAN SzaBo (Springer-Verlag, 
‘Berlin 1958). 498 S., 414 Abb.; DM 31 50. 
Der Verfasser des vorliegenden Lehrbuches, das eine Fortsetzung seiner Ein- 
fiihrung in die Technische Mechantk ist, gibt eine knappe Darstellung des grossen 
-Stoffgebietes, in der auch ein Grossteil der mathematischen Hilfsmittel entwickelt 
wird. 

Das erste Kapitel umfasst das Prinzip der virtuellen Arbeiten, das Prinzip von 
‘p’ALEMBERT, das Hamiltonsche Prinzip, die Lagrangeschen Gleichungen, einen 


kurzen Einblick in die Variationsrechnung, das Ritzsche Verfahren und Anwendungen 
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auf die Statik und Dynamik starrer und elastischer Korper. Im zweiten Kapitel 
wird der ebene und achsensymmetrische Spannungszustand, die Theorie der diinnen 
Platten, eine Einfithrung in die Schalentheorie und die Torsion von Staben und 
Balken behandelt. Das dritte Kapitel bringt Plastizitatstheorie und das vierte die 
Dynamik der idealen Fliissigkeiten und Gase und Ausfiihrungen liber zahe Flissig- 
keiten. An vielen Beispielen und 73 Aufgaben mit Lésungen wird die Anwendung 
der dargestellten Theorien gezeigt. Bei den meisten Problemen, die keine exakte 
Lésung zulassen, wird eine Naherungslésung angegeben oder der Weg zu einer 
solchen angedeutet. Ein Anhang mit 11 vermischten Aufgaben aus verschiedenen 
Gebieten des behandelten Stoffes dient teilweise als Erganzung der Theorie oder 
enthalt Lésungen von praktisch wichtigen Fragen. 

Mit diesem Buch wird zur Entwicklung beigetragen, in der sich die theoretische 
Mechanik und die technische Mechanik angleichen. Aus diesem Grund wird es sich 
sowohl fiir den Mathematiker und Physiker als auch fiir den mathematisch inter 
essierten Ingenieur als wertvoll erweisen. Cu. WEHRLI 


Stabilitatstheorie mit Erlauterungen zu DIN 4114. Von G. BURGER-— 
MEISTER unter Mitwirkung von H. Steup, Teil I (Akademie-Verlag, Berlin 1957). 
407 S., 306 Abb.; DM 35.50. 

Dieses Buch, das der Theorie und Praxis der elastischen und inelastischen Stabi 
litatsprobleme gewidmet ist, bringt einen erfreulichen Fortschritt auf diesem, in 
seinen Grundlagen noch immer nicht sehr fest gefiigten Gebiet. Zwar wird auch hier 
noch zwischen «statischer» und «dynamischer» Stabilitat unterschieden; die ver- 
schiedenen Spielarten, in denen die Stabilitatsprobleme auftreten, werden indessen 
sorgfaltig gegeneinander abgegrenzt, in ihren typischen Merkmalen untersucht und 
an Beispielen besprochen. Besondere Vorziige des Buches liegen darin, dass gelegent- 
lich auf die Grenzen der gebrauchlichen Verfahren hingewiesen wird und dass die 
iiblichen Vernachlassigungen ins rechte Licht geriickt und Schritt fiir Schritt_ 
behoben werden. Ein grosser Teil des Textes ist den Vorgaéngen nach Eintritt des 
Knickens, dem inelastischen Verhalten und der Ermittlung von Traglasten ge- 
widmet. 

Dieser erste Band (mit den Haupttiteln: Einfiihrende Betrachtungen zur Sta- 
bilitat des Gleichgewichts, Lésung von Stabilitatsproblemen nach der Gleichgewichts- 
methode, Lésung von Stabilitatsproblemen nach der energetischen Methode, Nahe- 
rungsmethoden zur Berechnung von Verzweigungslasten, Mehrteilige Druckstabe 
Biegedrillknickung und. Kippung) enthalt alle Probleme, die bei Staben und Balken 
auftreten, und ist gut dokumentiert. Ein zweiter Band soll sich mit Stabsystemen 
und Flachentragwerken befassen. H. ZIEGLER 


Grenzschicht-Theorie. Von H. Scuiicutine, 3. Aufl. (G. Braun, Karlsruhe, 
1953)2003!5., 375 Abbs; 32) Labs chra72-40) 


schaft angepasste Ausgabe vor. Wenn auch die Einteilung im wesentlichen unver- 
andert geblieben ist, sind jedoch wichtige Erganzungen zugefiigt worden. Insbeson 
dere das Kapitel iiber die kompressible laminare Grenzschicht und der Abschnitt 
Umschlag laminar-turbulent sind erheblich erweitert worden, einem Wunsch der 
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ersten Rezension Rechnung tragend. Damit ist der Unifang des Buches von 500 auf 
620 Seiten gestiegen. Ein wesentlicher Vorteil des Werkes ist der umfangreiche, die 
neuesten Arbeiten beriicksichtigende Literaturnachweis. 

Das Werk stellt eine ausgezeichnete Eimfithrung in die physikalischen und mathe- 
matischen Grundlagen der Stroémung reibungsbehafteter Flissigkeiten dar. Die zahl- 
reichen Vergleiche mit experimentellen Ergebnissen erlauben eine genaue Abgren- 
zung des Giiltigkeitsbereiches der Theorie und geben besseren Einblick in die physi- 
kalischen Vorgange. 
|} Fir den Studierenden sowie fiir den praktisch in Industrie oder Forschung 
tatigen Ingenieur wird diese Grenzschicht-Theorie lange ein Standardwerk bleiben. 

P. DE HALLER 


Die Wirbelschicht als Energieiibertragungsfliche. Von W. Marz 
(Springer-Verlag, Berlin 1958). 66 S., 27 Abb.; DM 10.50. 

Der Verfasser versucht, die Vorginge bei dem Energieaustausch durch eine 
vektoranalytische Betrachtungsweise zu beschreiben. Dadurch sollten allgemeine 
Gesetze ermittelt werden, die fiir mechanische, strémungsdynamische oder thermo- 
dynamische Energieiibertragung Geltung haben. Den Ausfihrungen ist ein 
gewisser formaler Reiz nicht abzusprechen, der praktische Wert dieser Theorie aber 
scheint dem Referenten etwas fraglich. Die Behauptung, dass mit dieser Betrach- 
tungsweise die Verlustquellen eines Verfahrens leichter zu erkennen seien, bleibt 
unbewiesen, beziehen sich doch alle aufgefiihrten Beispiele auf langst bekannte Vor- 
gange, deren Behandlung mit den klassischen Mitteln der Mechanik und der 
Thermodynamik ersché6pfend méglich ist. Immerhin tragen solche Verallgemeine- 
rungsversuche zum Verstandnis verwandter Probleme bei und sind in diesem Sinne 
zu begriissen. P. DE HALLER 


' Gas Dynamics. Von Kraus OswatitscH. Englisch von Gustav KUERTI. 


Band 1 von Applied Mathematics and Mechanics. Herausgegeben von F. N. FREN- 
KIEL (Academic Press, New York 1956). 610 S., 300 Abb.; $ 12.00. 

' Dies ist die englische Ubersetzung der in dieser Zeitschrift bereits ausfiihrlich 
besprochenen deutschen Ausgabe. Sie zeichnet sich durch eine klare, einfache 
Sprache aus sowie durch ein Minimum an mathematischen Exkursen. Es ist den 
Autoren gelungen, eine leichtverstandliche und logisch aufgebaute Einfiihrung in 
die Methoden der modernen Aerodynamik der hohen Geschwindigkeiten zu geben. 
Das Buch behandelt nach einem kurzen Uberblick der notwendigen thermodyna- 
mischen Begriffe zunachst die eindimensionale stationare und instationare Stré- 
mung, um nachher allgemeine Gleichungen, Integralsatze und Ahnlichkeitsregeln 
abzuleiten. Zwei- und dreidimensionale Str6mungen in Schallnahe, im Uber- und 
Unterschallgebiet werden ausfiihrlich behandelt, rechnerische und experimentelle 
Hilfsmittel sowie der Einfluss der Viskositat bilden je den Inhalt eines Kapitels. 
Wertvoll sind die numerischen Tafeln sowie das ausfiihrliche, die neuesten Erschei- 
nungen beriicksichtigende Literaturverzeichnis. P. DE HALLER 


_ Experimental Designs. Von W. G. Cocuran und G. M. Cox, 2. Aufl. (John 
Wiley and Sons, New York 1957). 617 S.; $ 10.25. . 

Die erste Auflage dieses Werkes erschien 1950 und wurde bald als ein Standard- 
werk anerkannt. Die Grundsatze fiir das richtige Planen und Auswerten von Ver- 
suchen gehen auf R. A. FisHer zuriick, der 1935 in seinem klassischen Werk The 
Design of Experiments die grundlegenden Ideen vortrug. CocHRAN und Cox kommt 
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das Verdienst zu, die bis 1950 entwickelten Versuchsplane tibersichtlich dargestellt 
und auf Grund zahlreicher Beispiele deren Auswertung veranschaulicht zu haben. 

Wie die gesamte mathematische Statistik, so ist auch das Teilgebiet des Planens— 
von Versuchen in kraftiger Entwicklung begriffen. Seit 1950 sind weitere Plane von 
bekannten Systemen neu gefunden und ihre Auswertung untersucht worden. Dies 
gilt besonders fiir Versuche mit mehreren Faktoren, bei denen nur ein Teil der még- 
lichen Kombinationen ausgefiihrt wird («fractional factorial designs»). COCHRAN = 
und Cox haben iiber diesen Gegenstand ein Kapitel beigefigt. 

Der zweite Gegenstand, dem die Verfasser ein neues Kapitel widmen, sind die 
von G. E. P. Box entwickelten Verfahren zum Auffinden der optimalen Kombination 
mehrerer Faktoren. : 

Das Buch von Cocuran und Cox wird in seiner neuen Form noch in vermehrtem 
Masse ein niitzliches Werkzeug in der Hand des Versuchsanstellers sein. 

A. LINDER™ 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLUGGE. Band 5, 1. Teil: Prinzipien der Quantentheorie (Springer-Verlag, Ber- 
lin 1958). 376 S., 7 Fig.; DM 90.-. 

Der 1. Teil des 5. Bandes behandelt die Prinzipien der Quantentheorie und 
enthalt einen Beitrag von W. Pauti tiber Die allgemeinen Prinzipien der Wellen- 
mechanik und einen zweiten von G. KALLEN tiber Quantenelektvodynamtk. 

Beim Artikel von W. Pautt handelt es sich im wesentlichen um einen Neu- 
druck des klassischen Artikels aus dem Band 24, Teil 1, des Handbuches von 
GEIGER und ScHEEL. Die Anderungen betreffen zur Hauptsache die relativistische 
Quantentheorie, wo der Abschnitt iiber das elektromagnetische Feld weggelassen 
wurde. An seine Stelle tritt der ausfiihrliche Artikel von G. KALLEN. (Die erste 
Fussnote des Bandes verrat iibrigens den urspriinglichen Plan des Herausgebers, 
der nur teilweise zur Ausfiihrung gelangen konnte.) Ebenso entfallen natiirlich die 
beritihmte Kritik an der Diracschen Theorie des Elektrons und einige Worte in 
der Diskussion der Weylschen zweikomponentigen Gleichung fiir ein Spin 1/2- 
Teilchen der Ruhmasse Null. Der Rezensent vermisst besonders die erste Stelle 
mit einiger Wehmut, denn sie illustriert, welch iiberraschende Lésungen in 
scheinbar hoffnungsloser Situation die Natur hin und wieder bereithalt. 

Im iibrigen braucht der Neudruck des Paulischen Artikels kein Wort der 
Rechtfertigung. Er ist und bleibt die meisterhafteste Behandlung des Gegen- 
standes. 

Der Artikel von G. KALLEN ist weitgehend aufs Sachlich-Technische ausge- 
richtet. Grundsatzliche Betrachtungen iiber die Messbarkeit des elektromagne- 
tischen Feldes werden nur gestreift. Andererseits wird auf alle Matzchen ver 
zichtet (was in der Feldtheorie nicht ganz selbstverstandlich ist). Mir erscheint 
diese Einstellung gerechtfertigt, handelt es sich doch um eine Theorie, deren 
Grundlagen weitgehend schwankend sind, deren Erfolge aber zum Bedeutendsten 
gehéren, was die Physik der letzten 10 Jahre hervorgebracht hat. Mit Recht 
werden daher die Anwendungen in gebiihrender Breite behandelt. Das gilt auch 
fiir die niedrigsten Naherungen fiir den Compton-Effekt, die Paarerzeugung, die 
Bremsstrahlung und die Moller-Streuung. 

Den Kern des Artikels bildet der VI. Abschnitt, in welchem die niedrigsten 
strahlungstheoretischen Korrekturen behandelt werden. Hier findet man neben 
den «klassischen» Effekten (magnetisches Moment des Elektrons, Lamb-Verschie- 
bung) auch einen Abschnitt iiber Positronium. Die Rechnungen fiir die hdheren 
Korrekturen kénnen natiirlich nur erwahnt werden. 
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Der VII. Abschnitt behandelt die allgemeine Theorie der Renormierung unab- 
hangig von der Stdrungsrechnung und gibt zum erstenmal eine buchmassige Dar- 
stellung des Kallénschen Ergebnisses, dass unter den Renormierungskonstanten 
wenigstens eine nicht endlich ist. 

Der Kallénsche Artikel diirfte auch als Einfiihrung in die Feldtheorie und die 
Quantenelektrodynamik eine willkommene Erganzung des Schrifttums dar- 
stellen. ; R. Jost 


| Nichtlineare Mechanik. Von H. Kauprrer (Springer-Verlag, Berlin 1958). 
684 S., 229 Abb.; DM 64.50. 

Genau genommen, sind die Differentialgleichungen, welche die in der Natur und 
insbesondere in der Mechanik auftretenden Vorgange beschreiben, nichtlinear. Auf 
den meisten Gebieten (Mechanik der Fliissigkeiten und Gase, der plastischen Korper 
und des Kreisels) gelingt ihre Linearisierung nur fiir ganz spezielle Falle. Im Gegen- 
satz hiezu ist es in der Schwingungslehre und in der Elastizitatstheorie méglich, die 
Differentialgleichungen im vornherein und ohne Riicksicht auf das spezielle Pro- 
blem zu linearisieren, und diese beiden Gebiete sind denn auch bis heute die ergiebig- 
sten Jagdgriinde des mechanischen Waidwerks geblieben. Vorstésse ins nichtlineare 
Gebiet sind nicht nur miihsam, sondern auch weniger ergiebig, da hier das Uberlage- 
rungsprinzip nicht gilt. Sie sind aber sehr interessant, da die Linearisierung vielfach 
nicht nur die Genauigkeit der Resultate affiziert, sondern feinere Erscheinungen 
uberhaupt unter den Tisch wischt. Seit den zwanziger Jahren sind bei den nicht- 
linearen Schwingungen grosse Fortschritte erzielt worden, wobei sich KAUDERER 
bleibende Verdienste erworben hat. In der Elastizitat haben MURNAGHAN, RIVLIN 
und andere eine geometrisch nichtlineare Theorie geschaffen; gleichzeitig ist von 
KAUDERER und seinen Mitarbeitern eine physikalisch nichtlineare Theorie ent- 
wickelt worden. 

, Das vorliegende Buch berichtet iiber den heutigen Stand der Forschung auf den 
beiden Gebieten, auf denen der Verfasser tatig ist: die nichtlinearen Schwingungen 
und die von ihm geschaffene physikalisch nichtlineare Elastizitat. Zahlreiche Ergeb- 
nisse sind in aller Stille herangereift und hier zum ersten Mal veréffentlicht. 

Der Stil des Buches ist vorbildlich, die Darstellung leicht verstandlich und von 
einer seltenen Prazision. H. ZIEGLER 


Physikalische Abhandlungen und Vortrage. Von Max PLanck (Friedr. 

Vieweg & Sohn, Braunschweig 1958). 3 Bande; DM 150.-. 
_ Mit der vorliegenden dreibandigen Ausgabe der Physikalischen Abhandlungen 
und Vortvdge von MAx PLANCK erfiillt der Verband der Deutschen Physikalischen 
Gesellschaften und die Max-Planck-Gesellschaft zur Forderung der Wissenschaften 

ine moralische Verpflichtung, die seit PLancKs Tod im Jahre 1947 bestanden hat. 
Teh denke, dass jeder Physiker den erwahnten Gesellschaften fiir ihre Tat dankbar 
ist. 
' Die beiden ersten Bande enthalten zusammen 121 wissenschaftliche Abhand- 
lungen aus den Jahren 1879 bis 1941. Weggelassen wurden einige wenige rein pole- 
mische Arbeiten und die Lehrbiicher. Letztere sind ja sonst leicht zuganglich. 
a Der 3. Band enthailt Reden und Aufsitze, die sich an einen breiteren Kreis 
wenden, dann Jubilaumsadressen und Nachrufe auf zeitgendssische Physiker, und 
= aufschlussreiche wissenschaftliche Selbstbiographie. Schliesslich Reden 
und Erwiderungen zu Max PLancks 80. Geburtstag, die Traueransprache von Max 
von LAvE und zwei Rundfunkansprachen von Orro Hann und Max von LavE 


zm 10. Todestag. 


DR NE Rs RE 


r 


112 | Buchbesprechungen — Book Reviews — Notices bibliographiques 


“Der nicht vorwiegend an Thermodynamik interessierte Leser wird sich besonders 
fiir PLancxs Arbeiten zur Hohlraumstrahlung und zur Quantentheorie interessie- 
ren. Dazu gehért jene Kette von Arbeiten, die zu den beiden umstiirzenden Mit- 
teilungen vom 19. Oktober und 14. Dezember 1900 vor der Physikalischen Gesell- 
schaft zu Berlin fiihrten. Was an diesen Mitteilungen, neben ihrer historischen Be- 
deutung, auch heute noch ergreift, ist ihre enorme Unmittelbarkeit. Der Schluss 
der zweiten Mitteilung (in welcher die Quantenhypothese zur Herleitung der Inten= 
sitatsverteilung der Hohlraumstrahlung zuerst eingefiihrt wird) gibt zum ersten 
Mal zuverlassige Gréssen fiir die Loschmidtsche Zahl und die Elektronenladung und 
schliesst mit der stolzen Bemerkung: : ; 

«All diese Beziehungen beanspruchen, wenn die Theorie tiberhaupt richtig ist, 
nicht annadhernde, sondern absolute Giiltigkeit... Ihre Priifung durch direktere 
Methoden wird eine ebenso wichtige wie notwendige Aufgabe der weiteren For 
schung sein.» = 

Nicht minder interessant aber ist PLancxs Einstellung zur weiteren Entwick- 
lung der Quantentheorie. Sie ist gekennzeichnet durch den fast verzweifelt anmu: 
tenden Versuch, an der klassischen Physik festzuhalten. on had 

Das folgende Zitat aus dem Jahre 1923 zur Bohrschen Atomtheorie mag hier 
zur Illustration dienen: . a 

«Zunachst muss es jedem Unbefangenen kaum begreiflich vorkommen, wie es 
iiberhaupt méglich war, dass ein derartiger Frontalangriff auf die so sicher fun- 
dierte, so sorgfaltig abgerundete und vielseitig erprobte klassische Theorie ... so 
rasch gelingen konnte, und es mag im Zusammenhang damit wohl sogar ein Zweifel 
sich regen, ob angesichts solcher Vorkommnisse die Physik iiberhaupt noch den 
Anspruch erheben darf, als die bestbegriindete unter den Naturwissenschaften zu 
gelten.» 

Dieser Zweifel wird zwar im folgenden einigermassen besanftigt, beherrscht aber 
doch PLancKs weitere Arbeiten zur Quantentheorie, welche zur Hauptsache aus 
einer Reihe von «klassizistischen» (und erfolglosen) Modifikationen der Quanten- 
hypothese und spater der Wellenmechanik bestehen. - fd 

Es ist tréstlich, aus der wissenschaftlichen Selbstbiographie zu héren, dass 
PLANCK in diesen vergeblichen, miihsamen Versuchen nichts Tragisches findet, 
sondern dass er ihretwegen schliesslich genau wusste, «dass das Wirkungsquantum 
in der Physik eine viel bedeutendere Rolle spielt, als ich anfangs geneigt war a 
zunehmen ». 7 

Die Reproduktion der Abhandlungen geschah zum gréssten Teil photomecha: 
nisch. Sie ist im allgemeinen adaquat. ; ; 

R. Jost 

= 

Microwave Transmission Design Data. Von THEODORE MoRENO (Dover 
Publications, New York 1958). 248 S., 200 Fig.; $1.50 (2. Auflage des 1948 
Sperry Gyroscope Company erschienenen Buches). » Se 

Der Fachmann wird es freudig begriissen, dass MorENos Microwave Transmis- 


-nenen ersten Auflage; dass keine Anderungen des Inhaltes not i urd: 
spricht fiir die Qualitat des Buches. aa ats “ 
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